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review recent developments in tlie pliysics of ultracold atomic and molecular gases in optical 
'-T^attices. Such systems are nearly perfect realisations of various kinds of Hubbard models, and as 
^J^ilch may very well serve to mimic condensed matter phenomena. We show how these systems 
^^rriay be employed as quantum simulators to answer some challenging open questions of condensed 
^^^taiatter, and even high energy physics. After a short presentation of the models and the methods of 
C _ t reatment of such systems, we discuss in detail, which challenges of condensed matter physics can 
"j-'^e addressed with (i) disordered ultracold lattice gases, (ii) frustrated ultracold gases, (iii) spinor 
j^attice gases, (iv) lattice gases in "artificial" magnetic fields, and, last but not least, (v) quantum 

S'nformation processing in lattice gases. For completeness, also some recent progress related to the 
ubove topics with trapped cold gases will be discussed. 
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Motto: 

There are more things in heaven and earth, Horatio, 
Than are dreamt of in your philosophy [T] . 



1 Introduction 

1.1 Cold atoms from a historical perspective 

Thirty years ago, atomic physics was a very well established and respectful, but ev- 
idently not a "hot" area of physics. On the theory side, even though one had to 
deal with complex problems of many electron systems, most of the methods and 
techniques were developed. The main questions concerned, how to optimize these 
methods, how to calculate more efficiently, etc. These questions were reflecting an 
evolutionary progress, rather than a revolutionary search for totally new phenomena. 
Quantum optics at this time was entering its Golden Age, but in the first place on the 
experimental side. Development of laser physics and nonlinear optics led in 1981 to 
the Nobel prize for A.L. Schawlow and N. Blocmbergen "for their contribution to the 
development of laser spectroscopy" . Studies of quantum systems at the single particle 
level culminated in 1989 with the Nobel prize for H.G. Dehmelt and W. Paul "for the 
development of the ion trap technique" , shared with N.F. Ramsey "for the invention 
of the separated oscillatory fields method and its use in the hydrogen maser and other 
atomic clocks" . 

Theoretical quantum optics was born in the 60ties with the works on quantum 
coherence theory by the 2005 Noble prize winner, R.J. Glauber [HIS], and with the 
development of the laser theory by M. Scully and W.E. Lamb (Nobel laureate of 
1955) and H. Haken [5]. In the 70ties and 80ties, however, theoretical quantum 
optics was not considered to be a separate, established area of theoretical physics. One 
of the reasons of this state of art, was that indeed the quantum optics of that time 
was primarily dealing with single particle problems. Most of the many body problems 
of quantum optics, such as laser theory, or more generally optical instabilities [S], 
could have been solved either using linear models, or employing relatively simple ver- 
sions of mean field approach. Perhaps the most sophisticated theoretical contributions 
concerned understanding of quantum fluctuations and quantum noise [Sl[7]. 

This situation has drastically changed in the last ten - fifteen years, and there are 
several seminal discoveries that have triggered these changes: 

• First of all, atomic physics and quantum optics have developed over the years 
quite generally an unprecedented level of quantum engineering, i.e. preparation, 
manipulation, control and detection of quantum systems. 

• Cooling and trapping methods of atoms, ions and molecules have reached regimes of 
low temperatures (today down to nanoKelvin!) and precision, that 15 years ago were 
considered unbelievable. These developments have been recognised by the Nobel 
Foundation in 1997, who awarded the Prize to S. Chu [5], C. Cohen- Tannoudji [5] 
and W.D. Phillips [Tn] "for the development of methods to cool and trap atoms 
with laser light" . Laser cooling and mechanical manipulations of particles with 
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light [TT] was essential for development of completely new areas of atomic physics 
and quantum optics, such as atom optics [12j . and for reaching new territories of 
precision metrology and quantum engineering. 

• Laser cooling combined with evaporation cooling technique allowed in 1995 for 
experimental observation of Bose-Einstein condensation (BEC) [T31II1], a phe- 
nomenon predicted by S. Bose and A. Einstein more than 70 years earlier. The 
authors of these experiments, E.A. Cornell and C.E. Wieman [15] and W. Ket- 
terle [T^ received the Noble Prize in 2001, "for the achievement of Bose-Einstein 
condensation in dilute gases of alkali atoms, and for early fundamental studies of 
the properties of the condensates" . This was a breakthrough moment, in which 
"atomic physics and quantum optics has met condensed matter physics" [H] . Con- 
densed matter community at this time remained, however, still reserved. After all, 
BEC was observed in weakly interacting dilute gases, where it is very well described 
by the mean field Bogoliubov-de Gennes theory [TJ • Although the finite size of the 
systems, and spatial inhomogeneity play there a crucial role, the basic theory of 
such systems was developed in the 50'ties. 

• The seminal theoretical works of the late A. Peres [T5], the proposals of quan- 
tum cryptography by C.H. Bennett and G. Brassard [H], and A.K. Ekert [5D|, the 
quantum communication proposals by C.H. Bennett and S.J. Wiesner |21] and C.H. 
Bennett, G. Brassard, C. Crepeau, R. Jozsa, A. Peres, and W.K. Wootters [52], the 
discovery of the quantum factorizing algorithm by P. Shor [23j . and the quantum 
computer proposal by J.I. Cirac and P. ZoUer [24j have given birth to experimental 
studies of quantum information [25ll26j . These studies, together with rapid devel- 
opment of the theory have, in particular in the area of atomic physics and quantum 
optics, led to enormous progress in our understanding of what quantum correlations 
and quantum entanglement are, and how to prepare and use entangled states as a 
resource. The impulses from quantum information enter nowadays constantly into 
the physics of cold atoms, molecules and ions, and stimulate new approaches. It is 
very probable that the first quantum computers will be, as suggested already by 
Feynman [27j , computers of special purpose - quantum simulators |28] , that will ef- 
ficiently simulate quantum many body systems that otherwise cannot be simulated 
using "classical" computers [25] ■ 

• The physics of ultracold atoms entered the areas of strongly correlated systems 
with the seminal 1998 paper of Jaksch et al. [5S] on the superfluid-Mott insulator 
transition in cold atoms in an optical lattice. Apart from stimulation from the con- 
densed matter physics |30| , the authors of this paper were in fact motivated by the 
possibility of realising quantum computing with cold atoms in a lattice |31j . Transi- 
tion to the Mott insulator state was supposed to be an efficient way of preparation 
of a quantum register with a fixed number of atoms per lattice site. The experi- 
mental observation of the superfluid-Mott insulator transition by the Bloch-Hansch 
group [32] (see Fig. [1]) marks the beginning of age of the experimental studies of 
strongly correlated systems with ultracold atoms [33]. Several other groups have 
observed bosonic superfluid-Mott insulator transitions in pure Bose systems [31] . 
in disordered Bose systems [35], or in Bose- Fermi mixtures [36l|37]. Very recently 
MI state of molecules have been created [35] , and bound repulsive pairs of atoms 
(i.e. pairs of atoms at a site that cannot release their repulsive energy due to the 
band structure of the spectrum in the lattice) have been observed |39| . 

Since 1998 the physics of ultracold atoms has made enormous progress in the stud- 
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Figure 1. Distribution of atoms among lattice sites vs. momentum distribution of atoms released 

from an optical lattice, (a): low lattice potential: the largest interference peak is at zero 
momentum. The system is a superfluid with fluctuating number of atoms per site, (b): high lattice 
potential case where a perfect Mott state with exactly one atom per site is achieved - the 
momentum distribution is blurred (from I33I V 

ies of strongly correlated systems. Number of the theory papers that propose to mimic 
various condensed matter systems of interest is hardly possible to follow, and the num- 
ber of experiments in which strongly correlated regime has been met grows also very 
significantly. Moreover, condensed matter physicists, in particular theorists, joint the 
efforts of atomic physicists and quantum opticians. Among those who have "commit- 
ted" a paper on cold atoms are Noble Prize winners: A.J. Leggett F. Wilczek |¥T| . 
D. Politzer [32] , or such personalities of condensed matter theory as M.P. A. Fisher [33] , 
or mathematical physics as M. Aizenman or E.H. Lieb [44] . 



1.2 Cold atoms and the challenges of condensed matter physics 

The physics of cold atoms touches nowadays the same frontiers of modern physics 
as condensed matter and high energy physics. In particular many of the important 
challenges of the latter two disciplines can be addressed in the context of cold atoms: 

• ID systems. The role of quantum fluctuations and correlations is particularly im- 
portant in ID. The theory of ID systems is very well developed due to existence 
of exact methods such as Bethe ansatz and quantum inverse scattering theory 
(c.f. [3S]), powerful approximate approaches, such as bosonisation. or conformal 
field theory [46j . and efficient computational methods, such as density matrix renor- 
malisation group technique (DMRG. c.f. [3Z])- There are, however, many open 
experimental challenges that have not been so far directly and clearly realized in 
condensed matter, and can be addressed with cold atoms (for a review see [IH])- Ex- 
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amples include atomic Fermi, or Bose analogues of spin-charge separation, or more 
generally observations of microscopic properties of Luttinger liquids [^^[501151] . Re- 
cent experimental observations of the ID gas in the deep Tonks-Girardeau regime 
by Paredes et al. [S2] (see also ) are the first steps in this direction. 

• Spin-boson model. A two state system coupled to a bosonic reservoir is a paradigm, 
both in quantum optics, as well as in condensed matter physics, where it is termed 
as spin-boson model (for a review see [ST])- It has also been proposed [55] that an 
atomic quantum dot, i.e., a single atom in a tight optical trap coupled to a superfluid 
reservoir via laser transitions, may realise this model. In particular, atomic quantum 
dots embedded in a ID Luttinger liquid of cold bosonic atoms realizes a spin-boson 
model with Ohmic coupling, which exhibits a dissipative phase transition and allows 
to directly measure atomic Luttinger parameters. 

• 2D systems. According to the Mermin-Wagner-Hohenberg theorem, 2D systems 
with continuous symmetry do not exhibit long range order at finite temperatures 
T. 2D systems may, however, undergo Kosterlitz-Thouless-Berezinskii transition 
(KTB) to a state in which correlations decay is algebraic, rather than exponen- 
tial. Although KTB transition has been observed in liquid Helium, its microscopic 
nature (binding of vortex pairs) has never been seen. A recent experiment of Dal- 
ibard's group [55] makes an important step in this direction. 

• Hubbard and spin models. Very many important examples of strongly correlated 
state in condensed matter physics arc realised in various types of Hubbard models 
[45ll60] . While in condensed matter Hubbard models are "reasonable caricatures" of 
real systems, ultracold atomic gases in optical lattices allow to achieve practically 
perfect realizations of a whole variety of Hubbard models [61j . Similarly, in certain 
limits Hubbard models reduce to various spin models; again cold atoms and ions 
allow for practically perfect realizations of such spin models (see for instance |65[ 
[64^.63.62] ). Moreover, one can use such realizations as quantum simulators to mimic 
specific condensed matter models, and to address various very particular and well 
focused questions. 

• Disordered systems: Interplay localisation-interactions. Disorder plays a central role 
in condensed matter physics, and its presence leads to various novel types of effects 
and phenomena. One of most prominent quantum signatures of disorder is Ander- 
son localisation [66j of the wave function of single particles in a random potential. 
The question of the interplay between disorder and interactions has been inten- 
sively studied. For attractive interactions, disorder might destroy the possibility of 
superfluid transition ("dirty" superconductors). Weak repulsive interactions play 
a delocalising role, whereas very strong ones lead to Mott type localisation [67] . 
and insulating behaviour. In the intermediate situations there exist a possibility of 
delocalized "metallic" phases. Cold atom physics starts to investigate these ques- 
tions. Controlled disorder, or pseudo-disorder might be created in atomic traps, 
or optical lattices by adding an optical potential created by speckle radiation, or 
several lattices with incommensurate periods of spatial oscillations [SSKnSKTD] . In 
an optical lattice this should allow to study Anderson-Bose glass and crossover 
between Anderson-like (Anderson glass) to Mott type (Bose glass) localisation. 
Very recently, Bose glass state has been realized experimentally by the M. Inguscio 
group [35] . The same group, as well as several others have initiated experimen- 
tal and theoretical studies of the role of interactions in Anderson localisation ef- 
fects for trapped Bose gases [7T1[7^[711[75] . According to theoretical predictions of 



Ref. [79 l[73ll77l l78 l[76ll75] , prospects of detecting signatures of Anderson localisation 
in the presence of weak nonlinear interactions and quasi-disorder in BEC arc quite 
promising. One expects in such system the appearance of a novel Lifshits glass 
phase [5D] , where bosons condense in a finite number of states from the low energy 
tail (Lifshits tail) of the single particle spectrum. 

Disordered systems: spin glasses. Since the seminal papers of Edwards and Ander- 
son [5T], and Sherrington and Kirkpatrick [55] the question about the nature of 
the spin glass ordering has attracted a lot of attention [531|H31[55] ■ The two com- 
peting pictures: the replica symmetry breaking picture of G. Parisi, and a droplet 
model of D.S. Fisher and D.A. Huse are probably applicable in some situations, 
and not applicable in others. Cold atom physics might contribute to resolving 
this controversy, and even add understanding of some quantum aspects, like for 
instance behaviour of Ising spin glasses in transverse (i.e. quantum mechanically 
non-commuting) fields [M[|87j . 

Disordered systems: Large effects by small disorder. There are many examples of 
such situations. In classical statistical physics a paradigm is the random field Ising 
model in 2D (that looses spontaneous magnetization at arbitrarily small disorder). 
In quantum physics the paradigmatic example is Anderson localisation, which oc- 
curs at arbitrarily small disorder in ID, and should occur also at arbitrarily small 
disorder in 2D. Cold atom physics may address these questions, and, in fact, much 
more (cf. the Ref. [88| , where a disorder breaks the continuous symmetry in a spin 
systems, and thus allows for long range ordering). 

High Tc superconductivity. Despite many years of research, opinions on the nature 
of high Tc superconductivity still vary quite appreciably [89| . It is, however, quite 
established (cf. contribution of P.W. Anderson in Ref. that understanding of 
the 2D Hubbard model in the, so called, J limit [501 fo^' two component (spin 1/2) 
fermions provides at least a part of the explanation. The simulation of these models 
are very hard and numerical results are also full of contradictions. Cold fermionic 
atoms with spin (or pseudospin) 1/2 in optical lattices might provide a quantum 
simulator to resolve these problems [90] (see also [91]). First experiments with both 
"spinless", i.e. polarized, as well as spin 1/2 unpolarized ultracold fermions [Mll^ . 
and Fermi-Bose mixtures in lattices has already been realized [5Sll37j . 

BCS-BEC cross-over. Physics of high Tc superconductivity can be also addressed 
with trapped ultracold gases. Weakly attracting spin 1/2 fermions in such situ- 
ations undergo at (very) low temperatures the Bardecn-Cooper-Schricffer (BCS) 
transition to a supcrfluid state of loosely bounded Cooper pairs. Weakly repuls- 
ing fermions, on the other hand may form bosonic molecules, which in turn 
may form at very low temperatures a BEC. Strongly interacting fermions un- 
dergo also a transition to the superfluid state, but at much higher T. Sev- 
eral groups have employed the technique of Feshbach resonances [551 [Ml 
to observe such BCS-BEC cross-over (for the recent status of experiments see 

[MiiHTiinaiMiiinaiiniiiinaiiMiiMiinsiiins] ) . 

Frustrated antiferromagnets and spin liquids. The "rule of thumb" says that every- 
where, in a vicinity of a high Tc superconducting phase, there exists a (frustrated) 
antiferromagnetic phase. Frustrated antiferromagnets have been thus in the centre 
of interest in condensed matter physics for decades. Particularly challenging here 
is the possibility of creating novel, exotic quantum phases, such as valence bond 
solids, resonating valence bond states, and various kinds of quantum spin liquids 
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(spin liquids of I and II kind, according to C. Lhuillier [107|,ll08j . and topological 
and critical spin liquids, according to M.P.A. Fisher [109p . Cold atoms offer also 
in this respect opportunities to create various frustrated spin models in triangular, 
or even kagome lattice In the latter case, it has been proposed by Damski et 
al. |110lllllj that cold dipolar Fermi gases, or Fermi-Bose mixtures might allow to 
realize a novel state of quantum matter: quantum spin liquid crystal, characterized 
by Neel like order at low T (see also |112j ). accompanied by extravagantly high, 
liquid- like density of low energy excited states. 

• Topological order and quantum computation. Several very "exotic" spin systems 
with topological order has been proposed recently [113llll4j as candidates for robust 
quantum computing. Despite their unusual form, these models can be realized with 
cold atoms [Mill 15] . Particularly interesting [116j is the recent proposal by Micheli 
et al. [115| . who propose to use heteronuclear polar molecules in a lattice, excite 
them using microwaves to the lowest rotational level, and employ strong dipole- 
dipole interactions in the resulting spin model. The method provides an universal 
"toolbox" for spin models with designable range and spatial anisotropy of couplings. 

• Systems with higher spins. Lattice Hubbard models, or spin systems with higher 
spins are also related to many open challenges; perhaps the most famous be- 
ing the Haldane conjecture concerning existence of a gap, or its lack for the ID 
antiferromagnetic spin chains with integer or half-integer spins, respectively (see 
for instance [BD|). Ultracold spinor gases |117| might help to study these ques- 
tions. Again, particularly interesting are in this context spinor gases in optical 
lattices [TT8l[TT9l[T20ill21[|122j . where in the strongly interacting limit the Hamilto- 
nian reduces to a generalized Heisenbcrg Hamiltonian. Using Fcshbach resonances 
and varying lattice geometry one should be able in such systems to generate variety 
of regimes and quantum phases, including the most interesting antiferromagnetic 
(AF) regime. Garcfa-Ripoll et al. [SI] propose to use a duality between the AF and 
ferromagnetic (F) Hamiltonians, Haf = —Hp, which implies that minimal energy 
states of Haf are maximal energy states of Hf, and vice versa. Since dissipation 
and decoherence are practically negligible in such systems, and affect equally both 
ends of the spectrum, one can study AF physics with Hf, preparing adiabatically 
AF states of interest. 

• Fractional quantum Hall states. Since the famous work of Laughlin [123] , there has 
been enormous progress in our understanding of the fractional quantum Hall effect 
(FQHE) |124] . Nevertheless, many challenges remain open: direct observation of 
the anyonic character of excitations, observation of other kinds of strongly corre- 
lated states, etc. FQHE states might be studied with trapped ultracold rotating 
gases |1251ll26j . Rotation induces there effects equivalent to an "artificial" constant 
magnetic field directed along the rotation axis. There are proposals how to detect 
directly fractional excitations in such systems [127] . Optical lattices might help in 
this task in two aspects: first, FQHE states of small systems of atoms could be ob- 
served in a lattice with rotating site potentials, or an array of rotating microtraps 
(cf. [128lll29j and references therein). Second, "artificial" magnetic field might be 
directly created in an lattices via appropriate control of tunneling (hopping) matrix 
element in the corresponding Hubbard model [130] . Such systems will also allow to 
create FQHE type states [TITIIT^IT^ . 

• Lattice gauge fields. Gauge theories, and in particular lattice gauge theories (LGT) 
[134| are fundamental for both high energy physics and condensed matter physics, 
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and despite the progress of our understanding of LGT, many questions in this area 
remain open. Physics of cold atoms might help here in two aspects: "artificial" 
non-abelian magnetic fields may be created in lattice gases via appropriate control 
of hopping matrix elements |135| , or in trapped gases using effects of electromag- 
nctically induced transparency |136j . One of the most challenging tasks in this 
context concerns the possibility of realizing generalizations of Laughlin states with 
possibly non-abelian fractional excitations. Another challenge concerns the possi- 
bility of "mimicking" the dynamics of gauge fields. In fact, dynamical realizations 
of U{1) abelian gauge theory, that involves ring exchange interaction in a square 
lattice [33] , or 3 particle interactions in a triangular lattice |1371ll38j have been also 
recently proposed. 

• Superchemistry. This is a challenge of quantum chemistry, rather than condensed 
matter physics: to perform a chemical reaction in a controlled way, by using pho- 
toassociation or Fcshbach resonances from a desired initial state to a desired final 
quantum state. Ref. [139] proposed to use MI with two identical atoms, to create 
via photoassociation, first a MI of homonuclear molecules, and then a molecular 
SF via "quantum melting" . In Ref. [140] , a similar idea was applied to heteronu- 
clear molecules, in order to achieve molecular SF. Bloch's group have indeed ob- 
served photoassociation of ®^Rb molecules in MI with two atoms per site |141| . 
while Rempe's group have realized the first molecular MI using Fcshbach reso- 
nances |38j . Formation of three-body Efimov trimer states was observed in trapped 
Cs atoms by Grimm's group [142] . This process could be even more efficient in 
optical lattices |143j . 

• Ultracold dipolar gases. Some of the most facinating experimental and theoretical 
challenges of the modern atomic and molecular physics concern ultracold dipolar 
quantum gases (for a review, see |144) ). The recent experimental realisation of the 
dipolar Bose gas of Chromium |145| . and the progress in trapping and cooling of 
dipolar molecules [146j have opened the path towards ultracold quantum gases with 
dominant dipole interactions. Dipolar BECs and BCS states of trapped gases are 
expected to exhibit very interesting dependence on the trap geometry [144] . Dipolar 
ultracold gases in optical lattices, described by extended Hubbard models, should 
allow to realize various quantum insulating "solid" phases, such as checkerboard, 
and superfluid phases, such as supersohd phase [14711148] . Particularly interesting 
in this context are the rotating dipolar gases (RDG). Bose-Einstein condensates of 
RDGs exhibit novel forms of vortex lattices: square, "stripe crystal", and "bubble 
crystal" lattices [149] . We have demonstrated that pseudo-hole gap survives the 
large N limit for the Fermi RDGs [150] , making them perfect candidates to achieve 
the stongly correlated regime, and to realise Laughlin liquid at filling v = 1/3, 
and quantum Wigner crystal aX v < 1/7 |151] with mesoscopic number of atoms 
N ~ 50- 100. 



Several of the above mentioned open questions and challenges are addressed in this 
review. However, before we turn to the discussion of how ultracold atomic gases 
can mimic condensed matter systems, let us discuss shortly the properties of optical 
potentials in general, and optical lattices in particular (for a review, see [152) 1. 
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1.3 Plan of the review 

This review is addressed to two kinds of readers. First of all, it gives for condensed 
matter and, perhaps, high energy physicists an overview of what is being done in 
atomic physics and quantum optics in the area of ultracold gases in optical lattices. 
The particular emphasis is put here on the problems that are directly related to open 
problems and challenges of condensed matter, or even high energy physics. We discuss 
how to mimic condensed matter, and even go beyond toward completely new areas and 
problems. Second, the review is directed to atomic and quantum optics community. 
For these readers it should give some basic information and basic literature about 
challenging problems of condensed matter physics that can be attacked with atoms, 
ions or molecules. 

The plan of the review is as follows. In Section 2, we review the most general 
type of Hubbard-type model that can be realized with cold gases, and also review 
some spin models that can be reduced from the Hubbard model in specific limits. 
In Section 3, we present some basic theoretical methods of treatment of Hubbard 
models. Most of the material here is standard in condensed matter theory, but we 
include also a subsection about very recent developments in numerical treatments of 
many body systems based on quantum information and quantum optics ideas. In the 
following sections we address some of the challenges and open question described in 
this introduction. Each section has its own short introduction with basic condensed 
matter references to the considered problems, and then focuses on results obtained 
within the atomic physics and quantum optics context. Section 4 treats disordered 
ultracold gases, section 5 - frustrated ultracold gases, while sections 6 and 7 make 
short overviews of spinor ultracold gases and ultracold gases in "artificial" magnetic 
fields, respectively. The final section 8 discusses relations between ultracold gases and 
quantum information. 

Since the review intends to strength the analogies between condensed matter sys- 
tems and cold gases in optical lattices, it is mainly focused on the strongly interacting 
regime. Nevertheless, there is a wide range of interesting phenomena that appears in 
the weakly interacting regime that are not cover here. 

This review has been written by theorists, and as such describes experiments only 
in aspects concerning the physical results, or the experimentally accessible ranges 
of parameters. We do not discuss here experimental techniques and methods. This 
review should be considered as complementary to the excellent review by Bloch and 
Greiner [153] . 



2 The Hubbard and spin models with ultracold lattice gases 

In this section, we begin by a short discussion on optical potentials (Subsec. 12. ip . 
and then review the most general Hubbard-type model that can be realized with 
cold gases (Subsec. 12. 2p . We report then (in Subsec. [2.3^ the spin models that such 
Hubbard models reduce to, in different limits. In Subsec. [2.41 we discuss the amount of 
control that we have in the parameters involved in the Hubbard model, when realized 
with cold atoms. We then focus (in Subsec. [^75|l on the the paradigmatic model of a 
system that exhibits a quantum phase transition, namely on the homogeneous Bose- 
Hubbard (BH) model [5011154] . The model undergoes the superfluid-Mott insulator 
(SF-MI) quantum phase transition. 
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Figure 2. Schematic pictures of optical lattice potentials; a) 2D square lattice of quasi ID traps; b) 

3D simple cubic lattice (from [33]). 

2.1 Optical potentials 

The basic tool to create ultracold lattice gases are optical potentials. An electron in 
an atom in the presence of oscillating electric field E(r, t) of a laser attains a time 
dependent dipole moment d. When the field oscillations are far off resonance (i.e. they 
do not cause any real transition in the atom) . the induced dipole moment follows the 
laser field oscillations, 

= a^j{ujL)Ej{r,t), (1) 

where di is the corresponding component of d(i — x,y,z), ujl is the laser frequency, 
and Qfy (wl) denotes the matrix elements of the polarizability tensor. The polarizabil- 
ity depends in general on the laser frequency, and on the energies of the non-resonant 
excited states of the atom. One of these states (with excitation energy, say, Ei = huji) 
is usually much closer to the resonance than the others; in such case the polarizability 
becomes inversely proportional to the laser detuning from the resonance, A = lol — ijJi- 
Electronic energy undergoes in this situation a shift, AE, which is nothing else than 
a AC-version of the standard quadratic Stark effect. The energy shift is proportional 
to 

AE{r)= a.A^Lm{r,t)E,{r,t))^I{r)/A, (2) 

where the bra-kct denotes the averaging of the product of electric fields over the fast 
optical oscillations, and /(r) is the laser beam intensity. 

The consequences of the above simple formula are enormous. The atom feels an 
optical potential Vopt{^) = AE{r), that follows the spatial pattern of the laser field 
intensity. This is the basis for optical manipulations and trapping of atoms! If the 
laser is red-detuned (i.e., laser detuning A < 0), the atoms are attracted toward the 
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regions of high intensity; conversely, blue detuned laser pushed the atoms out of the 
regions of high intensity. 

Adding two, or more laser fields of the same frequency leads in general to interfer- 
ences and the corresponding interference pattern of the intensity. In particular, two 
counter propagating laser waves of the same polarisation will create a standing wave, 
and thus a spatially oscillating potential for atoms. One can easily avoid interferences 
when adding more and more laser fields, and add corresponding intensities, i.e. optical 
potentials. To this aim one can use laser fields with orthogonal polarisations. Three 
pairs of counter-propagating laser beams with orthogonal polarizations will form then 
the 3D optical lattice represented schematically in Fig.l^b. An alternative way to avoid 
interferences on demand is to use slightly, but sufficiently different frequencies. In this 
case, time averaging over the "sufficient" different frequencies washes out the inter- 
ference effects. Similarly, one can use laser beams with polarisations oscillating at 
different frequencies to avoid interferences. 



2.2 Hubbard models 

We are interested in the Hubbard-type models that are realizable with cold atoms in an 
optical lattice. The simplest optical lattice is a 3D simple cubic lattice, as presented in 
Fig. [^h- It is formed by three pairs of laser beams creating three orthogonal standing 
waves with orthogonal polarisations. As we will discuss below, one can, however, 
create practically arbitrary lattices on demand using optical potentials. Also, as the 
intensity of one of the standing waves increases, the probability of hopping along this 
direction decreases rapidly to zero [29]. In effect we obtain an ID array of 2D square 
lattices. Consequently, an increase of the laser intensity of another of the standing 
waves, creates effectively a 2D array of ID lattices (Fig. [2^). 

Optical lattices provide an ideal (contain no defects) and rigid (do not support 
phonon excitations Q) periodic potential in which the atoms move. As it is well known 
from solid state theory [155j , single particle energy spectrum (in the absence of inter- 
actions) consists of bands, and the energy eigenstates of the Hamiltonian are Bloch 
functions. If the lattice potential is strong, the band-gaps are large, and the bands are 
very well separated energetically. For low temperatures regime, it is easy to achieve 
a situation in which only the lowest band is occupied (tight binding approximation). 
The Bloch functions of the lowest band can be expanded in Wannier functions, which 
are not the eigenstates of the single particle Hamiltonian, but are localised at each 
site. In the tight binding approximation, we project all of the atomic quantum field 
operators that describe the systems in question onto the lowest band, and then expand 
into the Wannier basis. Leaving then just the most relevant terms in the Hamiltonian 
(such as hopping between the nearest neighboring sites) leads directly to Hubbard- 
type Hamiltonian. 

Let us write the most general Hubbard model that may be created in this way. Let 
us assume that we have several bosonic and fermionic species (or bosons/fermions 
with several internal states), enumerated by a. The basic objects of the theory will 
be thus annihilation and creation operators of a-bosons, bia,b\^ at the site z, and 
analogously annihilation and creation operators of /3-fermions, Jip^jla at the site i. 



^ This statement has to be revised when the lattice is created inside of an optical cavity. As we discuss 
later, the presence of atoms may affect the cavity field. 
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Bosonic (fcrmionic) operators fulfill, of course, the standard canonical commutation 
(anticommutation) relations: 



{fia I fja' } 



^ fia 1 fjai ^ — ^aa' &ij j 



(3) 
(4) 



where 5ab denotes the Kronecker delta. 

The most general Hubbard type of hamiltonian that one can realize with cold atoms, 
assuming lowest band occupation onlyQ, consists of four parts; 



H 



Hubbard 



H. 



hop 



Hi, 



H 



pot 



H 



Rabi ■ 



(5) 



The hopping part describes hopping (tunneling) of atoms from one site to another. 
Since the hopping probability amplitude decreases exponentially with the distance, 
hopping is typically assumed to occur between the nearest neighboring sites, denoted 
{ij). Hopping, on the other hand, might be laser assisted, and thus might allow for 
transitions from one internal state to another (of course, or perhaps unfortunately, it 
cannot lead to a change of element, or even isotope!). 



Hhop 



h.c 



(6) 



Atoms interact in the first place via short range Van der Waals forces, which in 
the low energy limit are very well described by various kinds of the zero range pseu- 
dopotentials [T7lll57[fl58| . That means that the dominant part of the interactions 
is of contact type, i.e. occurs on-site. There are, however, situations in which inter- 
actions do affect neighbouring sites, or even have long range (such as dipole-dipole 
interactions). We thus write the interaction part as 



Hint ~ H„ 



He. 



(7) 



where 



H„ 



(8) 



In the simplest Hubbard models interactions depend only on the on-site atom num- 
bers, nf^ = bl^bia, and nf^ = fj^fia- In general, however, they may depend in an 
non-trivial way on internal states, or atomic species; this is for instance the stan- 
dard case for spinor gases |117| . The models with non-contact interactions are usually 
termed as extended Hubbard models, and include 



H, 



E 



Val3a'l3'ihj)bLbj^bj[3>b,c,' + V^f[^,p,{i,j)flj]pfj[3'fi, 



^Some authors go beyond this assumption. See for instance Ref. 11561 . 



15 



+ 2VX'P'ih3)blfy,0'hc.'\ . (9) 

Most typically, the non-contact interactions will depend on the distance between the 
sites, |ri — rj|, and will be of the density-density form (i.e. they depend only on n^, 
and nfp), but in general, again this does not have to be the case. Dipolar interactions, 
for instance, depend on the angles between the dipole moments and on the vector 

The last two parts of the Hamiltonian (Eq. ([5])) describe on-site single atom pro- 
cesses, and essentially have the same form as the tunneling part, except that they 
occur on-site. We do distinguish them since they arc well defined and controllable 
in experiments. Hpot combines effects of all potentials felt by atoms such as external 
trapping potential, possible additional superlatticc (i.e. additional lattice) potentials, 
possible disorder potentials, and last, but not least, chemical potential which is nec- 
essary if one uses the statistical description based on the grand canonical ensemble: 



Hpot 



EK"4&«+/^fa/l/»] . (10) 



The last part of the Hamiltonian, HuaM, describes possible coherent on-site transi- 
tions between the internals states of atoms; such transitions may be achieved using 
laser induced resonant Raman transitions, or microwave Rabi type transitions (we 
can write this part of the Hamiltonian as time-independent in the interaction picture 
with respect to the on site internal states Hamiltonian) 

HRaM = IY. [^^cpblhp + nf^pflu] . (11) 

Q,/3,i 



2.3 Spin models 

As it is well known (see for instance |60|). Hubbard models reduce to spin models in 
certain limits. Most of these limits, and even more, are accessible with cold atoms. 
Generally speaking, if bosonic atoms can occupy only 25* -I- 1 different states in a 
lattice site, then one can always map these states onto the states of pseudo-spin S. 
These states may even correspond to different number of bosons, and the dimension 
of the local, on-site Hilbert space might vary in space; in such case we will deal with 
inhomogencous models, where at each site there is, in general, a different spin. Similar 
construction might be done for fermionic atoms, with the remark that at a given site, 
the fermion number differences might attain only even numbers, since otherwise the 
fermionic character of particles cannot be eliminated. 

When constructing specific spin models two aspects play a role: lattice geometry 
(which we discuss in the next subsection), and the form of interactions, which includes 
Ising, XY, Heisenberg, XXZ, and anisotropic XYZ types, as well as ring exchange 
types. Below we list the most obvious constructions of spin models, that has been 
discussed in the literature on cold atoms. We restrict ourselves here to translationally 
invariant models. 

• Hard core bosons and XY models. Perhaps the simplest way to obtain a non-trivial 
spin model is to use the simplest Bose-Hubbard Hamiltonian for one component 
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"spinlcss") bosons 

, [6!6j+ /J.C.J + ^ 

iv) 



H = -tJ2[blbj + h.c^ +^^7i,(n, (12) 



where (ij) denotes sum over nearest neighbors. In the hard boson Hmit (i.e. when 
U ^ t,fj,) we may have at most 1 boson per site. We may encode the spin 1/2 
states as presence (t), or absence (|) of the boson at the site. The Hamiltonian 
reduces then to that of XY model in a transverse field, 

where ct,; = {ax,i + iay^i)/2 and <Jx,y,z,i denote the standard Pauli matrices at site i. 
This model has the advantage that in ID it is exactly solvable via Jordan- Wigner 
transformation |154| . One interesting application of this approach concerns the ID 
disordered chain studied in Rcf . |159| . The same approach was used recently in [55] 
to realize XY model in random parallel field. 

Spatially delocalised qubits. Somewhat similar idea considers two neighbouring 
traps, or potential wells ("left" and "right"), assumes 1 atom per double well, and 
encodes the spin 1/2 (qubit) as the presence of the atom on the "left", or "right", 
respectively |160| . In Ref. [SS] was proposed to encode one qubit by the presence 
or the absence of a whole string of neutral atoms yielding improved robustness. 
This system may be used for generation of maximally entangled many atom states 
(Schrodinger cat states) by crossing a quantum phase transition. 

Multi- component atoms in Mott states. Whenever we deal with a system of multi- 
component atoms, i.e. atoms with say 25* + 1 internal states, in the Mott insulator 
limit, the system will be well described by the appropriate spin model. The most 
prominent example is a two-component (or spin 1/2) Fermi gas |60] . which in the 
Mott state with one atom per site forms a perfect Heisenberg model. Several groups 
are planning experiments with ultracold spin 1/2 Fermi atoms heading toward an- 
tiferromagnetism in various kinds of lattices. Prospects for observing antiferromag- 
netic transition in such systems are quite good, especially since one expect to be 
able to employ interaction induced cooling (an analogue of the Pomeranchuk ef- 
fect, known in liquid Helium physics), |161j . or disorder induced increase of Tc [55] . 
One should stress, however, that although the Mott transition takes typically place 
when t <U , the low temperature physics of the resulting Heisenberg model requires 
temperatures of order /U ~ fc^T. Such temperatures arc often in the nanoKelvin 
range, i.e. still hard to achieve experimentally. There is, however, a lot of new pro- 
posals of cooling atoms in the Mott states |1621I163[|1641I165] . and hopefully the 
temperatures will not set any limitations on experiments with these kinds of spin 
models in the next future. 

The calculation of the effective Hamiltonians in pseudo-spin 1/2 Bose-Bosc or Bosc- 
Fermi mixtures is quite complicated, and was accomplished for the first time re- 
cently |166lll67j . The Bose-Bose case can be reduced to an XXZ spin model (for the 
case of XXZ model in random fields see for instance [55])- The effective Hamiltonian 
for the Fermi-Bose mixture in general cannot be reduced to a spin model, since it 
involves Fermi operators describing composite fermions, consisting of one fermion 
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paired with some number of bosons, or bosonic holes |167j . The Hamiltonian de- 
scribes a "spinless" interacting Fermi gas of such composites. It can, however, be 
transformed to an XXZ model in external fields in ID via Jordan- Wigner transfor- 
mation. 

• Spinor gases in Mott states. Of course, the above statements are particularly valid 
for spinor gases, which for atoms with spin F > 1/2 have effective Hamiltonians 
containing generalizations (a power series) of Heisenberg interactions. For F = 1 
and the Mott state with 1 atom per site, we deal with the so-called Quadratic- 
Biquadratic Hamiltonian [TT51[ngi[T^01ll21l[S^ . In the Mott state with two atoms 
per site, the pair can either compose a singlet state or a state with on-site spin 
S ~ 2. The resulting Hamiltonian contains then higher powers of the Heisenberg 
term. The situation is obviously more complicated for higher Mott states, and atoms 
with higher individual spin F. In Refs. [122lll68lll69j . the case = 2 is studied; 
here, already with one atom per site, the effective Hamiltonian is a polynomial of 
the fourth order in Heisenberg term. 

• Spin models in polymerized lattices. Yet another interesting way to obtain spin 
models with cold atoms, is based on the use of polymerized (dimerized, trimerized, 
quadrumerized etc.) lattices. These are lattices that can be easily realised with 
optical potentials, and have no analogue in condensed matter physics. A simple 
example of dimerized lattice in 2D is a square lattice of pairs of close sites; trimer- 
ized kagomc lattice, discussed in section 4, is a triangular lattices of trimes of close 
sites located on a small unilateral triangle; 2D quadrumerized square lattice is a 
square lattice of small squares of close sites, etc. When one considers ultracold 
gases in such lattices, one has to take into account first the lowest energy state 
in a dimer, trimer etc. If we deal with polarized "spinless" fermions in trimerized 
lattice, and we consider two fermions per trimer, fermions have to their disposal 
zero momentum state (which will be necessarily filled at low temperatures), and the 
two states with left and right chirality. The latter two are obviously degenerated, 
and can thus be used to encode the effective spin 1/2. This is the model discussed 
in Refs. [551lll0irillj . Note, that similar encoding is possible in a quatromerized 
lattice with 2 atoms per quadrumer. 



2.4 Control of parameters in cold atom systems 

Atomic physics and quantum optics offer many new types of methods to quantum 
engineer systems in question. Toutes proportions gardees, there are instances in which 
atomic physics and quantum optics not only meets, but rather "beats" condensed 
matter physics. That is one of the reasons why the physics of ultracold atoms attracts 
so many theorists from other disciplines. 

Let us list shortly what can be controlled in the experiments with cold atoms in 
optical lattices: 

• Lattice geometry and dimensionality. As mentioned above, practically any lattice 
geometry may be achieved with optical potentials. The method of superlattices (i.e. 
adding a new lattice on top of the exiting one) is very well developed. Changing 
of the lattice dimentionality does not pose any problem (compare Fig. [5]). Also 
periodic boundary conditions can be realized in ring shapeed optical lattices [170j . 

• Phonons. Optical lattices are rigid and robust: they do not have any phonons. An 
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interesting situation arises when the lattice is formed in an optical cavity: atom- 
light coupling might suffice then to shift the cavity resonance. Cavity will affect the 
atoms, but atoms will perform back action, and create "phonon" like excitations. 
For early works, see [17111172] . For more recent studies of superfluid-Mott insulator 
crossover, see |173| . and for the first attempts towards "refracton" physics (analogs 
of "phonons"), see |174| . Phonons, obviously, play a role in ion traps, where they 
provide the major mechanism for ion-ion interactions. 

• Tunneling. Tunneling can be controlled to a great extend using combination of pure 
tunneling, laser assisted coherent transitions, and lattice tilting (acceleration) tech- 
niques. The prominent example of such control describe the proposals for creating 
artificial magnetic fields [MllfmifmfminSj . 

• On-site interactions. These interactions are controlled by scattering lengths, which 
can be modified using Feshbach resonances in magnetic fields, [551IM1I55] . or optical 
Feshbach resonances (for theory see |175j . for experiments [1771ll76j l. On-site in- 
teractions can be set to zero in dipolar gases, by changing the shape of the on-site 
potential [147j . 

• Next neighbour and long range interactions. Effective models obtained by calculat- 
ing effect of tunneling in the Mott insulator phases, contain typically short range 
interactions of energies oc t^ /U . Stronger interactions can be achieved using dipolar 
interactions, such as those proposed in Refs. [147llll5lll78j . Dipolar interactions are 
of long range type, are anisotropic, and exhibit a very rich variety of phenomena 
(for a review, see 11441 ). They can also be achieved in trapped ion systems, where 
they are mediated via phonon vibrations of the equilibrium ionic configuration. 
This case is discussed in detail in the section [7] of this review. The group of T. 
Pfau has recently realized the first experimental observation of ultracold dipolar 
gas jl45| , by condensing bosonic Chromium. Dipolar interaction are mediated here 
by magnetic dipoles of the Chromium atoms; they are weak, but nevertheless lead 
to observable effects |179| . 

• Multiparticle (plaquette) interactions, ft has been also demonstrated how to gener- 
ate effective three-body interactions in triangular optical lattices |137| . These in- 
teractions result from the possibility of atoms tunneling along two different paths. 
Similarly, ring exchange interactions in square optical lattice can be generated em- 
ploying the correlated hopping of two bosons [15] . 

• Potentials. Various types of external potentials can be applied to the atoms, de- 
pending on the situations. One can use magnetic potentials whose shape can be 
at least cotroUed on the scale of few microns. Magnetic potentials with larger gra- 
dients can be created on atom chips (cf. |180] V The most flexible are, however, 
optical potentials. Apart from limitations set by the diffraction limit, they can 
have practically any desired shape and can form any kind of optical lattice: regu- 
lar, disordered, modulated, etc. Recently, J. Schmiedmayer demonstrated also great 
possibilities offered by the so called radio frequency potentials [181] . 

• Rahi transitions. Similarly, apart from limitations set by diffraction, they are in 
practice highly controlable. 

• Temperatures. The typical critical temperatures T of trapped ultracold condensed 
Bose gases are of order of nano-Kelvins. Using evaporative cooling one can reach, 
however, lower T (which in fact are not very well known, because of the lack of 
reliable temperature measurement methods. For recent advances see |182] ). Sim- 
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ilarly, the temperature of superfluid Fermi gases are in the range of tens of nK. 
One can thus say that temperatures in the range of tens of nK are becoming nowa- 
days a standard. There are many proposals for reaching even lower Ts employing 
additional cooling and filtering procedures |165| . SF-MI transition occurs in the 
regime of Ts accessible nowadays. Many of the strongly correlated phases occur in 
the regime when the tunneling t is much smaller than U and require temperatures 
of order ksT ~ t'^ /U, i.e. 10-20 nK, or even less. This is at the border of the cur- 
rent possibilities, but the progress in cooling and quantum engineering techniques 
allow us to believe that these limitations will be overcome very soon (for a detailed 
discussion see |183j ) 

• Time dependences. The time scales of coherent unitary dynamics of these systems 
are typically in the millisecond range. It implies that, in contrast to condensed 
matter systems, all of the controls discussed above can be made time dependent, 
adiabatic, or diabatic, on demand. Some of the fascinating possibilities include 
change of lattice geometry, or turn-on of the disorder in real time. 

The huge range of parameters which are experimentally controllable indicates the 
rich possibilities offered by cold gases in optical lattices to implement condensed 
matter models and beyond. The first proof came with the seminal paper of Greiner 
et al. j32j reporting the superfluid-Mott insulator transition with cold bosons in an 
optical lattice. 



2.5 Superfluid - Mott insulator quantum phase transition in the Bose 
Hubbard model 

Let us now consider the ideal homogeneous Bose-Hubbard (BH) model of the form 
H = -tY. (bib, + h.c.) + ^ ^ n,(n, - 1), (14) 

where («, j) indicates sum over nearest neighbors. We denote here the tunneling energy 
by t (both t and J are used in the literature). Below we discuss the superfluid (SF) - 
Mott insulator (MI) quantum phase transition separately for the case of transitions 
at fixed density, and for the, so-called, generic (density-driven) transitions, where the 
number of atoms changes: 

• Transitions at fixed integer density. We consider here transitions driven by a change 
of t/U ratio in a system with a fixed number of bosons. The phase transition occurs 
when the lattice filling factor n (the number of atoms per site) is exactly integer. 
For t/U < {t/U)c < 1 there is a Mott insulator (MI) phase, while for t/U > 
{t/U)c there is a superfluid (SF) phase. As discussed by Fisher et al. [3D], such a 
transition in a d-dimensional BH model lies in the universality class of the {d+ 1)- 
dimensional XY spin model. This result implies that the one dimensional Bose- 
Hubbard model undergoes a Kostcrlitz-Thoulcss phase transition. Also it permits to 
determine the critical exponents. The quantum phase transition happens, ideally, at 
absolute temperature equal to zero |154| , and its signatures are reflected in different 
quantities as discussed below. 

The superfluid fraction |184[|185] . which is defined through the response of the 
system to an externally imposed velocity field (equivalent to a twist in boundary 
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conditions) vanishes in the Mott phase. This was verified in a ID system numerically 
by a DMRG calculation [186j . and analytically in a simplified BH model subjected 
to the restriction of a maximal site occupation of two particles |187j . Additionally, 
a jump in the superfluid fraction was observed at the transition point |186| of a one 
dimensional system, a result expected at the Kosterlitz-Thouless critical point |188j . 
It is worth to stress here, that the superfluid fraction should not be confused with 
the condensate fraction [iD]. The latter one is equal to the highest eigenvalue of 

the single particle density matrix {b\bj)^ divided by the number of particles. These 
quantities, both equal 100% for an untrapped 3D Bose-Einstein condensate in the 
dilute limit jl85| , can be very different in the Bose-Hubbard model, see e.g. [55] 
for a simple example where the condensate fraction is 100%, while the superfluid 
fraction hits zero once a sufficiently strong on-site disorder is present. 

The excitation spectrum is gapless for the SF phase while gapped for the MI. In 
the MI neighborhood of the transition point, the gap scales exponentially in a one 
dimensional system as cx exp(~ [[t/U)c — t/U]~^^'^), where the proportionahty 
factor in the exponent is smaller than zero. In two and three dimensional models, 
it exhibits a power law behaviour oc [{t/U)c — t/Uy^ ^ where z and v are the 
critical exponents [1541 [5U] . These exponents for a two dimensional system are 
z ~ 1 and V k, 0.67 [189| . and for a three dimensional model they read as z = 1 
and V = 1/2 [3D]. 

Another quantity that shows a critical behaviour is the correlation length ^: 

e = l/2Y,r^b]bj+r + h.c.)/Y,{b]b,+r + h.c). (15) 

r r 

It is finite in the Mott phase, diverges at the critical point and stays divergent 
in the superfluid phase. In the neighborhood of the critical point it behaves as: 
^ ~ gap~^/^, where the critical exponent z = 1 in one, two, and three dimensional 
systems [30] . 

Additional insight into the superfluid-Mott insulator phase transition can be ob- 
tained investigating the relation between /i/J7, t/U , and the number of particles 
in the system. The chemical potential ^ is conveniently introduced during min- 
imization of {H — fi rii) leading to the determination of ground state with 
a /x-dependent total number of atoms. When t/U = 0, one obtains that for 
h — 1 < n/U < n {n is an integer number) the ground state is a single Fock 
state 

...), (16) 

When tunneling is nonzero, the range of ^/U describing the system with integer 
n gradually shrinks and finally disappears at some {{^/U)c, {t/U)c\- In this way, 
the famous lobes are formed: see Fig.[3| The fixed density transition happens when 
the system moves along the thin solid lines depicted in Fig. [3] This schematic 
plot illustrates also that a system with a non-integer filling factor never enters the 
MI lobes, i.e., stays always superfluid as depicted by the dashed lines. It is easily 
understood by considering a state with one particle added (subtracted) to (from) a 
system having integer filling factor. Such a particle (hole) can freely flow through 
the lattice, so that the system becomes superfluid for all values of t/U ratio. 

It is important to stress that though inside the lobes the flUing factor n is integer, 
the ground state is not a single Fock state ([T5|) for t/U ^ 0. To illustrate this fact 
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Figure 3. A schematic plot of the phase diagram of the Bose-Hubbard model. The lobes, 
surrounded by the superfluid sea, correspond to the Mott insulator islands with integer filling 
factor n. The thin solid lines represent the lines of constant, integer density. The dashed lines show 
trajectories of a system with fixed, non-integer, filling factor n (0 < e <IC 1). The dotted line 
presents an example of a trajectory leading to a generic phase transition when the system enters 
either Mott, or superfiuid phase by changing a total number of atoms. 

one can look at the expectation values of some operators at the transition point 
of a ID BH model at fi = 1. The ground state wave-function of this system has 
presumably pronounced deviations from (|16|) . For instance, one obtaines there that 
(i) the nearest neighbor correlation length function (bjbj+i) is approximately 0.8 
and (it) the variance of on-site number operator, [{nj) — (rij)^]^/^, equals about 
0.6 |190| . These quantities significantly differ from the predictions obtained from 
the ground state for t/U = ([TB)) . Notice also that they attain substantial values 
since they arc bounded by unity for any t/U. 

It is interesting to ask what are the critical points {t/U)c in different dimensions, 
and at different filling factors. Here we list the most accurate estimations up to 
date for the case fi = 1 which has been systematically studied in the past. In 
one dimension, the critical point was precisely determined by DMRG analysis: 
(t/U)c ~ 0.29 |191| . In two dimensional system, a recent quantum Monte Carlo 
studies estimate the position of the critical point, {t/U)c, to be around 0.061 [192j . 
In the three dimensional model, the perturbative expansion gives {t/U)c ~ 0.034 
|193| . The locations of the critical points for different small filling factors in different 
dimensional lattices were quite accurately calculated in Ref. [193] . 

• Generic transitions. A quick look at the dotted line in Fig. [3] reveals that the system 
can cross the superfluid - Mott insulator phase boundary through trajectories that 
do not correspond to a fixed filling factor. In such a case we say that the system 
undergoes the generic (density-driven) SF - MI phase transition. This transition 
does not belong to the universality class of the XY spin model, thus it is charac- 
terized by different critical exponents. In particular, one finds that zv = 1 during 
generic transitions in ID, 2D and 3D systems [30] . 

Since for this case the number of particles is not conserved, one should add the 
term —fJ-'^^ni to the Hamiltonian (|14p . and find a ground state with a number 
of particles depending on the chemical potential /x. The critical behaviour can be 
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observed in at least two quantities: the compressibility k = dp/d^i, (where p is atom 
density), and the superfluid fraction pg. The first one, diverges as one approaches 
the Mott lobes from a superfluid side (Fig. [3]), while the second goes to zero in 
this limit. Both the compressibility and the superfluid fraction stay zero inside the 
lobes. 

To illustrate these statements let us consider a one dimensional system, at fixed 
t/U ^ undergoing a phase transition induced by a change in the number of atoms. 
In this case, the theory of Fisher et al. [3D] predicts that (i) k ~ |/i — /ij"^'^"^) 
(/X is a chemical potential at phase boundary); (ii) ps ^ \p — pY~^ {p is a filling 
factor of a Mott lobe approached during transition). The early Quantum Monte 
Carlo simulations |1941I195] have verified these predictions giving the following 
estimations of the critical exponents: z ~ 2.04 and v « 0.48. 

Recently, there has been a lot of interest in generalizing MI-SF transition and 
Bose-Hubbard model to more "exotic" situations, such as atoms in optical cavities 
|174lll73j . or p band Hubbard model, where transverse staggered order occurs |196j . 
In triangular lattices, stripe order was predicted both in SF and MI phases [197] . 



3 The Hubbard model: Methods of treatment 
3.1 Introduction 

In this section we review at a rather basic level some of the standard theoretical tools 
used in condensed matter theory to treat many body systems of interest. We discuss 
also novel developments which - taking advantage of quantum information methods - 
provide an efficient way to calculate ground state properties and dynamical evolution 
of many condensed matter systems. The underlying philosophy of these new meth- 
ods is closely related to the well-established Density Matrix Rcnormalization Group 
(DMRG) method |231[|47j . and consist in truncating the dimension of the Hilbert 
space, which diverges exponentially with the size of the system, to a manageable size, 
considering entanglement properties of different bipartite partitions. 

Analytical and numerical methods often rely on the size and dimensionality of the 
system. Powerful techniques like bosonization pSll48l[^5^ , Bethc ansatz |233[|234lH5] , 
Jordan- Wigner transformation [154,235 , or the mentioned DMRG exist and allow to 
solve some paradigmatic one dimensional systems, such as for instance Hcisenbcrg spin 
1/2, XXZ spin chains, or ID Hubbard model of strongly interacting electrons. These 
methods are very well established to treat many body problems in one dimension, 
but often fail in higher dimensions. Finite size effects are also crucial in the study 
of strongly correlated systems, because quantum phase transitions occur only in the 
thermodynamic limit at zero temperature. It is thus important to know how finite 
size effects affect the statics and dynamics of strongly correlated systems, and the 
signatures of quantum phase transitions. 

We begin this section by focussing on methods of treatment of the ideal homoge- 
neous Bose Hubbard model. As we have discussed in the previous section, this system 
exhibits a phase transition: the superfluid-Mott insulator (SF-MI) quantum phase 
transition. Despite its simplicity, the interplay between tunneling and on-site inter- 
actions in the BH model is by no means trivial. Excitations in the limit of small 
interactions can be described with the help of Bogoliubov transformations. In the 
strong interaction limit, when the ratio between tunneling and on-site interactions is 
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much smaller than one, tunneling can be treated as a perturbation, and the original 
Hamiltonian can be replaced by an effective Hamiltonian. In some cases, it is also 
possible to use a mean field approach, like for instance the Gutzwiller ansatz, which 
assumes that the many body wave functions have a product-over-sites form, and is 
conceptually and numerically relatively easy to implement. As we shall see, it predicts 
quite correctly the critical points separating the Mott phase from the SF phase for 
3D, or even 2D lattices, but its accuracy decreases dramatically for ID systems. We 
discuss then briefly exact diagonalisation. Quantum Monte Carlo, and phase space 
methods. Subsequently, we discuss very shortly ID methods: Bethe ansatz, Jordan- 
Wigner and bosonisation. We analyze in more detail the novel approach to DMRG 
provided by Quantum Information. Finaly, we discuss some methods to treat Fermi 
and Fermi-bose Hubbard models. 



3.2 Weak interactions limit 

Here we want to illustrate how the Bose-Hubbard model can be solved in the limit of 
small interactions. Our discussion follows Ref. [236], where the Bogoliubov approach 
was developed. 

After adding the — A^Si'^i term to Eq. (jl4[) . the Hamiltonian can be conveniently 
written as: 

H = -tJ2 (^Ih + h-^-) + ^J2 ^\^\^^^^ -f^T^^lK (17) 

where /i denotes the chemical potential. We assume a regular d-dimensional lattice, 
consisting on M sites, and the distance between neighboring sites is "a" . In the limit 
of t/U oo, interactions between atoms are negligible. The system is completely 
condensed in the ground state, and A'^o (the number of condensed atoms) equals A^ (the 
total number of atoms). When interactions become non-negligible, atoms gradually 
leave the condensate. To describe this process, it is convenient to work in momentum 
space: bj ~ -^^'^j:exp{—ikxj)aj:, where Xj points into j-th lattice site, and k is 

discretised over the first Brillouin zone. Using the identity e'^^'^^^'^^-'^* ~ ^^^ki k^ 
one obtains: 



H = [—£(&) — u ata,- H (5,- , r r ( a- ar ar , (18) 

Z-^ [ ^ ' ' J fc 2M ^ ki+k2M+k4 ki k2 '=3 fc4 ' ^ 
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where e{k) = 2tJ2i=i'^os{kia). As t/U oo, the ground state converges towards 

~ Oq^IO). The Bogoliubov approach relies on the transformation qq — s- \/No + oq, 
where the new operator oq is responsible for fluctuations of the number of condensed 
atoms. Substituting the above expression in Ea. p8|) . one finds, up to the quadratic 
terms: 



H ^ N< 



i^zt - f-i+ -^no^ + v/Ao(C/71o -tz- ^)(ao + Oq) 



^(-e(fc) - M)ata£ + ^ E^H^fe + + ""h-k^' ^^^^ 
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where no = Nq /M is the condensate density and z = 2d. Setting the chemical poten- 
tial to = Un^ — zt removes the linear part while the quadratic one is diagonalized 
by the Bogoliubov transformation: = Uj^aj: + vj:a} Notice that 



\urA - \vt 



= 1 



from the requirement that [cg, ct] 1. After a simple algebra, one obtains that within 
the quadratic approximation the Hamiltonian reduces to: 



H ^ - 



UnoNo 



k 



2^' 

k 



(20) 



if 



(|Ufe|^ + |wj:|^) UnQ-e{k)+zt - Uno{vj:u*j^ + v'j:Uj:) = hujj: 



^k^k = 0- 



{ui + v^)Uno-2 Uno-e{k) + zt 
Assuming that and vj: are real, one easily obtains from these equations: 



\J [zt - e(fc)]2 + 2Uno[zt - e(fc)], 
zt — e{k) + Uno 



k k 2 



^k 



(21) 
(22) 



These results reveal that the excitation spectrum is gapless in the thermodynamic 
limit M, TV ^ oo at M/N being fixed. Indeed, for the long wavelength (phonon) 

modes (|fc|a ^ 1) we find that: 



^k 



\k\aVi^Jt\k\^a^ + 2Uno 



i.e., the energy of a single excitation can be arbitrarily small: an expected result 
showing that the Bogoliubov approach does not work in the Mott phase. 

At zero temperature, there are no excitations in the system, so that the ground 
state is a Bogoliubov vacuum |vac), such that cg|vac) = 0. At finite tempera- 
ture, say T, excitations are present, and occupations of different modes satisfy 

(ctcjr) = exp (J^^ — 1 , in accordance with the Bose-Einstein statistics. Us- 
ing these properties and the solutions ((2T|) and ([22]) one can easily calculate different 
quantities (e.g., correlation functions, number of condensed atoms, etc.) both at zero 
and finite temperatures. It has to be remembered, however, that reliable predictions 
can be obtained as long as No{T) ^ N, which can be self-consistently verified within 
this approach. It is also worth to stress that the Bogoliubov approach can be applied to 
time- dependent problems without further complications. Time dependent Bogoliubov- 
de Gennes method, together with variational approach and the Kibble-Zurek mech- 
anism has been recently used to show the scaling behaviour of the time-dependent 
correlations [237j . 
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3.3 Strong interactions limit 

Let us now study the Bose- Hubbard model in the hmit of stong interactions, i.e., 
when the system is in the Mott phase. A systematic approach for studies of the Mott 
insulator phase is provided by the strong coupling expansion, i.e., a perturbative 
expansion in t/U of Eq. p^ . 

To perform strong coupling expansion one splits the Hamiltonian into two parts: 
Hq, whose eigenstates are exactly known at t/U = 0, and treat —t/UJ2{i j) (^I^i+'^-C-) 
the tunneling part, as a perturbation. Within this approach, the expectation values 
of some operators are expressed as a series of the form J^i '^i {t/Uy . Expansions up 
to 14th-ordcr have been calculated, which guarantees in some cases, a high accuracy 
of perturbative predictions. High order calculations can be performed symbolically on 
a computer, so that the expansion coefficients (ai's), can be obtained exactly (not as 
the double precision numbers). The theoretical background for these calculations was 
set up in [23811239] . where this method was applied to spin systems. Below we review 
the relevant results in the context of the Bose-Hubbard model. 

To start with, one can use the strong coupling expansion to determine the phase 
boundaries on the (/i/t/,t/C/) plane (Fig. O [^[TMIM^IMlfTM] . In this case, the 
unperturbed bare hamiltonian is = i ni{ni — l) — ii/lf^^ n,;, and one calculates 
perturbatively (i) the energy of the ground state with exactly n atoms per site; (ii) the 
ground state energy of the system with one particle added (subtracted) to (from) the 
system with filling factor h. Setting the energy difference between (i) and (ii) cases to 
zero, one obtains the value of the chemical potential at the upper (lower) boundary 
between insulator and superfluid phases. These calculations can be performed for any 
dimensional lattice, and the order of expansion can be as large as 13th j240] . For one 
dimesional systems |242ill91] . the predicted structure of Mott insulator lobes is in 
a very good agreement with the numerical results obtained via DMRG calculations. 
The perturbative expansions can also be used together with different extrapolation 
methods leading to the determination of the critical exponents z and ly |193[|240j . 

Very recently this method has been applied to calculate the lobes for a modified 
ID Bose-Hubbard model describing atoms in an optical lattice created by pumping 
a laser beam into the cavity [173] . The major difference to the standard case is that 
the intensity of the cavity field depends on the number of atoms present, since the 
atoms shift collectively the cavity resonance. In effect the coefficients t and U become 
very complex functions of all of the relevant parameters; cavity detuning, intensity 
of the pumping laser, N, etc. Moreover, quantum fluctuations of the resonance shift 
induce long range interactions between the atoms. The phase diagram, as a function 
of the dimensionless parameters fi/U and k/t], where k is the cavity width and 77 is 
the pumping laser strength, is shown in Fig. 01 ( 1 / ^nph = k/t], where Uph is the 
number of photons in the cavity.) The striking effect is the overlap of different Mott 
phases, which is the consequence of the fact that the expressions for t and U for 
no = 1, 2, 3, . . . Mott phases are different. 

The strong coupling expansion has also been used to calculate the correlation func- 
tions {h\bi+j) and the structure factor [240[|190| . A typical prediction of the strong 
coupling expansion reads as 
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Figure 4. Overlapping MI phases for the Bose-Hubbard model in an optical cavity (from 11731 ). 
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which was obtained in a one dimensional system at unit fiUing factor in |190j . The 
differences between such analytical result and the numerical calculation turn out to be 
hardly visible for t/U < 0.3, i.e., in entire Mott phase. The strong coupling expansion 
has been also employed to determine density-density correlations, {mrii+j), and the 



variance of on-site atom occupation [{r 



[Ho]. All these quantities should 



be directly measurable in an ongoing experiment in a homogeneous, one dimensional 
lattice pjg] . 



3.4 The Gutzwiller mean-field approach 

The Gutzwiller mean-field approach [245U244[[^ has been used in numerous papers 
devoted to the Bosc-Hubbard model. In its simplest version, it is based on the approx- 
imation of the many-body wave function by the product over single site contributions 

i^>=nE/»^i")- (23) 

i n—0 

where \n)i denotes the Fock state of n atoms in the i-th lattice site, rimax is a system 

size-independent cut off in the number of atoms per site, and fn^ corresponds to the 
amplitude of having n atoms in the i-th lattice site. The amplitudes are normalized 

To see what Gutzwiller approach predicts for quantum phase transitions, we focus 
now on the case of an homogeneous system having an integer number of particles 
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per site n. In this case, one obtains that for t/U < {t/U)g, where (t/U)g denotes a 

critical value, the Gutzwiller amplitudes are fn'^ = 6nn, so that the wave function 
becomes a single Fock state ([T6|) . Therefore, the Gutzwiller wave function exactly 
reproduces the system wave function when t/U — s- 0. It is also possible to argue 
that the differences between exact result and (|23l) arc negligible in the limit of a large 
lattice and t/U oo |246j . Therefore, the common expectation is that ([^5]) reasonably 
interpolates between superfluid and Mott insulator limits, which we discuss below 
pointing out the advantages and disadvantages of the Gutzwiller method. 

The critical point according to the Gutzwiller approach, {t/U)g, is located at — 
for n = 1 [246], where z is the number of nearest neighbors. Comparing this result 
to the more reliable findings from Sec. 12. 5[ one observes that the agreement improves 
with the system dimensionality: the Gutzwiller result is poor for a ID system (z = 2), 
and satisfactory for a 3D one (z = 6). 

The Gutzwiller method is quite straightforward in the numerical implementation. In 
the static case the amplitudes fn^ are real numbers and can be found by minimization 
of — ^J'Y^i niY^), where H is given by (fH)) . and ^ is a chemical potential used 

to enforce a desired number of atoms in the ground state. The minimization can be 
done in a standard way, e.g. with the conjugate gradient method |247| . and faces no 
problems as long as the system is homogeneous, or the external potential imposed on 
it is quite regular, e.g. harmonic. In the later calculation in an experimentally 

realistic 3D configuration consisting of 65^ lattice sites has been recently done |248| . 

The extension of the Gutzwiller approach to time-dependent problems is simple 
|139| . (Alternative dynamical mean field approach has been formulated in Refs. [249] .) 
Indeed, the stationarity of — ^^|^), where H is given by p^ . leads to the 

equation 

'1^" ^ = ^"(« - 1)/^^^ - ti>- V^/^'ii - i<i>.V^/l'2i, (24) 

where $, = E(ij)(^i) = E(,;j) E„ \/^/n-i/« ^ (the first sum goes over all j being 
nearest neighbors of i). Numerical integration of ()24p is straightforward. 

One should also appreciate the simplicity of Gutzwiller approach extensions to 
different systems, e.g., mixtures of bosonic gases [1391ll40j . Bose-Fermi mixtures |183l 
I210( l87j. etc. Particularly interesting are extended Hubbard models such as those 
involving dipolar interactions, where one expects the appearance of supcrsolid and 
checkerboard-like phases at low filling factors [1471I178[|250| . Despite all these nice 
features, there are also problems associated with the Gutzwiller ansatz. 

The picture that the Mott insulator corresponds to a single Fock state (fT6)) for the 
filling factor n is obviously incorrect - see a more detailed discussion of the Mott 
phase in Sec. 12.51 Another drawback of is that the correlation functions between 
different sites factorize into products of single site contributions, e.g., (\l/|6|6j|\l/) = 
for i 7^ j. As a result, there is lack of dependence of correlations 
on the distance between lattice sites. Also, formula (|23p does not correspond to a 
well defined number of particles. This problem can be solved by a proper projection 
of the wave function onto the subspace with fixed number of atoms |1391l244j . 
but the subsequent calculations become complicated. Finally, the Gutzwiller approach 
underestimates finite size effects: it predicts a "quantum phase transition" in systems 
of any size due to decoupling into the product of single site contributions (^5)1 . The 
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true quantum phase transition, however, requires a large system. 

The performance of the GutzwiUer ansatz can be, to a hmited extent, perturbatively 
improved |251| . These corrections significantly modify the GutzwiUer wave function for 
t/U smaller than {t/U)g. As a result, both the variance of an on-site atom occupation 

and the correlation functions {b\bj^i) become nonzero for t/U > 0, which is a progress 
with respect to the traditional GutzwiUer approach. 

Finally, we mention that the GutzwiUer approach can be supplemented by other 
mean-field-like calculations exploring the properties of Green functions |256ll2521l255l 
12531125411257] . They predict virtually the same transition points in different dimensions 
as the GutzwiUer method, but allow for determination of the excitation spectra and 
finite temperature calculations. 



3.5 Exact diagonalizations 

Exact diagonalizations of the Bose-Hubbard model can be done for small systems 
only. The problems result from the enormous size of the Hilbcrt space, given by 

where A'' and M stand for the number of atoms and the number of lattice sites, respec- 
tively. To illustrate predictions of (j25p we consider the simplest system undergoing a 
quantum phase transition in the thermodynamical limit, i.e., the M = N case. For in- 
stance, for TV = 8, 10, 12 one obtains nS{N,N) = 6435, 92378, 1352078 respectively. 
Moreover, one can easily verify that 

nS{N + l,N + 1) Ar»i^ 

ns{N, N) ^ ' 

which quantifies the fast increase of the Hilbert space size with system size. This 
shows that a standard diagonalization where all matrix elements are stored, and no 
symmetries are employed, can face problems already for N ^ M > %. Additionally, 
from the exponential increase of the Hilbert space size with the system size, any 
significant progress due to improvement of computer resources is unlikely. To over- 
come to some extend the above limitations, one can take into account the following. 
First, one can use numerical routines that store non-zero matrix elements only, e.g., 
ARPACK |258j . Then, diagonalization of A^ = 7\/ = 12 system faces no problems on 
a computer with about 1Gb of memory provided that one looks for a limited number 
of eigcnstates instead of the full spectrum. Second, one can cut the Hilbert space by 
restricting maximal site occupation to K atoms. Such a choice can be justified by a 
quadratic increase of the interaction energy with the site occupation number and is 
present in DMRG and Quantum Monte Carlo schemes. The size of the Hilbert space 
(see Appendix [B|) is then 



nSiN,M)W= f (-ir ('' + ^,/_7(^ + ^)) ), (26) 
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where [a;] stands for the largest integer number not greater than x. For K = N, 
([^S]) reduces to ([25]) . A typical choice leading to well converged results in a wide 
range of Bose-Hubbard model parameters is K = 4 |186[|242j . It results in the size 
of Hilbert space equal to 5475, 72403, 975338 for iV = 8, 10, 12, respectively. These 
numbers suggest that rather a slight progress can be achieved by cutting the Hilbert 
space size this way. The most powerful simplifications are possible when the problem 
under consideration has some symmetries, i.e., when there exists a set of operators 
{Oi} commuting with the Hamiltonian. Then, one can split a Hilbert space into 
subspaces composed of states with well defined eigenvalues of the Oi operators. In 
this way, the problem of diagonalization of the full Hamiltonian reduces to a few 
independent diagonalizations of smaller matrices. It allows for reduction of memory 
requirements and helps in getting a relatively large number of excited eigenstates. 
An example relevant for this review can be found in |110lllll| . where a system of 
cold fermions placed in a kagome lattice is considered. There exists a translational 
symmetry generated by an operator having S different eigenvalues. The Hilbert space 
splits into S subspaces, and diagonalization of the full Hamiltonian, say Lx L matrix, 
reduces to S independent diagonalizations of ^ L/ S x L/ S matrices. A more detailed 
description of exact diagonalization procedures for many-body quantum systems can 
be found in ^55] . 



3.6 Quantum Monte Carlo 

The Quantum Monte Carlo (QMC) method was successfully used for studies of the 
Bose-Hubbard model |1941l260ll261H2621ll95j . after publication of the seminal paper 
of Fisher et al. [30] ■ The computer resources and QMC algorithms allowed, at that 
time, to consider systems of the size of a few tens of lattice sites/atoms in one and 
two dimensional models. Nowadays, systems composed of 10'^ sites and atoms in any 
dimension can be routinely studied. It allows to consider 3D configurations, which 
are quite realistic from an experimental perspective [263,192J. The progress comes 
obviously from better computer resources and more efficient algorithms (see e.g., 
|2641ll09ll265p . Interestingly, efficient numerical codes for QMC simulations are now 
publicaly available |266| . 

The QMC approach allows to calculate different properties of a system being in 
equilibrium at finite temperature. During studies of the Bose-Hubbard model these 
temperatures can be chosen so low that the simulation describes essentially the zero- 
temperature physics of the system. The quantities usually calculated within the QMC 
approach are: the superfluid fraction, the chemical potential (//), the density of atoms 
(p), the variance of the on-site occupation, and the compressibility (k = dp/d^). 

Early QMC calculations [194, 260, 26L 262, 195j focused on determination of the 
phase diagram in the {p,/U,t/U) plane (Fig. [3]), and on the infiuence of disorder on 
it. In particular, they provided estimations of the position of the critical points, and 
they have verified critical behaviour of different quantities, e.g., scaling properties of 
compressibility and superfluid fraction, in the neighbourhood of the phase boundary, 
predicted in [3D] (see Sec. 12. 5p . 

Among the recent QMC studies of the Bose-Hubbard related to physics of cold 
atoms in optical lattices we would like to focus on investigations of harmonically 
trapped systems |2671I192[|2691I268] . Below we briefly review findings of these works, 
stressing the differences between the harmonically trapped model and the homoge- 
neous one. 
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The Bose- Hubbard Hamiltonian in the presence of the harmonic trap takes the form 
(fH)l with an additional term — J^i IJ'f"^^^'ni. The fj}°'^'^^ = — V^f is a local chemical 
potential, r.i is the distance from the trap centre, and V is the strength of a harmonic 
trap. The state diagram for such a problem in one dimension was calculated in [267] . In 
a generic situation one finds that there are plateaus, characterized by integer density 
of atoms and small atom number fluctuations, surrounded by regions of space where 
atom number fluctuations are large and density is non-integer. The former (latter) are 
identified as Mott insulator (superfluid) domains. Amazingly, images of these domain 
shells have been recently observed experimentally |270[|271| . It should be stressed, 
however, that there are significant differences between superfluid and Mott insulator 
"phases" in harmonically trapped and homogeneous models. 

The harmonically trapped system is gapless even in the presence of Mott do- 
mains j268| . It is also always globally compressible |267| . Both these properties are 
in a striking difference to the homogeneous case. Inhomogeneity of a trapped system 
suggests that the critical behaviour might be recovered in local quantities, e.g., the 
local compressibility ~ '^^^j ^j) / 8^°'^^'^ (see |192[|272j for a systematic discus- 

sion in both harmonic and quartic trapping potentials). This quantity, however, does 
not show the critical behaviour at the border between Mott and superfluid domains. 
All this leads to the conclusions that: i) in a harmonically trapped system, instead of 
a true quantum phase transition, there is rather a crossover, and ii) it will be very 
interesting to perform the experiments in homogeneous systems in a "box" poten- 
tial. Such potentials are currently being realized, for instance, by the group of M. 
Raizen [243] . 

Another difference between homogeneous and trapped models is observed in the 
visibility of interference patterns measured after releasing the atoms from the external 
potentials |273[|269| . While in a homogeneous system the visibility of interference 
fringes is a smooth monotonic function of t/U ratio, in the trapped model one finds 
both kinks (sudden changes of the slope) and non-monotonic behaviour caused by the 
presence of correlations between disconnected superfluid domains. 



3.7 Phase space methods 

Phase space methods have been introduced to quantum mechanics relatively early 
by Wigncr |275il274 ] and Moyal [276] . but their rapid development started in the 
60-ties with the applications to quantum optics by Glauber ^ and Sudarshan |277| . 
To a great extent, quantum optical studies of quasi-probability distributions, such as 
Wigner functions, or P- or Q- distributions, contributed enormously to our modern 
understanding of quantum noise [7l[Q . These methods usually map the density matrix 
of the considered system onto a function fulfilling a generalized Fokkcr-Planck type 
of equations, and try to replace this equation by systems of Langevin-like equations 
that can be simulated using classical Monte Carlo methods. The necessary condition 
is that the corresponding quasi-probability, or phase space quasi-distribution must be 
a reasonable probability measure, which not always is the case. Glauber-Sudarshan 
P-representation often is a highly singular distribution, and Wigncr function might 
take negative values. 

All these problems can be sometimes overcome, and several groups have started 
to use phase space methods to simulate many body problems. Pioneering work in 
this direction for bosonic gases has been done by P. Drummond's group |278| . These 
methods were extended to fermions in Ref . [279j . The authors study various signatures 
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of strongly correlated ultra-cold fermions in optical lattices, performing collective 
mode calculations, where a sharp decrease in collective mode frequency is predicted at 
the onset of the Mott metal-insulator transition. They have also looked at correlation 
functions at finite temperatures, using a new exact method that applies the stochastic 
gauge technique with a Gaussian operator basis. 

Somewhat similar approach has been developed by Y. Castin group [280j . This 
approach allowed, for instance, for precise determination of the fluctuations of the 
number of condensed atoms in an interacting Bose gas. It has been recently extended 
to describe lattice Hubbard problems [281] , and to Fermi gases |2821I283) . 

It is worth stressing that stochastic, or phase space methods become very efficient 
when the fluctuations become classical, i.e. at high, or at least moderate high tem- 
peratures, which for the BEG might mean 0.1-0.2rc. Truncated Wigner approach, or 
simulations of the Gross-Pitaevskii equations with random initial conditions, that 
mimic the initial thermal equilibrium are easy to implement and very accurate. 
Several groups have used such approach with great success. Hannover group used 
such method to describe phase fluctuations in quasi-lD BEC's [^Ml [^551 [^5511^57) . 
Barcelona-Hamburg collaboration applied these approach to study dynamics of spinor 
condensates [288]. 

The so-called classical field method was developed by K. Burnett j289j and K. 
Rzfizewski [290] groups. Both groups put some emphasis on understanding the concept 
of temperature in the microcanonical ensemble, whereas the latter recently realized 
the modeling of the decay of unstable vortex states. 



3.8 ID methods 

In this section we give a very short guide of ID methods, illustrating the subsequent 
steps from Bose-Hubbard model to Luttinger liquid theory. We start by analyzing the 
Bose-Hubbard Hamiltonian in ID (Eq. (|12|) ). i.e. a chain of sites with open ends and 
we perform then the following steps: 

a. Hard core bosons. In this limit (i.e. when U 3> /i), as descibed in Section [^?^ 
the Bose Hubbard hamiltonian (Eq. ([12])) reduces to the XY model in a transverse 
field (Eq. 

b. Jordan- Wigner transformation. Jordan- Wigner transformation is a way to 
"fermionize" the, otherwise, bosonic system |154] . One defines: 

f^ = l['yz,J<y^, (27) 

j<i 

fl-Il^^A^ (28) 

so that fffi = <j\(Ji = (1 + (T^,i)/2. It is easy to check that such nonlocal operators 
are fermionic, i.e. 

{/.,/,} = {///]} = 0, {/4} = %. (29) 
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The inverse relations are also simple: 



n(2//.A - 1)/, 



(30) 



n(2///. - 1)//. 



(31) 



With this transformation the Hamiltonian becomes the one of "spinless" fcrmions: 



and the model becomes exactly solvable. Had we started from an extended Hubbard 
model with, say, nearest neighbour interactions of the form V '^^a\a\_^-^ai+iai, we 

would end up with "spinless" interacting fcrmions with the term ft fl+ifi+ifi- 

Jordan- Wigncr transformation is essentially a ID transformation. There were very 
many attempts of generalizing it to higher dimension (for a review see |235| . for a 
discussion in the context of applications for atoms see [183j ). Only recently, new ideas 
related to the concept of PEPS, allowed Verstraete and Cirac to propose an interesting 
efficient generalization of the Jordan- Wigner transformation to 2D _291J. 

c. Bosonisation, and Luttinger liquid theory. Let us consider now the interact- 
ing "spinless" Fermi gas, and try to formulate the low energy effective theory for this 
model. We assume that fcrmions do not fill the lowest band, so that in the absence 
of interactions, the Fermi level is somewhere in the middle of the band. We expect 
that at low temperatures interesting physics will occur close to the Fermi energy. It 
is thus reasonable to linearize the fermionic dispersion relation at the Fermi energy, 
e{k) = e{±kp)±cp{k±kF), with the " sound velocity" cp = de{k)/dk\kp- In ID there 
arc two values of momenta ifc^? where the Fermi energy is reached: one corresponds 
to left, and another to right going fcrmions. Note that all of the interesting low energy 
physics happens close to the Fermi surface (points). That implies that the states with 
momenta far from ikp will practically never participate in any relevant physical pro- 
cess, and will remain deep in the filled Fermi sea. It is thus reasonable treat the left 
and right going fcrmions close to ±kp as independent and introduce two fermionic 
species, described in the momentum representation by L{k),L^{k), and i?(fc), R^{k), 
respectively; k is here momentum relative to zkkp, so that it attains values from — oo 
to -|-oo. The Hamiltonian, becomes 




(32) 



H = +^[-CF{k)L\k)L{k) + CF{k)R\k)R{k)] 



k 




(33) 



This is a Luttinger model [292[|46j . One introduces now the operators 




(34) 



k 
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pd^q) = Y.^\k)L{k + q) = p{{q), (35) 

k 

Pfl(<7) = ^i?^(fc + g)i?(fc), (36) 

k 

PR{~q)^Y.^^^^)^^^ + ^)^Pl^^)^ (37) 
k 

and observes that [p_r(— g), Pi?(+<z')] ^ X^A-F-g' R^{k + q — q')R{k). In the Luttinger- 
Tomonaga approximation [292j one replaces now these and similar commutators by 
their averaged values over the ideal Fermi sea. The result is 

[pR{~q).PR{+q')] = ^5{q-q'l (38) 

implying that for g > 0, we can introduce 6j = p R{q) 27^ / qL ^ hq = p r{~ q) 2ti / qL, 
which may be regarded as bosonic annihilation and creation operators, respec- 
tively, since bq,b^^, = 6{q — q'). Similar construction is done for the operators 
c^_q = pL[—q)yj2'K/qL, and c_g = pL{q)yj2'K / qL. These operators have the inter- 
pretation that h]j takes a particle from state k and puts it into k + q. It creates thus 

fermion-hole pairs when k < kp and k + q > kp- Analogously, the operator cl^ takes 
a fermion from an unoccupied state —kp < k + p to the unoccupied state k < —kp. 
The whole Hamiltonian becomes 

H^Y. 9Cf[6^?^, + clqc.q] + ^ Y.^V{q)clqbl + 1/ (?) * 6,c_,] . (39) 

q>0 q>Q 

The problem has been thus reduced to an exactly solvable system of interacting har- 
monic oscillators, describing linear ID hydrodynamics of a, so called, Luttinger liquid. 
Bosonisation theory |46] assures that such description can be found for most of the ID 
bosonic and fermionic systems with local interactions. For the very recent and very 
complete description of the interacting Bose gases in quasi-one dimensional optical 
lattices see the review article of Cazalilla et al. |293j . 



3.9 Bethe ansatz 

Bethe ansatz is an analytical method to find exact eigenstates and eigenvalues of 
some strongly correlated one-dimensional models (although some times it has to be 
complemented with numerical analysis). By exploiting appropriately the symmetries 
involved in the Hamiltonian it is possible to diagonalise the Hamiltonian exactly. 
However, its power relies in the fact that Bethe ansatz characterises all eigenstates 
by a set of quantum numbers which enumerate the states, according to their physical 
properties. This method was originally developed by Bethe in 1931 [294j to solve a ID 
array of electrons with uniform next neighbour interactions, i.e., spin 1/2 Heisenberg 
model: 

N 

i^ = -J^Ti-ai+i, (40) 
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where a = {(Jx,cry,az) are the standard PauU matrices and N denotes the number 
of sites. The parametrisation of the eigenvectors, i.e. the Bethe ansatz has become a 
fundamental tool with which many other ID quantum systems have been shown to 
be solvable. From the eigenvectors one can compute easily the quantities of interest 
by calculating the expectation values of the desired operators. In what follows, we 
summarize the original formulation of Bethe following the pedagogical work of Kar- 
bach et. al. |295H296j . Since, the Heisenberg Hamiltonian commutes with the total 
spin along the z-dircction = J^i i [-^^i ^t] = 0, then: 



h=-jj: 



(41) 



Thus, eigenstates of cr^ are also eigenstates of H . Bethe ansatz is, as discussed below, 
a basis transformation. 

We denote by |0) = | t) the eigenstate of ct^ with value +1/2 and by |1) = | |) the 
eigenstate with eigenvalue -1/2. Sorting the basis according to the quantum number 
N/2—r, where r is the number of flipped spins, is all that is needed to block diagonalise 
the Hamiltonian. The block r — corresponds to a single state with all spins up 
which we call ferromagnetic: H\F) = H\ ^ ... f) = —JN/4. The case r = 1 (one spin 
down) has N invariant vectors, which are labeled by the position of the flipped spin: 
\n) = (T~\F). To diagonalise this block, which has size NxN,we take into account that 
the Hamiltonian possesses translational symmetry. Therefore a translational invariant 
basis can be constructed in the subspace with r = 1, 

1 ^ 

l'/') = 7f E^^'"l")' (42) 

^ n— 1 

with wave numbers k = 2Trm/N, m = 0, 1, ...iV — 1. The vectors {ip) correspond to a 
complete spin alignment of the ferromagnetic ground state \F) which is periodically 
disturbed by a spin wave with wavelength A = 27r/fc. These states are called magnons, 
or spin-waves. 

The subspace r = 2 cannot be solved by applying any further symmetry and here 
is where the full power of the Bethe ansatz appears. Solutions (|42p . for r ^ 1, can be 
also obtained starting from the following ansatz 

N 

|^) = ^a(n)|n) (43) 

ri=l 

and finding the values of the coefficients a{n) that satisfy the eigenvalue equation 
Hltp) = E\ip) using periodic boundary conditions. If one applies the same procedure 
to the subspace r = 2, i.e. one looks for the coefficients a(ni,n2) that determine the 
eigenstates, one finds: 

N 

IV') = X] a(ni,n2)|ni,7i2). (44) 

I<rii<ri2<-/V 

There are N{N — l)/2 eigenstates. These are characterised by pairs of Bethe (integer) 
quantum numbers Ai, A2, which define the "momenta" via the relations Nki = 2ttXi + 
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9, Nk2 = 27rA2 - 6, such that 



(45) 



whereas the angles 6*12 = -6I21 = 9 fulfill 2cot(6l/2) = cot(fci/2) - cot(A:2/2). De- 
pending on the relative values of the pair Ai, A2, different type of solutions exist. In 
general, the Bethe ansatz for an unrestricted number r of flipped spins reads: 



N 

E 

l<)ii<...<rv<JV 



a{ni...nr)\ni...nr) 



(46) 



and the eigenvectors span a subspace of dimension N\/{N — r)\r\. The coefficients 
a{ni...nr) have an analogous form to and a generalised solution of the resulting 
transcendental equations can be easily found numerically. 

It is also worth to mention that the Bethe ansatz can be applied to the one di- 
mensional Bose-Hubbard model - see [187] for a complete discussion. In this case, 
however, the method faces fundamental problems since one has to assume that the 
maximal occupation of each lattice site is not larger than 2. As a result, the ansatz 
is never exact though it gives some valuable and correct, predictions, e.g., vanishing 
superfluid fraction in the Mott phase. 



3.10 A quantum information approach to strongly correlated systems 

The density matrix renormalisation group (DMRG) [2311IT7] is a variational method 
that has had an enormous success in describing ground states of some strongly in- 
teracting ID systems with rather modest computational effort. The underlying phi- 
losophy of all DMRG oriented algorithms is that many body systems can be treated 
almost "exactly" if one is able to truncate the full Hilbert space by removing the 
degrees of freedom that are not involved neither in the ground state, nor in the dy- 
namical evolution of the system. The difficulty and glory of the method relies on how 
reliable the truncation is done. Very recently |297[|298ll300[|299[l301| . quantum infor- 
mation theory has provided a new perspective on the following questions: (i) how to 
perform an efficient truncation of the Hilbert space, (ii) which quantum systems can 
be efficiently simulated, (iii) how to simulate dynamical evolutions of strongly corre- 
lated systems, (iv) how and when DMRG-oriented methods can be implemented to 
investigate ground states of 2D and 3D systems, (v) how classical concepts like corre- 
lation length, which diverge on the critical points is linked to entanglement |302il30lj , 
etc. In general, this approach is shedding new light in our understanding of complex 
many body physics. Anders et al. have recently proposed an approximate method to 
calculate ground states of quantum many-body systems, based on quantum informa- 
tion concepts of "weighted graph states" [303j . 

Real space renormalisation methods are iterative methods to describe accurately 
mainly ID many body systems. The starting point is a small ID chain system which 
can be exactly diagonalised. Then, the size of the chain is iteratively increased by 
adding a new lattice site at each step. The Hamiltonian that coupled the new site to 
the chain is renormalized by disregarding all physically irrelevant couplings and the 
resulting Hamiltonian diagonalised. 

There are, at least, 3 different approaches to perform the truncation [i7] : 
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• By optimizing expectation values |304j . In this approach the value of some bounded 
operator (energy, magnetization or density) determines which are the most relevant 
states contributing to it. 

• By optimizing the wave function [231| . The truncation is done such that the renor- 
malizcd wave function minimizes the Schmidt norm distance to the exact one, i.e. 

Ill*) -|*>Z5Mi?G|P^0 

• By optimizing entanglement properties. It was shown by G. Vidal |297| that the 
efficiency in simulating many body systems is directly related to their entanglement 
behaviour. Efficient simulation is possible if the entanglement of a subsystem with 
respect to the whole is bounded, or grows at most logarithmically with its size. Also, 
if entanglement grows linearly in time and block size, simulation of time evolution 
may be not efficient even in some non-critical ID systems |299| . 



3.10.1 Vidal 's algorithm. In the following we briefly review the algorithm 
of G. Vidal [M] (also called Time Evolving Block Decimation (TEBD)) to cal- 
culate ground states and dynamics using the Schmidt decomposition to truncate 
the Hilbert space. Implementing operations on individual sites and on 2 neighbor- 
ing sites requires only local updating of the expansion, rending the calculations 
of ground states, expectation values, and 2-body correlations efficient. We also 
comment how the expansion is linked to the Matrix Product states (MPS) ansatz 
used in DMRG, and why MPS describe so well ground states of quantum spin 
systems |300j . 

Assume that our system consists of a set of atoms loaded in a ID optical lattice 
with M-sites, whose Hamiltonian is well described by a Bose-Hubbard Hamiltonian 
(see Eq. (HI])) and we are seeking for the ground state 1^") of the system. We choose 
as ansatz the most general description of a state in such Hilbert space which is 
given by 

d d 

l*) = E- E C^^i....m,Ni)®...Na/) (47) 

10=0 iM=0 

where {|jj)} denotes an orthonormal basis, M is the number of sites, TV is the 
number of atoms, and d indicates the maximal occupation number of site j (max- 
imal (a|aj)). Note that by taking d — N, wc can be sure that Eq. (|T7|) gives an 
exact representation of the state. Let us now split the ID chain into two blocks: 
block A, which consists of a single site, and block B, consisting of the remaining 
M — 1 sites. For this (or any other bipartite) splitting {A, B), there exist always a 

bi-orthonormal basis denoted by |<&[^') and | respectively such that: 

I*> = Ea„I*L^i>I*L''i>, (48) 

a 

where XA ^ d, i.e. the number of coefficients (Schmidt coefficients) of the expansion 
is bounded by the dimension of the smallest of the two subspaccs of the partition. 
This simply indicates that the number of degrees of freedom that can be entangled 
between A and B, is the maximum number of degrees of freedom of the smaller 
subsystem. If there is a single coefficient different from zero in the expansion, then 
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the two parties are in a product state, i.e. not entangled. The coefBcients of the 
expansion are unique (up to the degeneracy of the reduced density matrices) , real 

and Ea^a = 1- Notice also that pA = rrB(|*)(*|) = Ea I^L"^') (^L"^' I and 

= TrA(|^')(*|) = Ea-^al^a^'X^a^'l- The Schmidt decomposition is nothing 
else than the usual Single Value Decomposition (SVD) of the coefficient matrix C 
corresponding to the decomposition of \^) in an arbitrary orthogonal basis \i)A 
(b>B): 

\^)^Y.^^,,V)^\^)B■ (49) 

Ci,j is a d X d^^^^ matrix. Applying the Single Value Decomposition to C one 
obtains that C = U DV , where U (V) is a unitary matrix of dimensions d x d 
j-^M-i ^ f^^-i) and D is a diagonal matrix of dimensions d x d^^^^. Therefore, 

i,j a I I 

= EAa|4^I)|0L^I). (50) 

An iterative treatment of the above procedure (i.e. splitting now part B into a block 
of 1 lattice site and block of AI — 2 sites and finding its Schmidt decomposition, 
and again for the block of M — 3 sites) leads to the following expression: 

(7- . - V rl^l'iAl^lrl^l^^ \[2]p[3]i3 p[M]i„ 
ai-.-CcM-l 

The tensors F take into account the correlations arising from the splitting, and 
are straightforwardly obtained by performing at each bipartite partition the cor- 
responding SVD. Notice that except for the first and the last lattice site, for all 
the other lattice sites, the corresponding tensor T depends on two a indices, cor- 
responding to the two involved partitions, and takes into account the correlations 
with the block at its left, and the block at its right. It can be shown that for 
translationally invariant ID systems with short range interactions, the coefficients 
of the Schmidt decomposition behave as Aq ^ cxp(— a), and the truncation of the 
Hilbert space is performed by removing all small coefficients. 

By grouping the F tensors with the coefficients A the expression of the ansatz takes 
the more compact form: 

d 

E Tr[AW-'^..^M'^-]|*i,...*M), (52) 

ii,...iM=0 

where A'^''^^'^ . . .AW''' are Dk x Dk complex matrices with Dk+i — < D, and 
d is the occupation number on site i. Operations on a single site i involve only 
updating the value of ^[''I'*'= , which basically implies of the order of operations, 
^ {7^[fe],ifc^ In a similar way, operations involving two neighbouring sites 
correspond to updating the matrices of the corresponding sites and the correlations 
between them and again can be efficiently performed. Thus, the ground state of 




10 9 8 7 

/'O'N /a^^ /'O'N /'O'v /cy\ 

( 04 K) OH hoi K) OH fo ! -^gOj hO tO Oj tO (p ,1 



11 



12 13 14 15 



■(O OH (O OH (O ) qO) I.CUP! 



20 19 18 17 



16 



01 10 Ol K> Ol K) 04- 



w '^J 



(^^ (O*^ (O*^ 0^ 

Figure 5. Schematic representation of a) MPS for a ID system, b) "inefficient" MPS in 2D, and c) 
PEPS for a 2D system. Each pair of nearest neighbour sites is connected via maximally entangled 
state from an auxiliary Hilbert space. In a) and b) matrices, A projects onto the physical space. In 
c) matrices have to be replaced by 4th rank tensors in the bulk (from [306]). 

the Bose-Hubbard Hamiltonian, which includes only next neighbour interactions 
plus on site collisions can be calculated variationally, whereas time evolutions can 
be simulated by means of a Trotter expansion. (For the recent work of dynamical 
response of BH model at the MI-SF transition, see ISOS] .) 
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3.10.2 Matrix product states. Already Ostlund and Rommer [30711308] re- 
alized that DMRG leads to an special type of ansatz known as Matrix Product 
States p09l[3T0] . Matrix Product States (MPS) and their gencrahzat ion to higher 
dimensions, Projected entangled-pair states (PEPS) (see Fig. (5]) corresponds to 
translationally invariant states, whose general expression is of the form [308] : 

d 

|*)= Tr[^W'*i...^M'*«]|zi,...ZM), (53) 

ii,...ii/ = l 

which describes a chain of AI spins of dimension d. The ^I'^l''-- are complex matri- 
ces with dimension < D x D matrices. The reason why this type of states (MPS) 
can reproduce very well the ground states of quantum spin systems with nearest 
neighbor interactions, is because the ground states (that have also the transla- 
tional invariance of the Hamiltonian) and low energy excitations are completely 
determined by their two-body reduced density operators as shown by Verstraete 
and Cirac [311] . An intuitive but handwaving argument to understand the reason 
why MPS parametrize ground states of gapped quantum spin systems, comes from 
Hastings theorem [312] . which states that for gapped systems all correlations de- 
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cay exponentially. This means that since {OaOb) — {Oa){Ob) — cxp(—^AS /'Ccorr); 
blocks A and _B, separated by Iab, can be described by a product state, when Iab 
is much larger than the correlation length. 

MPS and PEPS have already proven to be an enormously useful tool. In particular, 
they can be easily extended to calculate time evolution. They have also been used 
to characterize quantum phase transitions in spin chains |300] . simulate infinite 
systems |313| . and establish a closer relation between criticality, area law, and 
block entanglement [314] . 



3.11 Fermi and Fermi-Bose Hubbard models 

Most of the methods that are used for Bose-Hubbard models (some of which were 
discussed in the previous subsections of this section) can be carried over to fermionic 
systems, although typically, treatment of fcrmions is more difficult. For example, 
mean field BCS method, valid for weakly interacting fermions, requires much more 
effort for the trapped gases than the bosonic Bogoliubov-de Gennes approach. 
Quantum Monte Carlo methods (QMC) suffer notoriously the, so-called, "sign 
problem" (negative probabilities). Gutzwiller ansatz in dimensions higher than 1, 
does not take the fermionic anticommutation relations into account properly, due 
to non-existence of standard Jordan- Wigncr transformation (for a discussion, see 
for instance |183j ). 

When we deal with Fermi-Fcrmi, or Fermi-Bose mixtures, it is always useful to 
study some limiting cases, that moreover, can be realized in experiments. One limit 
would be the weakly interacting regime, where Hartree-Fock, BCS or Landau Fermi 
liquid theory could be applied. In the strong coupling limit of Mott insulating states, 
a very useful method is to construct an effective Hamiltonian. Such Hamiltonian 
will often correspond to a spin model, if wc deal with commensurate filling factors 
(as we discussed in the Introduction). If the Mott states involve states differing by 
odd number of fermions (for instance presence or absence of a fermion), then the 
effective model will necessarily involve composite fermions. 

Fermionic Hubbard models. Perhaps the most well-known are the effective mod- 
els obtained in the strong coupling limit of the electronic Hubbard model for spin 
1/2 fermions, described by the Hamiltonian, 

H = -tY. iflfi- + ^■^■) + ^Y. fhflhf^^ - E /t/- (54) 

In the limit of large repulsive interactions, U ^ t, there can be one, or zero 
particles at a site. In this limit the model is termed in the literature "i — J" model 
(cf. [5D]). When the system is at full filling, i.e. there is 1 fermion at each site 
(in the Mott limit), the effective Hamiltonian reduces to that of Heiscnberg quan- 
tum antiferromagnet with the magnetization fixed by the corresponding numbers 
of spin-up and spin-down fermions (or external "magnetic" field). In particular, 
if these numbers are equal, the system has a zero net magnetization (vanishing 
"magnetic" field). Allowing for certain number of unoccupied sites, is termed dop- 
ing. Doped Mott insulators provide a key to understand high Tc superconductiv- 
ity of cuprates (for a recent review, see |198| : for the recent contribution in the 
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atomic context, see [199j ). P. Torma's group has investigated in recent years, the 
role of spin-density imbalance, both in trapped [200] , and superfluid lattice Fermi 
gases |201[|202j . Asymmetric Hubbard model with light and heavy fermions has 
been studied by Ziegler's group, who has, in particular, shown quantum phases 
with self-induced disorder ppgll^l^ . 

Hubbard models for particles with higher spin, as well as Fermi-Fcrmi mixtures 
can be realized with ultracold atoms. There has been thus considerable interest 
in studying such systems (for a review of SU{N) symmetric Hubbard models of 
fermions and bosons see |206] ). In Ref. |122| effective spin Hamiltonians for spin 
3/2 ultracold spinor gases, described by the spinor version of the Hamiltonian (|54p 
have been derived and investigated. In ID the same system has been studied using 
the bosonisation approach [207j . and exact Bethe ansatz method [208] . 

Fermi-Bose Hubbard models. The spectacular advances in loading atomic sam- 
ples in optical lattices have allowed for the realization of systems which are well 
described by a Fermi Hubbard model, as well as mixtures of fermions and bosons 
(FB) described by Fermi-Bose Hubbard models. In the latter case, and in the limit 
of strong atom- atom interactions, such systems can be described in terms of com- 
posite fermions, consisting of a bare fermion, or a fermion paired with one boson 
(bosonic hole), or two bosons (bosonic holes), etc. [209] . The quantum phase di- 
agram displayed by those systems is amazingly rich and complex. The physics of 
Fermi-Bose mixtures in this regime has been studied recently in [1671ll83ll210| . 
using perturbation theory up to second order to derive an effective Hamiltonian. 
There are a number of recent studies of FB mixtures in optical lattices Q, and 
also of strongly correlated FB mixtures in traps |1661l218j . In particular, the va- 
lidity of the effective Hamiltonian for fermionic composites in ID was studied us- 
ing exact diagonalization and Density Matrix Renormalization Group method in 
Ref. |219| . whereas trap effects were analyzed in Ref. |i220i . Very recently, exact 
Bethe ansatz solution for the special case of FB Hubbard model has been found 
and intensively studied |221[|222) . Also, intensive studies of quantum phases in the 
multiband Fermi-Hubbard model (that models a Fermi gas subject to Feshbach 
resonance) have been undertaken |223ll224 [ |225j . Recently, particular interest has 
been stimulated by the experiments of the Hamburg and ETH groups [227[|226| . 
in which presence of fermions attracting bosons was seen to decrease the lattice 
potential threshold for MI state. This is quite paradoxical, since DMRG calcula- 
tions predict the opposite effect j228| The physics of disordered FB mixtures will 
be reported in subsection 14.71 

Here, we focus on the homogeneous case and explain how to derive the corre- 
sponding effective Hamiltonian. We consider a homogeneous mixture of ultracold 
bosons (b) and spinless (or spin-polarized) fermions (f), for example ^Li-^Li or 
®^Rb-^°K, trapped in an optical lattice. As for the bosonic case, in the tight- 
binding regime, it is convenient to project the wave functions onto the Wannier 
basis of the fundamental Bloch band, corresponding to wave functions well local- 



^The studies of trapped FB gases concerned in particular FB phase separation I211II2121 . the phase 
diagram |213j , novel types of collective modes |212I214] , Fermi-Fermi interactions mediated by bosons 
121211215] , the collapse of the Fermi cloud in the presence of attractive FB interactions [216) , or the 
effects characteristic for the ID FB mixtures |217l . 
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izcd in each lattiee site [ 15511229] . This leads to the Fermi-Bosc Hubbard (FBH) 
Haniiltonian [Sni.il54..29.,gB] : 



U 

—ni{ni - 1) + Vriimi - Ui - fi^rrii 

(55) 

where 6^, bi, j} and fi arc bosonic and fermionic creation- annihilation operators 
of a particle in the i-th localized Wannier state of the fundamental band, and = 
b\bi, nrii = fjfi are the corresponding on-site number operators. The FBH model 
describes: (i) nearest neighbor (n.n.) boson (fermion) hopping, with an associated 
negative energy, —tb (— t/); (ii) on-site boson-boson interactions with an energy 
U, which we will assume to be repulsive (i.e. positive), (iii) on-site boson-fermion 
interactions with an energy V, which is positive (negative) for repulsive (attractive) 
interactions and finally (iv) on-site energy due to chemical potentials in grand 
canonical description. Following [230.,87>il551l210| . let us review how to derive an 
effective Hamiltonian to second (or higher) order int = tf = tb (we assume here the 
same tunneling for bosons and fermions for the sake of simplicity). Generalization 
to the case tf tb is straightforward. We follow here the notation of Ref. |167| . 
In the limit of vanishing tunneling (t — 0) with finite repulsive boson-boson inter- 
action U, and in the absence of interactions between bosons and fermions {V = 0), 
the bosons are in a Mott insulator (MI) phase with exactly n = + 1 bosons 
per site, where fi'' = fJ^ /U and [x] denotes the integer part of x. In contrast, 
the fermions can be in any set of Wannier states, since for vanishing tunneling, 
the energy is independent of their configuration. The situation changes when the 
interparticle interactions between bosons and fermions, V ^ are turned on. In the 
following, we define a — V/U , and consider the case of bosonic MI phase with fi 
bosons per site. We also assume that the temperature is small enough (T < V) in 
order to neglect the particle-hole excitations. It follows that tunneling of a fermion 
is necessarily accompanied by the tunneling of — s bosons (if s < 0) or opposed- 
tunneling of s bosons (if s > 0). The dynamics of the Fermi-Bosc mixture can thus 
be regarded as the one of composite fermions made of one fermion plus — s bosons 
(if s < 0) or one fermion plus s bosonic holes (if s > 0). The annihilation operators 
of the composite fermions are [167] : 



[fl — S}' / 4.\^ 

Fi ~ y - — r-p^ [blj fl for s bosonic holes (56) 

= 1/7"^ TT (^*) /» bosons. (57) 

y (n - s)! 

These operators are fermionic in the sub-Hilbert space generated by |n — ms,m) 
with m = 0, 1 in each lattice site. Note that within the picture of fermionic com- 
posites, the vacuum state corresponds to MI phase with n boson per site. At this 
point, different composite fermions appear, depending on the values of a, n and 
p,'', as detailed in Fig.[6l 

Because all sites are equivalent for the fermions, the ground state is highly 
[iV!/iV/!(A^ — A^/)!]-degenerated, where TV denotes the total number of atoms and 
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Figure 6. Quantum phase diagrams of Fermi-Bose mixtures in an homogeneous optical lattice as 
functions of fi^ and a = V/U , for pf = 0.4 and t/U = 0.02. Roman numbers denote the total 
number of particles that form the composite and a bar means that the composite is formed by 
bosonic holes rather than bosons, (a) Diagram of composites where the filled small (blue) dots 
symbolize fermions, large (red) dots symbolize bosons and empty (red) dots, bosonic holes. The 
subindex A (R) indicates attractive (repulsive) composites interactions, (b) Detailed quantum 
phase diagram of fermionic composites. The subindices denote here different phases: DW (density 
wave), FL (fcrmi liquid), SF (superfluid) and FD (fermionic domains). The strongly correlated 
phases for small but finite t are surrounded by characteristic lobes 11831 . beyond which bosons 
become superfluid. Therefore, there are thin regions of bosonic superfluid between the various 

composite phases (from 11671 ). 

Nf the number of fermions. Hence, the manifold of ground states is strongly cou- 
pled by fermion or boson tunneling. We assume now that the tunneling rate t is 
small but finite. In the homogeneous case, one uses standard second order pertur- 
bation theory, to derive an effective Hamiltonian [60j for the fermionic composites 
(for the disordered case, see Appendix [K\ |230[(57] ) : 

Heff = -deff {f}f, + h.c.) + K,ff M,M, -KffJ2 (58) 

where Mi = PjPi and /le// the chemical potential, which value is fixed by 
the total number of composite fermions. The nearest neighbor hopping for the 
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composites is described by —deff and the nearest neighbor composite-composite 
interactions is given by Keff, which may be repulsive (> 0) or attractive (< 0). This 
effective model is equivalent to that of spinless interacting fcrmions. The interaction 
coefScicnt K^f f originates from 2nd order terms in perturbative theory and can be 
written in the general form |183| : 



Keff = — 



{2n - s){n + I) - s{n - s) ~ > 



1 + s — a 1 — s + q; sa 

(59) 

This expression is valid in all the cases, but when s = 0, the last term (1/sa) 
should not be taken into account. c?e// originates from (|s| + l)-th order terms 
in perturbative theory and thus presents different forms in different regions of the 
phase diagram of Fig. [51 For instance in region I, dg// = t, in region // deff = 21^ /V 
and in region II, deff ~ At^ /\V\. 

The physics of the system is determined by the ratio K^f f / d^f f and the sign of 
Keff- In Fig. [6lja), the subindex A/R denotes attractive {K^ff > 0) / repulsive 
{Keff < 0) composites interactions. Fig. [5]Jb) shows the quantum phase diagram 
of composites for fermionic filling factor pf = 0.4 and tunneling t/U = 0.02. For 
large values t/U, a transition to the SF state takes place (for the corresponding 
3D lobes over the different regions in Fig. [SI see Ref. |183j ). 



4 Disordered ultracold atomic gases 



4.1 Introduction 

In 1958, Anderson [66j reported for first time the quantum localization phe- 
nomenon. He pointed out that the extended wave functions of electrons (Bloch 
waves) in a crystal, for strong enough disorder, become localized with an exponen- 
tially decaying envelope {ipl ~ exp{\f — rll/S,), where ^ is the localization length. 
This occurs at the single particle level by coherent back scattering from random 
impurities. In one dimension, it has been rigorously proven that infinitesimally 
small disorder leads to exponential localization of all the eigenfunctions by re- 
peated backscattering |315[I316 [ [67] and it is known that the localization length is 
of the order of the backscattering mean free path. Ref. [317] proposed the concept 
of mobility edge, which separates the localized states from the extended ones. The 
scaling theory has been one of the basic tools to deepen the understanding of the 
Anderson localization phenomenon by considering the conductance as a function 
of the system size or of other scale variables pT51[5Tgi[5^l32 1113221(5^ . In ID the 
localization length is a function of the ratio between the potential and the kinetic 
(tunneling) energies of the eigenstate and the disorder strength. For the case of dis- 
crete systems with constant tunneling rates and local disorder distributed according 
to a Lorentzian distribution (Lloyd's model, cf. j324j ) the exact expression for the 
localization length is known. In general, an exact relation between the density of 
states and the range of localization in ID has been provided by Thouless }319j . In 
2D, following the scaling theory, it is believed that localization occurs also for arbi- 
trarily small disorder, but its character interpolates smoothly between algebraic for 
weak disorder, and exponential for strong disorder. There are, however, no rigorous 
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arguments to support this belief, and several controversies aroused about this sub- 
ject over the years. In 3D, scaling theory predicts a critical value of disorder, above 
which every eigenfunction exponentially localizes, and this fact has found strong 
evidence in numerical simulations. 

Since the discovery of Anderson [66| . the field of quantum disordered systems 
has been a very active research field in condensed matter physics (for reviews 
see |325[l3261l3271l328j ). Disordered systems appear in some of the most challenging 
open questions concerning many body systems (see for instance j3291l3301l331p . In 
particular, quenched (i.e., frozen on the typical time scale of the considered systems) 
disorder determines the physics of a large variety of phenomena: transport, con- 
ductivity, localization effects and metal-insulator transition (cf. [325,326,327,355]), 
spin glasses (cf. f551l332lll54p . neural networks (cf. |333| ). percolation ^334ji335j . high 
Tc superconductivity (cf. [50]), or quantum chaos (cf. |324| ). 

The theoretical studies of disordered systems imply severe difficulties. One of the 
main ones is that in order to characterize the systems independently of the partic- 
ular disorder realisations, one needs to average over the disorder. Such quenched 
(frozen over the typical time scales of the system) average of physical quantities, 
such as for instance free energy, requires usually the use of special methods, such as 
the replica trick (cf. [53]) or supersymmetry method [336j . Interestingly, quantum 
information approach might help to calculate averages using quantum superposi- 
tions, see |337j . Regarding numerical approaches, the characterization of the system 
demands cither simulations of very large samples to achieve "self-averaging" or nu- 
merous repetitions of simulations of small samples with different configurations of 
the disorder. Obviously, this difficulty is particularly important for quantum dis- 
ordered systems. Another difficulty that arises in the study of disordered systems 
is the possible existence of a large number of low energy excitations. Similar limi- 
tations arise in frustrated systems (cf. |107| )fsec section (5]). Last but not least, the 
interplay of the disorder and the interactions has been and still is a challenging 
problem, as it is reported in the following subsections in the case of bosons and 
fermions. 

Disorder effects are present in many condensed matter systems, and in particular 
in electronic systems. Nevertheless, these systems do not provide the best possible 
scenario to test the theoretical predictions mainly because of two factors: (i) the 
disorder is not controllable, but fixed by the specific realization of the sample; and 
(ii) the electrons interact via the long-range Coulomb interaction, which is fixed by 
Nature. It is thus desirable to ask whether atomic, molecular physics and quantum 
optics may help to understand quantum disordered systems, for instance by real- 
izing in experiments disordered models using cold atoms in optical lattices. The 
advantages of atomic disordered systems are clear: (i) the control of the generation 
of random potentials that induce disorder in the system (see subsections 14.31 and 
14. 5p : (ii) the presence of interactions that can also be controlled by means of Fesh- 
bach resonances; (iii) the type of quantum statistics, that can be chosen by using 
ultracold bosons or fermions, and (iv) the full control of the trapping potentials 
and therefore the effective dimensionality of the system. Obviously, standard con- 
trol parameters, such as for example the temperature of the system, can also be 
controlled, but in these aspects atomic system do not differ much form the con- 
densed matter one. Of course, there are also disadvantages: atomic systems are 
"small". Typical experiments with ultracold lattice gases involve lattices of sizes 
up to 100 X 100 X 100 and up to few times iV = 10^ atoms. 
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4.2 Disordered interacting bosonic lattice models in condensed matter 

The theory of disordered interacting bosons is complex and there are essentially 
no exact solutions, not even in one-dimension. In this section we review different 
numerical and approximate results existing in the condensed matter literature. 

A system of bosons in a lattice with short-ranged repulsive interactions in the 
presence of external random potentials was considered in |30| where, by treating the 
tunneling as a perturbation, the phase diagram of the system was worked out. The 
three possible ground states predicted were: (i) an incompressible Mott insulator 
with a gap for particle-hole excitations; (ii) a gapless Bose-glass insulator with finite 
compressibility, exponentially decaying supcrfluid correlations in space and infinite 
superfluid susceptibility; and (iii) a superfluid phase with the usual off-diagonal 
long ranged order. It was predicted that the gapless Bose-glass intervenes between 
the Mott and the superfluid phases, so that the superfluid transition always occurs 
from the Bose-glass. The critical properties of this transition were characterized 
by three exponents: a dynamical exponent that equals the dimensionality of the 
system, a correlation length exponent that is bounded from above and an order 
parameter exponent that is bounded from below. 

Previously, the onset of superfluidity in a random potential was studied in |338| 
and [339j . In |338j . a system of strongly disordered hard-core bosons was consid- 
ered in the framework of a mean field theory including quantum fluctuations. A 
renormalization group approach was developed to study a one-dimensional system 
of interacting bosons in a random potential [339j . In this work it was shown the ex- 
istence of a localized-superfluid transition and universal power laws for correlation 
functions on the transition line were found. The case of zero temperature interact- 
ing bosons at commensurate density (one atom per site) on disordered lattices in 
one and two dimensions was addressed in |340j by using a real-space renormaliza- 
tion group. The results showed that when weak disorder is introduced, a transition 
directly from the Mott insulator to the superfluid occurs, so that infinitesimally 
weak disorder does not stabilize a Bose-glass at commensurate filling. The Bose- 
glass is found beyond a threshold disorder, in contradiction with the arguments 
in (SO]- The critical exponents for the supcrfluid-Bose-glass phase transitions at 
zero temperature for hard-core bosons in one-, two-dimensional [341j and three- 
dimensional [342| disordered lattices were calculated using a quantum real-space 
renormalization group method once the system was mapped onto a quantum spin- 
1/2 XY model with transverse random field. From these calculations |341[|342| . 
it was concluded that randomness is always relevant in one-dimensional systems 
while in two and three dimensions, there is a critical amount of disorder below 
which the superfluid phase is stable. Moreover it is also stated that there is only 
one universality class for the superfiuid-Bose-glass transition. Starting from the 
mapping of a system of hard core disordered bosons onto a quantum spin-1/2 XY 
model with transverse random field and generalizing it to a system of spins with 
arbitrary magnitude, a perturbative study was applied to get more insight into the 
low-energy excitations in a weakly disordered bosonic system [343] . 

In the early 90's, disordered Bose condensates were studied using the Bogoliubov 
approximation [344ll345ll346| . In [344] . the screening of the random potential in a 
two-dimensional dense Bose gas in a lattice due to the short-range repulsive inter- 
actions was addressed. This screening is effective when the healing length is short 
enough so that the condensate can adjust to the variations of the random potential. 
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The effects of the fluctuations deplete the condensate nonuniformly and lead to a 
spectrum of collective, phonon-like excitations. In |346j was reported that weak 
disorder hardly affects the condensate fraction or the superfluid density although 
the condensate distorts to screen the imposed random potential. In the strong dis- 
order limit, the condensate fraction and the superfluid density tend to vanish as 
the disorder increases and a constant density of states appears at low energy. In 
case of three-dimensional bosonic systems in random external potentials the de- 
pletion of the condensate and the superfluid density were explored in [345| . The 
conclusion obtained from this analysis was that disorder is more active in reducing 
superfluidity than in depleting the condensate. In this scenario, a formalism based 
on dispersive quantum hydrodynamics at zero temperature was applied to investi- 
gate the propagation of phonons in the system [347] . In particular, the shift of the 
sound velocity and its damping were calculated. In |348j . building on the contin- 
uum model of |345j and [347] , a systematic diagrammatic perturbation theory for a 
dilute Bose gas with weak disorder at finite temperature below the superfluid tran- 
sition temperature was developed and the disorder-induced shift of the superfluid 
transition temperature was derived. 

Monte Carlo techniques (world-line algorithm) have been also applied to study 
the interplay between interactions and disorder in bosonic systems. In |262] was 
reported the first convincing evidence of a second insulating phase, the Anderson- 
glass phase, in a one-dimensional lattice system. The Anderson-glass phase appears 
at weak couplings where the interactions compete with the disorder and tends to de- 
localize the bosons, contrary of what happens in the Bose-glass phase, that appears 
at strong couplings where disorder and interactions cooperate. Moreover, these two 
phases differ substantially on the nature of the boson density distribution: in the 
Bose-glass, the density is reasonably uniform while in the Anderson-glass, boson 
density correlations are expected to decay exponentially |262| . The existence of two 
such separate insulating phases was conjectured previously in [339j . The same year 
1991, in |349| . path- integral Monte Carlo techniques were used to study also the 
superfluid-insulator transition but in a two-dimensional square lattice. Using this 
technique in two dimensions, three phases were predicted: superfluid, Bose-glass 
and Mott-insulator and it was stated that at commensurate density, the system 
seems to undergo a direct transition from the superfluid phase to Mott insulator 
contradicting the picture of disorder of [30| and [339| . but in agreement with the 
results obtained using rcal-spacc renormalization group by [340] . The direct tran- 
sition from superfluid to Mott insulator without intervening Bose-glass at weak 
disorder was also reported for instance in [186] by using density-matrix renormal- 
ization group in one dimension, in [350j using quantum Monte Carlo simulations 
in two dimensions, or in (351^ using a mean field theory. Contrarily, the transition 
via the Bose-glass phase was predicted, for instance in |3521I353| by using Monte 
Carlo simulations, in |354| by using the renormalization group method, in [193j 
using a strong-coupling expansion for the phase boundary of the Mott-insulator or 
in [355j by an improved application of the density matrix renormalization group 
with respect to |186| . It is worth noticing that in the results of [355] there is no in- 
dication of a qualitative difference between the glass phase at small or large values 
of the repulsion i.e., between Anderson and Bose-glass. Thus, a complete under- 
standing of the phase diagram of interacting bosons in the presence of disorder 
and in various dimensions is still under debate. In particular, the possibility of a 
direct Mott-insulator to superfluid transition in the presence of disorder remains a 



47 



controversial issue. 



4.3 Realization of disorder in ultracold atomic gases 

In order to perform a detailed analysis of the properties of the disordered interacting 
Bose lattice gases it would obviously be very useful to have an experimentally 
accessible system, that could be studied in a controlled way. As we will see below, 
ultracold atomic lattice gases provide such opportunity. 

As has been discussed in section 1, since the experimental realization of Bose- 
Einstein condensation there has been an enormous progress in the studies of the 
ultracold gases: first in the weakly interacting regime and more recently also in 
the strongly interacting regime. Nowadays there exist a complete control of the 
generation and manipulation of ultracold bosonic, or fermionic g well as 

their mixtures. Among all the techniques of control, the transfer of these ultracold 
samples to optical lattices offers an unprecedented possibility to study disorder re- 
lated phenomena. In fact, by superimposing laser beams from different directions 
and with different frequencies, it is possible to generate a huge variety of lattice ge- 
ometries in a very controlled way. For instance, it has been proposed the use of two 
colour superlattices [551ISM7D] . i.e., the superposition of two standing- wave lattices 
with comparable amplitudes and with different wavelengths, as a form of quasidis- 
order. The so-called quasicrystal optical lattices in two and three dimensions have 
also been explored in [356II357II358II359| . These systems present long range order but 
not translational invariance. An example of a laser configuration [359] that gives 
rise to a quasicrystal lattice consists on Ni, laser beams arranged on the xy plane 
with iVfj-fold symmetry rotation. The polarization ej of laser j with wavevector kj 
is linear and makes an angle aj with the xy plane. The optical potential is in this 
case: 



where r = [x,y], < Sj < 1 stand for eventually different laser intensities and 
ipj are the corresponding phases. In [SS] a two-dimensional local quasidisordered 
potential was numerically simulated by using a main optical lattice, Vi (f) , to which 
a secondary lattice, K-(r), with much smaller amplitude and with frequency incom- 
mensurable with respect to the main lattice, is superimposed. The quasidisorder 
potential in this case reads: 



V{f) ^ Vi{t^ + Vr{r) = Vo[coa^{kx) + cos^{ky)] + Vi[coa^{kir) + cos^{k'2f)] (61) 



where f = [x,y]. The secondary lattice is responsible for the introduction of the 
(quenched) pseudodisorder which is determined by the ratio between the wave- 
lengths of the main and additional lattices, ki/k and k2/k. The same kind of 
potential has been also studied in one-dimensional geometries |360| . 

All the potentials discussed so far in this section are not strictly speaking disordered, 
but quasidisordered. Truly random potentials can be achieved by using a speckle 



pattern ^^^^'^W^JT^J^. The speckle field is a light field with highly 
disordered intensity and phase distributions but stationary and coherent. Such a 
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(60) 



48 




Figure 7. Theoretical prediction of Bose-Anderson glass in a ID optical lattice; top: superfluid 
fraction as a function of disorder Vi (thickline), respectively quasi-disorder V2 (thinline); middle: 
disordered speckle potential (thickline), and quasi-disordered (thinline) potential formed by 
incommensurable supcrlattices; bottom: occupation numbers of the ground state in the presence of 
disorder (black dots), and quasi-disordcr (open dots) (from \tiH\ ). 

speckle field can be easily generated experimentally, for instance, by introducing a 
diffusor in the path of a laser plane wave or by reflecting a laser beam in a surface 
that is rough on the scale of the laser wavelength. Specifically, the random intensity 
of an speckle field follows an exponential statistical distribution with a standard 
deviation given by the average intensity, P{I) = exp(— //(/)). In addition, the 
intensity correlation length or disorder correlation length, given by the half-width 
of the autocorrelation function, is limited by optical resolution being of the order of 
at least few ^m (see Fig.[7]for examples of ID random, and quasi-random potentials 
leading to Bose-Anderson glass). 

One can also create a disordered potential by using atoms of a second species or 
internal state quenched randomly at the nodes of an optical lattice |3641l365j . These 
randomly distributed atoms act as point-like scatterers for the atoms of the first 
species, that are not trapped by the optical lattice. One possibility to randomly trap 
the atoms of the second species is by rapidly quenching them from the superfluid 
to the Mott insulating phase. In this approach, the correlation length is very short 
since the distance between the lattice sites is of the order of 0.5/im. It is possible that 
such quenched random scatterers of fermionic composites (consisting of fermions 
bounded to 6-7 bosons) have been recently realised in an experiment with Fermi- 
Bose mixtures by the Hamburg group [36| . It is also worth pointing out that with 
this kind of disorder, the effective Hamiltonian for composites will involve disorder 
with bimodal, or at most trimodal distribution. This kind of disorder is particularly 
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interesting to study spin glass transition in the 2D Ising model [366| . 
Recently, also the generation of disorder in ultracold atoms in optical lattices 
through on-site interparticle interactions has been addressed j367il368j . The idea 
is to exploit the significant modifications of the scattering properties of two atoms 
that can be induced near Fcshbach resonances [MlIM] by slight modifications of the 
magnetic field. In this scenario, spatial random variations of the local interatomic 
interactions arise. The system is placed at the verge of a Feshbach resonance by 
means of an offset magnetic field in the presence of a spatially random magnetic 
field. This spatially random magnetic field appears for instance in magnetic micro- 
traps and atom chips as a result of the roughness of the underlying surface [369| . 
The intrinsic disorder that appears in magnetic microtraps and atom chips has also 
been addressed in the absence of an optical lattice |370[I371) . Controllable disorder 
on the level of next-neighbors interactions can be generated by means of tunneling 
induced interactions in systems with local disorder [86| . It is also worth mentioning 
a very recent attempt to create controlled disorder using optical tweezers meth- 
ods [37g. 



4.4 Disordered ultracold atomic Bose gases in optical lattices 

Previous section 14.31 reports on experimental feasible realizations of random and 
pseudorandom potentials for cold atoms in optical lattices. All these possibilities 
lead to experimental realizations of the disordered Bose-Hubbard model [551I69II70] , 
where the interplay between interactions and disorder could be explored. The Bose- 
Hubbard hamiltonian in the presence of disorder reads like Eq. p2[) but with site 
dependent tunneling rate,tij, and local chemical potential,^;: 

^ = - E l^^Mbj + h.c] + J2 ^Mn^^ - 1) + E A*^"'' (62) 

<ij> i i 

Although tunneling coefficients in the most general disordered case should be site 
dependent, it has been shown [68j that in optical lattices the hopping disorder is 
suppressed against on-site disorder included in the term /i,;. The last term in (|62p 
gives the on-site single-particle energy which originates from the external poten- 
tials and the on-site part of the kinetic energy. Therefore it accounts for an external 
harmonic trapping potential plus the inhomogeneities produced by the speckle pat- 
tern, or by the superlattice. Eventually, also contains the chemical potential in 
the grand canonical description. In [68j . the pseudorandom on-site energies are 
calculated by using: 

= J d^rw*{r~ fl)Vr{f)w{r- fl), (63) 

where w(f—T\) are the Wannier functions in the lowest Bloch band and Vr{r) is the 
superlattice potential introduced in or the potential induced by a numerically 
generated speckle pattern characterized by its mean value and the average speckle 
size. In [BniED] there was considered a sinusoidal variation of the in the range 
[—A, 0] recovering the regular case for A = 0. In [B^, the dynamical generation of 
the Bose-glass (strong interactions and non-integer filling factor) and the Anderson 
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Figure 8. Numerically simulated dynamical transition to the Bose glass state; first (in inset) a 
superfluid at high value of lattice potential is formed, and then (in the main figure) the disorder is 
turned on gradually. Condensate (solid line) and superfiuid (dashed line) fractions tend to zero 

(from leSl ). 



glass phases (weak interaction regime) in a two-dimensional ultracold bosonic gas 
was calculated using the Gutzwiller ansatz method (see Fig. [5]), while in ^69ii70] the 
ground state of the system was determined by solving the eigenvalue problem of 
the one-dimensional Bose-Hubbard hamiltonian numerically. In this last case, also 
both the Anderson-glass and the Bose-glass phases were found (see Fig.[S|). 

Very recently [159] . the problem of disordered one-dimensional bosonic systems 
with hard core boson interactions has been exactly solved via a Jordan- Wigner 
transformation (see subsection 13. 8p . This mapping establishes a connection to non- 
interacting disordered fcrmions and allows to take advantage of many known results 
on Anderson localization. In this scenario, the correlation functions of the particle 
density and the local density of states for the interacting bosonic case coincide with 
the non-interacting fermionic ones, because only depend on the modulus square of 
the wave-function. This is not the case for the momentum distribution, but the 
mapping allows for a simple calculation of the disorder-averaged boson momentum 
distribution [159| . These calculations show the complete destruction of quasi-long 
range order by disorder and the flattening of the momentum distribution for suffi- 
ciently strong disorder. Note that in the homogeneous case, the momentum distri- 
bution is well known to possess a singularity [3731137411375] . In |159j it was also shown 
that the Bose-Fermi mapping can also be established via the effective low-energy 
theory [339] . Moreover, it was pointed out that a similar mapping is also available 
for arbitrary interaction strength, but involves interacting fermions with nonstan- 
dard contact interactions [376] . For strong (but finite) repulsive bosonic interac- 
tions, the weak fermionic interactions can be treated perturbatively. In Ref. |378] . 
specific effects of quasi-disorder, in contrast to real disorder, in ID lattices were 
discussed. 
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Figure 9. Contour plot of the square of the largest coefficient, C^^xi of the expansion of the state 
of the system described by a Bose-Hubbard model in the number basis, as a function of disorder A 
and on-site interaction V , in units of the hopping energy J. The results correspond to exact 
diagonalization of a one-dimensional system with 8 sites and 8 bosons (from |69) . Note that V 
accounts for U and J for t in the figure). The labels identify the predicted different phases: 
superfluid, Mott insulator, Bose glass (BG), and Anderson localization (AL). 



Just before submission of this review, Krutitsky ct al. [377] calculated the mean 
field phase diagram of disordered bosons at finite T . They made the important 
observation that the distinction of MI and Bose glass at finite T in atomic lattice 
gases should be possible by looking at the density of states of low energy excitations. 

The experiments on disordered ultracold gases in optical lattices has just started, 
and most of them deal with trapped gases and weak interaction limit (we discuss 
them in the next subsection). Nevertheless, spectacular results have been achieved 
in the area of disordered lattice gases, already. The Florence group has recently 
applied to a basic lattice formed by light with wavelength Ai = 830nm another 
superlattice created by A2 ~ 1076nm light, creating a quasi-disordered potential 
(see Fig. fTU)) . They have applied the method of Bragg spectroscopy developed in [5S] 
and measured effects of the lattice modulation. This allows essentially to measure 
the low energy excitation spectrum; broadening of this spectrum was identified in 
Ref. [3S] to be the signature of the Bose glass phase (see Fig. [Tl|) . 

In the experiment of the Hamburg group |36| (see also |37| ) , an unexpected shift of 
the SF-MI transition point towards weaker lattice potentials has been observed in 
a Fermi-Bose mixture with attractive interspecies interactions. It has been specu- 
lated that this effect might have to do with the formation of composite fermions. 
Composite fermions, consisting of a fermion and few bosons, are hardly mobile and 
should remain in the conditions of the mentioned experiments immobile and play 
a role of random scatterers for (still relatively movable) bosons. 
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Figure 11. Experimental signatures of the MI— Bose glass transition: Broadening of the excitation 
spectrum as disorder (52) grows. One observes discrete linbes in the gapped MI, and significantly 
broadened and structureless spectrum in the BG phase (from 1350 . 



4.5 Experiments with weakly interacting trapped gases and Anderson 
localization 

During the last year the first experimental results on Bose-Einstein condensates in 
random potentials were reported [7Tl[7^l751F7i] . In the first experiments [7T] , static 
and dynamic properties of an harmonically trapped ^^Rb condensate in the pres- 
ence of a stationary, but spatially random potential created by imaging an speckle 
pattern onto the condensate were addressed. In this experiment the smallest length 
scale of the speckle was L^is = 10/im while the size of the quasi one dimensional 
condensate was 110/xm in the axial direction and ll/um in the radial. Absorption 
images after free expansion for different values of the speckle height, Vg, reveal dif- 
ferent regimes of behavior: (i) for strengths of the disorder much smaller than the 
chemical potential of the condensate in the harmonic trap, practically no change 
with respect to the standard Thomas- Fermi profile is observed; (ii) for intermediate 
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values of Vs, stripes appear in the expanded density profile, which can be a signa- 
ture of the development of different phase domains across the condensate, like the 
ones associated with phase fluctuations in highly elongated condensates |284'. !379j , 
or they can be explained through the presence of different momentum components 
due to the growth of instabilities [380j ; and (iii) when Vg is larger than the chem- 
ical potential of the harmonically trapped condensate, the tunneling between the 
different minima of the random potential is strongly suppressed and the system 
enters the tight binding regime, where a broad Gaussian profile is seen in the ex- 
pansion as expected for the expansion of randomly spaced condensates isolated 
in the individual speckle wells. The transport properties of the system were also 
studied in |71| through collective excitations, dipolar and quadrupolar modes. The 
dipolar mode was excited by abruptly displacing the magnetic trap in the axial 
direction obtaining: (i) for Vs much smaller than the chemical potential, slightly 
damped oscillations were obtained with the same frequency as in the undamped 
case in the absence of the speckle; (ii) for intermediate Vg, strong damping of the 
oscillations is observed; and (iii) for high values of Vg, the atomic cloud docs not 
oscillate and remains localized on the side of the magnetic trap. The quadrupolar 
mode was excited with a resonant modulation of the radial trapping frequency ob- 
serving that the anharmonicities of the random potential result in frequency shifts 
that are not correlated with the dipolar frequency like it happens in the absence 
of speckle. Moreover, it is shown that both the sign and the amplitude of the shift 
depend on the exact realization of the speckle potential. 

Following these first experiments, the suppression of the one-dimensional transport 
of an interacting elongated condensate in a random potential was reported nearly 
simultaneously by [721174] . In these experiments, one-dimensional expansion of a 
Rb condensate along an optical [71] or magnetic [73] guide in the presence of a 
speckle potential was studied. Without the speckles, the condensate freely expand 
and the growth of the root mean squared (rms) radius (rms= < > + < y"^ >) 
is self-similar and linear in agreement with the theory [38111382] . In the presence of 
a random potential with high enough amplitude, but without entering in the tight 
binding regime, the expansion dynamics changes completely and both the expansion 
and the centre of mass motion arc inhibited. Although this strong suppression 
of expansion corresponds to disorder-induced trapping of the BEG, it does not 
correspond to Anderson like localization. The reasons are mainly: (i) the screening 
of the disorder potential played by the interaction energy; (ii) the resolution of the 
optics of the disordered potential that fixes the correlation length of the disorder, 
Ldis, which is much larger than the healing length, Iheai = 1/ VSirna where n is the 
density and a the atomic scattering length. The healing length provides the typical 
distance over which the order parameter of the condensate recovers its bulk value 
when it is forced to vanish at a given point by an impurity for instance. Therefore, 
the density profile tends to follow the modulations of the disordered potential; 
and (iii) the fixed order of magnitude of the distance between the speckle sites 
imposes that the typical axial size of the system, L, is only ten times larger than 
the correlation length of the disorder. In order to achieve Anderson localization, the 
correlation length of the disorder has to be smaller than the size of the system, which 
seems to be difficult to achieve with speckles. Anderson localisation of elementary 
excitations in ID BEG seemingly sets a little less rigid requirement, but is still not 
easy to achieve with the presently achievable speckle patterns [771[7H]. Recently, a 
very detailed study shows that by expanding a quasi 2D cloud to a size of lew?, it 
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Figure 12. Density profiles (absorption images) 
harmonic trap plus random speckle (left column), 
lattice, and random speckle (right column). The 
third row shows the result of the n 



of quasi ID BEC released from the combined 
and from the combined harmonic trap, optical 
second row shows the column density and the 
umerical simulation (from [73]). 



should be possible to see effects of weak localization and perhaps even strong ones, 
using available speckle potentials [77]. 

In Ref. the first realization of an ultracold disordered lattice gas was reported, 
so the effect of a speckle pattern was superimposed to a quasi- ID condensate in 
a regular ID optical lattice. In this work, it was shown that the fragmentation, 
already reported in [7T] , also appears in the presence of an optical lattice (see Fig. 
[T2|) . More interestingly, in [TSj it has been explored the crossover from the Anderson 
localization in the absence of interactions, where the ground state wave function 
is characterized by an exponential localization, to the screening regime, where the 
number of localization centres is so high that one can no longer distinguish the 
individual localized states and the signature of nontrivial localization vanishes (see 
Fig. I13p. On one hand, the interactions can be reduced by reducing the number of 
atoms, lowering the trap frequencies or tuning the scattering length via Feshbach 
resonances and on the other hand, the limitations in the ratio between the corre- 
lation length of the disorder and the healing length and between the correlation 
length of the disorder and the size of the system, that arise when the disorder 
potential is created using an speckle pattern, could be overcome using the quasidis- 
order created by several lasers with incommensurable frequencies. Measurement of 
the superfluid fraction in an accelerated optical lattice is proposed [73| as a way 
of detection of the localization rather than the usual measurement of the density 
distribution after ballistic expansion. 
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Figure 13. Calculated condensate wavefunction for a BEC in a combined potential formed by 
harmonic trap, and two incommensurate standing waves. The sequence (a)-(c) shows the effects of 
the increasing nonlincarity (i.e., increasing the number of atoms) (from 1731 ). 

4.6 Disordered interacting fermionic systems 

So far wc have discussed disordered interacting bosonic systems and the experimen- 
tal possibihties of testing the still controversial points of the theory using ultracold 
bosonic systems in optical lattices. Fermionic gases have also been cooled to the 
quantum degeneracy regime using sympathetic cooling of two fermion species, or 
boson- fermion mixtures (for the first experiments see |383 !!384"385"386l! 387j ) . More- 
over, fermionic ^"^K atoms have been also loaded into a three-dimensional optical 
lattice [Mj . Several proposals of using ultracold fermionic atoms in optical lattices 
to study various condensed matter models have already been formulated. In [50] 
for instance, it was discussed how fermionic atoms in optical lattices allow for the 
realization of antiferromagnetism or high-temperature superconductivity. More re- 
cently, it has been shown ^388] that the fermionic SU(N) Hubbard model on a 2D 
square lattice can be realized with ultracold fermions in an optical lattice. Regard- 
ing the disordered case, very recently it has been reported [389| that fermionic 
atoms in an optical superlattice exhibit strongly correlated phases from Kondo 
singlet formation to magnetism of localized spins. 

Let us now briefly review the literature on the physics of disordered interacting 
electronic systems with the aim of identifying the still open questions that could be 
addressed exploiting the ultracold fermionic atoms in disordered optical lattices. 

As we have discussed in subsection 14. 11 originally disorder in non- interacting elec- 
tronic systems attracted a lot of interest, but the effects of electron-electron in- 
teractions received full attention only recently. It was believed that weak disorder 
should not modify essentially the Fermi liquid picture of Landau. Later a Fermi 
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liquid theory for electrons in disordered solids was formulated, and termed the 
Fermi glass theory [3901 1391] . According to this approach the Landau's quasiparti- 
cle description is valid but now the quasiparticle wave-functions are not extended, 
as in the translationally invariant case, but they have no long-range coherence and 
may be even localized. Each quasiparticle is viewed as a single entity moving in 
the self-consistent field of all the other quasiparticlcs, and the resulting theory has 
the form of a mean-field theory for the quasiparticlcs. In [391] , the phase transition 
into a Fermi glass state was discussed. In |392|j393ll394] , and independently in |395| , 
it has been shown, applying perturbation theory to the lowest order in the inter- 
action strength, that even weak disorder leads to surprisingly singular corrections 
to the electronic density of states near the Fermi surface, and to transport prop- 
erties. In 1994, Shcpelyansky [396] stimulated further discussion about the role of 
interactions by considering two interacting particles in a random potential, arguing 
that there exist an interaction-induced enhancement of the two-particle localization 
length compared to the noninteracting case. The lines of this work were developed 
further [5^71139811395111001111111 1111^ . A Fermi liquid approach for finite densities was 
suggested by [397] . and later developed by [400tl402j reducing the problem to the 
study of the delocalization of few quasiparticlcs above the Fermi sea. 

Several groups have tried to study effects of interplay between disorder and (repul- 
sive) interactions in more detail in the regime, when Fermi liquid becomes unstable 
as the Mott insulator state is approached by increasing the interactions. In other 
words, the crossover between Fermi glass (Anderson localized states populated 
with the restriction of the Pauli principle) and Mott insulator (state where repul- 
sion dominates over kinetic energy and disorder) has been studied. Density matrix 
renormalization group studies were performed for spinless fermions with nearest 
neighbor (n.n.) interactions in a disordered mesoscopic ring |403| . and for spin 1/2 
electrons in a ring described by the half-filled Hubbard- Anderson model [404] . Spin- 
less fermions with Coulomb repulsion (reduced to n.n. repulsion) in 2D [40511406] 
were also studied. This collection of works shows that as interactions become com- 
parable with disorder, delocalization does takes place. In a ID ring it leads to the 
appearance of persistent currents. In 2D, the delocalized state exhibits also an or- 
dered flow of persistent currents, which is believed to constitute a novel quantum 
phase corresponding to the metallic phase observed in experiments for instance 
with a gas of holes in GaAs hctcrostructurcs for the similar range of parameters. 

Another intensive subject of investigation concerns metal (Fermi liquid)-insulator 
transition driven by disorder in 3D. Theoretical description of this phenomenon goes 
back to the seminal works of Efros and Shklovskii [4070 408,409 and Mac Millan 
[410] . In this context, recent results of experiments on disordered alloys |411[|412] 
allowed to determine the critical exponents that govern the conductivity dynamics 
on both sides of a quantum phase transition in a disordered electronic system. 
Weakly doped semiconductors provide a good model of a disordered solid, and 
their critical behavior at the metal-insulator transition has been intensively studied 
(cf. [41311414] ). Results concerning various forms of electronic glass: from Fermi 
glass (with negligible effects of Coulomb repulsion) to Coulomb glass [415II416II417] , 
dominated by the electronic correlations were obtained in the group of M. Dressel 

mm- 

Very recently, the ground state phase diagram at half filling, for arbitrary inter- 
action and disorder strength, has been calculated [420] by applying the Dynami- 
cal Mean-Field Theory (DMFT) with the geometrically averaged local density of 
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states. It was shown that the presence of disorder increases the critical interac- 
tion at which the Mott-metal transition occurs and turns the sharp transition into 
an smooth crossover. Regarding the critical disorder strength for the Anderson- 
localization transition, it was reported that it increases for weak interactions, and 
it is suppressed by strong interactions [420] . In the phase diagram obtained with 
the same method, but with an arithmetic average of the local density of states the 
Anderson transition is missing |421j . 

Again, we would like to point out that many of the questions discussed in this 
section can be addressed using ultracold atomic lattice gases. 

4.7 Disordered Bose-Fermi mixtures 

As mentioned in the subsection l4.61 in order to cool fermionic gases to the quantum 
degeneracy regime it is necessary to use sympathetic cooling of two spin species or 
boson-fermion (B-F) mixtures [MUSHUMillMllMTj . In the latter case, the final 
phase of the system is a quantum degenerate B-F mixture. The loading of a B-F 
mixture in an optical lattice has been also recently reported |422[[36l[37] . 
In the absence of disorder, and in the limit of strong atom- atom interactions such 
lattice B-F systems can be described in terms of composite fermions consisting 
of a bare fermion, or a fermion paired with 1 boson (bosonic hole), or 2 bosons 
(bosonic holes), etc. |167j . The physics of Fermi-Bose mixtures in this regime has 
been studied recently in a series of papers [MllIBllSinilllHlilUlll, where it 
has been shown that the low temperature dynamics of the fermionic composites 
is described by an effective Hamiltonian (see also Sec. 13.111 and Appendix [A]) . de- 
scribing a spinless interacting Fermi gas. The validity of the effective Hamiltonian 
for fermionic composites in ID was studied using exact diagonalization and the 
DMRG method [219j . The effects of inhomogcneous trapping potential on lattice 
mixtures was for the first time discussed by Cramer et al. |220] . while the disor- 
dered case in the strong coupling limit was studied in [551112 • In the presence of 
disorder, degenerate second order perturbation theory cannot be applied to derive 
the effective Hamiltonian, as it was used in |166lll67j . since even for zero hopping 
rates there exists a well defined single ground state determined by the values of 
the local chemical potentials. Nevertheless, in general, there will be a manifold of 
many states with similar energies. The differences of energy inside a manifold are 
of the order of the difference of chemical potential in different sites, whose random 
distribution is bounded. Moreover, the lower energy manifold is separated from the 
exited states by a gap given by the boson-boson interaction. Therefore, one can ap- 
ply a form of quasidegenerate perturbation theory by projecting onto the manifold 
of near-ground states, as described in Appendix [Al |230U87| . As in the homogeneous 
case [167] ■ composite fermions behave as a spinless interacting Fermi gas, but in 
the presence of local disorder they interact via random couplings and feel effective 
random local potential. The effective Hamiltonian that describes their physics can 
be written as follows: 



where Fj are the annihilation operators of the composite fermions (Eqs. ([55)1 and 



Heff 




(64) 
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(|57p ). and Mi — F^Fi. The nearest neighbor hopping for the composites is described 
by — dij, the nearest neighbor composite-composite interactions is given by Kij, 
which may be repulsive (> 0) or attractive (< 0), and are the on-site energies. 
The exphcit calculation of the coefficients , Kij and Jl^ depends on the concrete 
type of composites fermions |86[[57] . For fermion-bosonic hole composite and for 
fermion-boson composite, it has been shown that the hopping amplitudes dij are 
always positive. Depending on the ratio a = V/U between the boson- fermion and 
the boson-boson interactions, the effective interactions between composites /Vy may 
be either repulsive or attractive for all the values of disorder, or for certain values of 
a the qualitative character of the interactions is controlled by the inhomogeneity. 
These two types of composites have been studied in two limits: (i) the small disorder 
limit, where the contributions of the disorder to the interactions and hopping are 
neglected {Kij ~ K and dij = d), and only the leading contributions in the on- 
site energies arc kept; and (ii) the large disorder limit or spin glass limit, where 
the tunneling becomes non-resonant and can be neglected in Eq. (j64|) . while the 
couplings Kij fluctuate strongly. This situation corresponds to the (fermionic) Ising 
spin glass model (see Section for details). 

For the case of disorder applied only for the bosonic component of the mixture, and 
in the small disorder limit (when disorder does not affect the composite formation), 
the following quantum phases can be achieved: 

• Fermion + bosonic hole composites 

When K/d ^ 1, i.e., when the interactions are negligible, the system is in the 
Fermi glass phase, i.e. Anderson localized (and many-body corrected) single par- 
ticle states are occupied according to the Fermi-Dirac rules. For large repulsive 
interactions, K/d ^ 1 and K > 0, the ground state will be a Mott insulator and 
the composite fermions will be pinned for large filling factors. The value of K/d 
is, however, bounded from above by 2, which seems not enough to achieve the 
Mott insulator state. For intermediate values of K/d, with K > 0, delocalized 
metallic phases with enhanced persistent currents are possible. For attractive 
interactions {K < 0) and l-ftTj/d < 1 one expects competition between pairing of 
fermions and disorder, i.e., a "dirty" superfluid phase while for \K\/d^ 1, the 
fermions will form a domain insulator, that is a state in which fermionic compos- 
ites will stick together to form a rigid immobile cluster. In Ref. [5^ the crossover 
from the Fermi gas to the Fermi glass phases has been studied numerically by 
means of the dynamical Gutzwiller ansatz method (Fig. [14]); similarly, dynamics 
of a transition from the fermionic domain insulator to a disordered insulating 
phase was investigated (Fig. [T5|) . 

• Fermion -|- boson composites 

The regimes where K <^ d lead to a non-interacting Fermi glass while the regimes 
of strong effective repulsive interactions, where K d and K > 0, would cor- 
respond to Mott insulator or checkerboard phase if the filling factor is 1/2. In 
this case, no strong attractive interactions regime occurs since K/d reaches a 
minimum of ~ —0.07. Therefore, the domain insulator phase does not appear, 
and even the "dirty" superfluid phase may be washed out. 

The summary of the possible phases, described in the case of fermions plus bosonic 
hole composites, and in the case of fermion plus boson composites, is presented 
in the schematic diagram of Fig |16l In the case of bare fermion composites, it 
has been shown that for finite boson-fermion interactions, the fluctuations of the 
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Figure 14. Numerically calculated dynamical transition from the Fermi liquid to the Fermi glass of 
fermionic composites (fermion + bosonic hole composites) in a 2D lattice with Af = 10 X 10 sites. 
In (a) , the decrease of the variance of the number of fermions per lattice site is shown as a function 
of the amplitude of the disorder. In (b), the probability of having one composite at each lattice site 
in the absence of disorder is given, and (c) gives the same as (b) after adiabatically ramping up the 

disorder (from \87\ ). 



effective composite interactions may be large, and the dynamics of this type of 
composites resembles quantum bond percolation. One can assume, in a somehow 
simplified view, that the interaction parameter Kij takes either very large, or zero 
values. The lattice decomposes then into two sub-lattices: a "weak" bond sub-lattice 
(corresponding to Kij <^ dij) in which fermions flow as in an almost ideal Fermi 
liquid, and a "strong" sub-lattice (corresponding to Kij ^ dy ), where only one 
fermion per bond is allowed. Additionally, for the case of lattices with different 
types of sites (i.e. sites in which disorder affects the formation and character of 
composites), it has been predicted that physics of quantum site percolation will 
become relevant. 
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Figure 15. Numerically calculated dynamical transition from the Fermi domain to the Fermi glass 
of fermionic composites (one fermion and zero boson) in a 2D lattice with A'^ = 10 X 10 sites. In (a), 
the decrease of the variance of the number of fermions per lattice site is plotted as a function of the 
amplitude of the disorder. In (b), the probability of having one composite at each lattice site in the 
absence of disorder is given, and (c) gives the same after adiabatically ramping up the disorder 

(from [87]'). 



4.8 Spin glasses 

Spin glass is a phase that appears in spin systems interacting via random couphngs 
that can be positive (anti- ferromagnetic), or negative (ferromagnetic). Such varia- 
tions of the couphngs lead typically to frustration i.e., if there are only two possible 
spin orientations and the interactions are random, no spin configuration can simul- 
taneously satisfy all the couplings (sec section [5] for more details). Thus, spin glass 
behavior requires two essential ingredients: quenched disorder and frustration. 

In a somewhat oversimplified picture, this kind of systems are characterized by two 
order parameters: (i) the order parameter for magnetic ordering, the magnetization 
given by M :=< o7 > and (ii) the Edwards- Anderson order parameter for spin glass 
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Figure 16. Schematic predicted phase diagram of (a) fermion plus bosonic hole composites and of 

(b) fermion plus boson composites 1871). 



ordering, Qea '■— < <Ji where (•) denotes the average over disorder and < (•) > 
the thermal average. Experimental studies have identified the phase diagram of the 
spin glasses consisting on three phases |425| : (i) At high temperature and small 
average spin exchange, the system exhibits a paramagnetic phase with M = 
and Qba = 0; (ii) for K < and large, the system exhibits a ferromagnetic phase 
with M 7^ and Qea 0', and (iii) for weak K and small temperatures, a spin 
glass phase appears with M = and Qea 7^ 0. Qea 7^ signals that the local 
magnetization is frozen, but may vary from site to site, so that the disorder prevents 
long range magnetic order. Spin glass phase is an example of "order in the presence 
of disorder" ^5] , 

The physics of spin glasses is up to now not well understood and remains as one 
of the challenges of statistical physics. In particular, the question of the nature 
of their order is still open. There exist two competing theories: the Mezard-Parisi 
(MP) model, and the droplet model. The MP picture is a mean-field theory based 
on the replica method [53] that predicts that the spin glass phase consists of a 
large number of low-energy states with very similar energies. The applicability of 
the MP picture for short range spin glasses (like the Edwards- Anderson model [81]) 
is very questionable. The rival theory, the droplet model [5^[551H^I427ll332j is a 
phcnomcnological theory based on scaling arguments and numerical results. The 
droplet model predicts that there are two ground states related by spin-flip symme- 
try, and that excitations over the ground state are regions with fractal boundary, 
the droplets, in which the spins are inversed with respect to the ground state. 

One of the most remarkable results of Refs. [551157] is the possibility of the realiza- 
tion of a fcrmionic Ising spin glass model by using a disordered Bose-Fermi mixture. 
The spin glass limit is obtained from the composite fermionic model (Ea. (|64p ) when 
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the interactions between fermions and bosons are of the same order, but shghtly 
smaller than the interactions between bosons, and in the limit of large disorder. 
In this situation, the hopping vanishes due to strong site-to-site energy fluctua- 
tions and the nearest neighbor interactions, Kij, fluctuate around mean zero with 
random positive and negative values. Replacing the composite number operators 
with a classical Ising spin variable cr^ := 2Mi — 1 = ±1, one ends up with the 
Hamiltonian: 

He- A = 4 X! ^io^^^i ^ 2 ^ '^^^^ 

(u> 

It describes an (fermionic) Ising spin glass [428|,l429j . which differs from the stan- 
dard Edwards- Anderson model |430i[55] in that the model fEq. (p5|) ) includes an 
additional random magnetic field Jl^, and, moreover, it has to satisfy the constraint 
of fixed magnetization value, m = 2Nf/N — 1, as the number of fermions, Nf, 
in the underlying BFH-model is conserved. N is the number of sites. It shares, 
however, the basic characteristics with the Edwards- Anderson model as being a 
spin Hamiltonian with random spin exchange terms Kij. In particular, this pro- 
vides bond frustration, which in this model is essential for the appearance of a 
spin glassy phase. Due to the mentioned differences, it is necessary to reformulate 
slightly the Mezard-Parisi mean field description of the system by adapting the 
Sherrington-Kirkpatrick-Parisi calculations [SH[S3j to this specific case [57] . 

The experimental study of this limit in ultra-cold Bosc-Fermi mixtures could 
present a way to address various open questions of spin glass physics concern- 
ing the nature and the ordering of its ground- and possibly mctastable states (the 
Mezard-Parisi picture [S^ versus the droplet picture [84l[85l l426, 427. 332J). broken 
symmetry and dynamics in classical (in absence of hopping) and quantum (with 
small, but nevertheless present hopping) spin glasses |1541l431j . Particularly inter- 
esting are prospects for studying questions that have been addressed only recently: 
the existence of the celebrated de Almeida-Thouless line (i.e. the line separating 
the replica symmetry broken spin glass phase from the replica symmetric phase 
in the magnetic field in the Parisi picture) in short range spin glasses |4321I433] . 
disorder chaos in spin glasses (i.e. ground state sensitivity to small changes of 
disorder, [IM]), spin glass transition in 2D and its dependence on the bond distri- 
bution [435j . universality classes |436j etc. (for a recent review see [366] V 

Usage of atoms offers unprecedented possibilities for detection of the spin glass 
properties. Recently for instance, it has been proposed [437j how to create replicas 
of disordered systems, i.e. systems with identical disorder landscape and to the 
measure of correlations between the replicas. These correlations can provide infor- 
mation about the ground state phases of the system. The replicas can be created 
in a three dimensional lattice with incommensurate frequencies by quenching the 
hopping between different planes obtaining then a set of two-dimensional incom- 
mensurate lattices, with the same realization of disorder in each plane. Another 
possibility is to use localized impurity atoms. In this case, the procedure would be 
first to prepare a disordered distribution of localized impurity atoms in such a way 
that only two or no atoms per lattice site are allowed. This can be achieved, for in- 
stance, using repulsive bound pairs, observed recently by the Innsbruck group |39] . 
Then, periodicity of the lattice in one direction should be adiabatically doubled by 
stretching the lattice in that direction. One would in this way obtain an array of 
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Figure 17. (a) Antifcrromagnetic Ising model (J > 0) on a triangular lattice: Orientation of a spin 
cannot be determined by satifying all the couplings simultaneously, (b) Square lattice with all 
ferromagnetic Ising interactions, except one, in which the interaction is antifcrromagnetic. The 
interactions are indicated on the sides of the lattice. The model is again frustrated. 

replicated pairs of 2D random landscapes. That would be, obviously, an ideal tool 
to study the questions concerning existence of finite T spin glass transition in Ising 
model with bimodal bond distribution |366j . 

Refs. [5S1I57] demonstrate that ultracold disordered Bose-Fermi mixtures in optical 
lattices may serve as a paradigm fermionic system to study a variety of disordered 
phases and phenomena: from Fermi glass to quantum spin glass and quantum 
percolation. No doubts, one can mimic condensed matter and even go beyond! 



5 Frustrated models in cold atom systems 
5.1 Introduction 

Frustration appears when all the constraints imposed by the Hamiltonian cannot be 
simultaneously fuUfilled and it is an inherent property of some strongly correlated 
systems. It introduces various interesting features, such as a rich phase diagram. 
Particularly fascinating effects appear when a system is both disordered and frus- 
trated (see, e.g., [8311154] ). In this section, however, we will study several frustrated 
models in regular lattices and discuss their possible realizations. Antifcrromagnetic 
(AF) models in a regular lattice are frustrated if the geometry of the lattice is 
sufficiently complicated (such as triangular or kagome lattice). 
A physical system described by the Ising Hamiltonian 



is AF, if J > 0, where {ij) denotes the summation over nearest neighbors. This 
model, when considered on a two-dimensional triangular lattice, is an example of a 
frustrated system |438] . It is clear from Fig. [T7](a) that it is not possible to minimize 
the energy on such a triangular lattice, satisfying all the bonds simultaneously. In 
this case, wc obtain a degeneracy in the ground state due to the frustration. It is 
not due to the invariance under global spin flips of the Hamiltonian, which is the 
case for the ferromagnetic Ising model (J < 0) on the same lattice. 
Frustration can be met even in a lattice, which has a relatively simple structure 
(e.g., a 2D square lattice), if some interactions are ferromagnetic while others are 
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AF. In the ferromagnetic Ising model, e.g., if one replaces an odd number of ferro- 
magnetic bonds by AF ones, one obtains a frustrated model from a nonfrustrated 
one. The square lattice in Fig. I17r b) is an example of such a model. In general, 
we say that the model is frustrated, when the orientations of its spins cannot be 
obtained by satisfying all the couplings simultaneously to the aim of minimizing the 
energy in the model [83| . Further discussion on this subject can be found in |107| . 
Detailed studies of quantum antiferromagnets can be found in the books of Auer- 
bach [SD] and Sachdev [154] . Frustrated models have been reviewed very recently 
by Misguich and Lhuillier |107| (see also |108[|109| '). They have been also discussed 
in the context of spin glasses (see f430j ). The concept of frustration has been also 
studied from the perspective of high temperature superconductors (e.g., see the 
reviews g39] and [440j ). 

The plan of this section is as follows. Wc will present a short review on quantum 
antiferromagnets in the succeeding subsection. The next two subsections (|5.3I and 
15. 4p will be devoted to the kagome lattices: we will review some recent results and 
discuss possibilities of experimental verifications of the theoretical predictions. In 
the last subsection, we will discuss the possible realization of quantum magnets 
with cold atoms, ions, and molecules. 



5.2 Quantum antiferromagnets 

In this section we discuss some basic properties of the quantum antiferromagnets. 
The interest in this subject is motivated mainly by two things: (a) frustrated quan- 
tum antiferromagnets are related to most challenging open questions of condensed 
matter physics; (b) there exist already several theoretical proposals for experimen- 
tal realization of quantum magnets in cold atom/ion systems. We expect that these 
experiments will verify existing theoretical predictions, and point out new direc- 
tions in theoretical investigations of the quantum magnetism. 
According to reviews |1071ll08lll09] , quantum antiferromagnets exhibit at low tem- 
peratures one of four generic behaviors (see below). We add to the list a fifth 
possibility pTOllTTT] . 

(i) Nccl order. This is a standard up-down-up-down type of ordering that occurs 
often in regular lattices such as square in 2D or cubic in 3D. Generalized 
(planar) Ncel order might, however, occur also in frustrated antiferromagnets, 
for instance in the triangular lattice [44 1| . Neel order breaks rotational and 
translational symmetry, exhibits long range order, and leads to gapless spin 
wave excitations. 

(ii) Valence bond solids (VBS). This order breaks the lattice translational sym- 
metry and consists of an ordered covering of the lattice by singlets, called 
also dimers. There is a long range order in dimer correlations, but the system 
is gapped and correlation functions decay exponentially. A famous example 
of a VBS is the Affleck-Kennedy-Lieb-Tasaki (AKLT) state, although it does 
not break the lattice translational symmetry [HD] . Excitations in 2D VBS are 
gapped (confined spinous) [109] . 

(iii) Spin liquid of the I type. This order has no apparent symmetry breaking, 
and decaying correlations. Particularly interesting are topological spin liquids 
[10911442] , for which degeneracy of the ground state depends on the topology 
of the underlying lattice. 
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(iv) Spin liquids of the II type. This order has no symmetry breaking, is gapless, 
and has large density of states of low energy excitations. Heisenberg antifer- 
romagnet in trimerized kagome lattice, discussed below is an example of such 
behavior, representing in a resonating valence bond state (RVB). The concept 
of resonating valence bond liquid was first introduced by Anderson |443| , in 
the context of superconductivity. When a valence bond solid starts melting 
due to quantum fluctuations introduced in the system, the situation gives 
rise to a new phase called resonating valence bond liquid. Numerical Monte 
Carlo simulations suggest that the liquid does not show any long-range cor- 
relations [444] , and is a products of singlets [109j . 

(v) Quantum spin-liquid crystal. This order combines Neel order with spin liquid 
type II behavior. Numerical studies reveal [110illll[|112j that it has a large 
number of low energy excitations as in the spin-liquid type II, and presumably 
no gap. 



5.2.1 The Heisenberg model. Let us first focus on spin systems that do not 
show frustration. Introduction of such models will help us to explain many results 
in the frustrated models. Let us consider the Heisenberg model 



where a is the vector of spin-1/2 Pauli operators. 

The ferromagnetic case. For the sake of completeness, we first look at the 2D 
ferromagnetic model (J < 0) in a square lattice. The exact ground state is given 
by the state [155j 



where \Q)i and represent respectively the up and down eigenstates of cr^ at the 
i-ih site. 

The antiferromagnetic case. In the case of 2D antiferromagnets (J > 0) on 
the square lattice, the exact ground state is hard to find though the system is not 
frustrated. Nonetheless, the classical ground states, obtained after replacement of 
the (T operators by vectors, can be explicitly written down. To this aim the lattice 
is divided into two sublattices A and B in such a way that all nearest neighbors 
of sublattice A belong to the sublattice B, and vice- versa. Such a splitted lattice 
is called a bipartite lattice. The classical ground state can be obtained by taking 
some orientation (e.g., in the z-direction) of all spins in one sublattice and opposite 
orientation of all spins in the other sublattice. In this way, one can find the classical 
ground state of a square lattice [154[|445j : 




(67) 



(68) 



(69) 
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This classical ground state is known as the Neel state. Note here that any pair of 
spins in this classical ground state is either parallel or antiparallel, which means 
that the ordering is collinear. The ordering of classical ground states in frustrated 
systems, however, is noncoUinear, even in bipartite lattices |154j . 

In the case of one-dimensional AF model, the exact ground state of the spin-1/2 
antiferromagnet in the thermodynamic limit is known (Bcthe solution) |2941I447[ 
I4461I5D] . For comparing the properties of integer and half-integer spin systems [448) 
(cf. [449] ) . and the famous Haldane conjecture, the systems may be realized with 
atoms either by using Fermi-Fermi mixtures, or by employing the approach of 
Ref. EH. 

For quantum antiferromagnets, made up of spin-1/2 particles, and placed in the 
infinite square lattice, the ground state energy can be approximated by that of the 
Neel state (Ea. (|Ml) ). despite the fact that it is not an eigenstate of the quantum 
Hamiltonian (Eq. ([571) ). 



5.2.2 The Ji — J2 model. Let us now move on to the Ji — J2 model that 
can be frustrated for some values of interaction couplings, and is a paradigmatic 
example of VBS. The Hamiltonian of this model can be written as 

if J, J, = 2Ji J2 ■ + 2-^2 ■ '^3^ 

where {{ij)) denotes the next nearest neighbors. This is a frustrated model (even 
in the case of a linear chain), when both J12 7^ and at least one of them is 
positive. As will be discussed in Sec. 15. 3i this Hamiltonian can also be obtained 
using Fermi-Fermi mixtures in an optical lattice. Its importance comes from its 
usefulness in explanation of magnetic properties of Li2VOSi04 and Li2VOGe04 
compounds [4501 . Below, we will discuss the different phases of this model. 

The one- dimensional system. The Hamiltonian (Eq. fTDl) ) for the case of the 
ID lattice, with an even number of sites, periodic boundary conditions, and Ji = 
2J2 > is known as the Majumdar-Ghosh model |451j . The ground state space of 
the model is spanned by the two dimersQ 

N/2 

\^Mg)± = n (|0>2.|1)2.±1 - |l)2.|0)2.±l) /V2, (71) 
i=l 

where N is the number of sites in the lattice. We also call these singlet states 
as "valence bond states". Note that unlike the AKLT state, this model breaks 
translational symmetry. A "valence bond solid" is formed from these valence bond 
states after they order between themselves. 



^The Majumdar-Ghosh model is one of the few frustrated spin models for which the exact ground 
states are known. See also, e.g. | 452ll453ir455ll454| . 
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Figure 18. The Ji — J2 model with dimer configurations in a square lattice. The double lines 
indicate singlets between the two corresponding spins. The interaction strength corresponding to 
the diagonals of the squares are J2 , while that of the sides are Ji . 



The two dimensional model. The 2D case is far more complicated. The quantum 
phase diagram of this model is not clear even for the square lattice, while the 
classical phase diagram is quite well- understood |107j . In the latter case, it is known 
that when J2/J1 is very small, the system is Neel ordered, and in the opposite 
extreme (i.e., J2 S> Ji), it has coUinear ordering (different than Neel) characterized 

by 

i*d,..»7,)=n n n (72) 

i odd j even j 

where i and j are respectively the indices of rows and columns of the 2D lattice. 
However, when 0.4 < J2/J1 < 0.6 (the strongly frustrated regime), no such order- 
ings exist. 

These classical predictions are in qualitative agreement with semi-classical and 
fully quantum calculations based on: series expansion [45611457114581145911460] (see 
[46111462] for an introduction to this method), spin wave theory [4631I464| (first 
introduced by Anderson [465] . and then extended to higher order by Kubo |466| 
and Oguchi |467| ) , exact diagonalizations [468il469ll470j , and Quantum Monte Carlo 

In the semi-classical limit, spin wave theory predicts a first-order phase transition 
from the Neel ordered phase to the coUinear ordered phase (in the highly frustrated 
region at J2/J1 ~ 0.5). However, when quantum fluctuations are introduced in this 
region, a new phase will appear to separate these two phases. From the exact 
diagonalizations [4561I472[|473] one obtains that there are two phase transitions. 
One is at J2/J1 ~ 0.38 and another one is at J2/J1 ~ 0.6. The first one is a 
second order phase transitions from the Neel state to a spin liquid valence bond 
state. The second transition is a first order phase transition from the spin liquid 
state to a coUinear state. In the regime 0.38 < J2/J1 < 0.6, many calculations 
[171I1751I17B, 477, 478,47511150] suggest that the ground state may have the VBS 
dimer configuration with long range ordering, as shown in Fig. 1181 

Despite all these efforts, the phase diagram of this model is still not completely 
understood. In particular, many candidates for the ground states are proposed for 
this highly frustrated regime. For example, the variational approach suggests [481] 
that a ground state is a spin liquid resonating valence bond state for a spin-1/2 
system. Sushkov et al. |482| (see also j468j ) found two new second-order phase 
transitions by using series expansion method: one of them around J2/J1 = 0.34 ± 
0.04, while the other in the strongly frustrated regime (at J2/J1 = 0.50 ± 0.02). 



Figure 19. (a) (Left) The kagome lattice. Two parent triangles superimposed on each other. Each 
of the parent triangles (one of which is inverted) contains three triangles at its three corners. The 
triangles in one of the parent triangles (say, in the non-inverted one) have interaction couplings J, 
while those in the other have J', between all their vertices. The kagome lattice consists of a 2D 
array of such superimposed parent triangles, (b) (Right) The honeycomb lattice for a spin-1/2 
system, whose Hamiltonian is given in Eq. I I78I I. The three types of links are depicted in the figure. 

5.3 Heisenberg antiferromagnets and atomic Fermi- Fermi mixtures 
in kagome lattices 

In the preceding subsection, we have discussed quantum magnets either in a hnear 
chain or in a square lattice. In this section, we move to a more comphcated lattice: 
the kagome lattice. The Heisenberg antiferromagnet in this lattice is frustrated. 
A schematic diagram of the kagome lattice is depicted in Fig. [121 During the last 
15 years, extensive work have been done on this model with Ising-type nearest 
neighbor interactions |483p438j . or with Heisenberg- type nearest and next-nearest 

us a lot of information about the ground and excited states of these models. There 
are still, however, many questions to be answered. 



5.3.1 Heisenberg kagome antiferromagnets. The Hamiltonian, in this 
case, is given by 



where J and J' are the couplings (see caption of Fig. [T9)) . Both J and J' are 
positive. 

Numerical simulations of the spin-1/2 system on the kagome lattice suggest that 
the energy gap between the ground state and the lowest triplet state, if any, is 
very small (of order J/20). This gap is filled with low- lying singlets, whose number 
scales with the number of spins, N, as 1.15^ [484] . All these results suggest that 
this model may be described by the resonating valence bond states. 
In the trimerizcd limit, i.e., when the ratio J' /J is very small, Mila and Mambrini 
|4921l493j have found the number, the form, and the spectrum of singlets by using 
mean field approximation. It corresponds to short-ranged resonating valence bond 
states. This method also predicts a gap of 2 J'/3 between singlet and triplets, and 
the gap remains in the thermodynamic limit. In the fully trimerizcd limit (i.e., 
J'/ J = 0), the ground state in the subspace of short-ranged resonating valence bond 




(73) 
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states, is one singlet in each of the triangles. Although there are many theoretical 
predictions about this model, there are no clear experimental confirmations. We will 
now consider a possible way of verifying these theoretical results by using ultracold 
atoms [55] . 



5.3.2 Realization of a kagome lattice by Fermi-Fermi mixtures. Con- 
sider a Fermi-Fermi mixture at 1/2 filling for each species |63llllj . The Hamiltonian 
in this case is the spin- 1/2 Hubbard model, given by 

(ij) ' 

where n.^ = ///; (n^ = fjfi), and the operators and // (/; and //) are the 
annihilation and creation operators for the two species. Here takes the value Iq 
for intratrimer, and t'o for intertrimer hopping. In the strong coupling limit, toj ^ 
V {t~ J model) [5U] , Hpp reduces to the Heisenberg antiferromagnet Hamiltonian 
Hkag (Eq.dZSl)), where J = Atl/V, and J' = M'^ jV , and a = {a'' ,a^), with 

n-n = 2a\ P f = a"" + ia^ , P f = a"^ - ia^ . 

The total spin in the trimer takes the minimal value 1/2, and there are four de- 
generate states having tr^ = ±1/2 and left or right chirality. The spectrum of the 
system in the singlet sector consists of a narrow band of low energy states of a 
width of the order J', separated from the higher singlet (triplet) bands by a gap of 
the order 3 J/4 (2J73). 



5.4 Interacting Fermi gas in a kagome lattice: Quantum spin-liquid 
crystals 

In this section we show that an interacting Fermi gas placed in a trimerized kagome 
lattice behaves as a very special quantum magnet, a quantum spin-liquid crystal^ 
possessing both antiferromagnetic order and an exceptionally large number of low- 
energy excitations [llOlllllj . 

The trimerized kagome lattice can be created by a proper superposition of standing 
laser beams [63lllll| . Such a lattice consists of trimers (i.e., sets of three closely 
packed potential minima) arranged in a perfect triangular pattern (see Fig. l20f all. 
The ultracold fcrmions that are loaded into the lattice stay localized around poten- 
tial minima, and the Fermi-Hubbard model can be shown to capture all essential 
properties of this system. The density of fcrmions is assumed to be such that there 
are two fermions per each trimer. It is also assumed that there is a long-range 
interaction between fcrmions at nearest neighbouring sites in nearest neighbour 
trimers. This long-range interaction can be achieved, e.g., in a gas of polarized 
dipolar fermions. 
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(a) 



3 



(b) ^-V-^. 



Figure 20. (a) Enumeration of intcrtrimer (intratrimer) nearest neighbours; (b) Classical 120° 

Neel state with left chirality. 

5.4.1 The quantum magnet Hamiltonian. The spinless interacting Fermi 
gas in the trimcrised kagome lattice is described by the following extended Fermi- 
Hubbard Hamiltonian 

-f^FH = - ^ {tabflfb + h.C.) + ^ UabnaTlb, 
(a.b) (a.b) 



where (a, b) denotes nearest neighbors, a = {a, i} with a referring to intra-trimer 
indices and i numbering the trimers. The tab and Uab take the values t and U 
for intratrimer, and t' and U' for intertrimer couplings, Ua = flfa and fa is the 
fermionic annihilation operator. The sites in each trimer are enumerated as in Fig. 
[20la). We denote the 3 different intra-trimer modes by /''^ — (/i,; + /2,i + /3,i)/V3 

(zero momentum mode), and f^^ = {fi^i + z±f2,i + z'^f5s)/V^ (left and right 
chirality modes), where z± = exp(±27ri/3). 

In the limit of weak coupling between trimers, the extended Hubbard Hamiltonian 
for the problem of two fermions per trimer becomes equivalent to a quantum magnet 
on a triangular lattice with couplings that depend on the bond directions [55111101 

J N 6 
i=l j = l 

where N denotes number of trimers, J — 4[/'/9, and the nearest neighbours are 
enumerated as in Fig. ISOT a'). In Eq. (|75|1 we have cri{(f)) = cos(0)crx'' -I- sm((l))ay\ 
where the spin-1/2 operators are defined as: ctx'' = {f+^\f-^ + ./-■'V+'')/2, cr^*^ = 
_,(jWt jCOt jW)/2. The angles arc: ^ 0,^^6 - 0, 0,_.2 - <^«^3 = 27r/3, 
0i_»4 = 4'i-*5 = -27r/3, 4>i^i = 4>2^i = -27r/3, 4>3^i = 4>A^i = 0, 05^^ = 
06^i — 27r/3. The physical picture behind mapping the Fermi-Hubbard onto spin 
model is the following. There are two fermions in each trimer. In the ground state 
configuration one of them occupies the zero momentum mode, while the second 
one has to choose between either right or left chirality modes, so it stays in the 
superposition of these two modes that are further identified as spin-1/2 states. 



5.4.2 Classical analysis. We first discuss the classical theory of the model 
(Eg. ([75]) ). i.e., the large S (spin) limit. In addition to being translationally invariant. 
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0^25 -OnS 0135 -0.067 -0.071 -0.085 

0.052 -0.02 -0.023 -0.028 

Figure 21. (a) Classical 120° Neel state with right chirality. Dots and triangles show 
'^'x^ o'x^^ + '^^y' '^'y^^^ where \a\ = 1/2 and the central spin defines the x axis; (b) Spin-spin 
correlations: ((t^'^ (t^"'^'^^ + o-y'^cr^^'^'). The upper (lower) set of values corresponds to kT = 
(kT = 10~^ J/2). In both plots = 21 and i = 10 at the central site. 



the model (Eq. (j75[) ) is invariant under the point group Zg = Z3 • Z2. The generator 
of Z3 is the combined rotation of the lattice by the angle 47r/3, and of the spins by 
the angle 27r/3, while the generator of Z2 is the spin inversion in the lattice plane. 
There exist three ordered classical states with small unit cells that are compatible 
with this point-group symmetry of the model: a ferromagnetic state and two 120° 
Neel type structures with left (Fig. l^Of b)) and right (Fig.l^TJa)) chiralities formally 
defined in |110[|111] . The energies per site of these states are (i) —35*^ J/4 for 
ferromagnetic and right-handed Neel states; (ii) 3S^J/2 for left-handed Neel state. 
Hence, for J < the state with left-handed chirality will be the ground state. For 
J > the situation is more complicated since the state with right-handed chirality 
and the ferromagnetic state are degenerate ground states. 

The analysis of classical ground states can be supplemented by a numerical study 
of the 12-spin cell done by fixing the direction of every spin to ??7r/3 (n = • • • 5), 
and checking the energies of the resulting configurations. This analysis has 
revealed that for J < there are 6 ground states (Zg symmetry of (|75p ) each of 
them exhibiting the left chirality Neel order. For J > the results are dramatically 
different: there are 240 degenerate classical ground states in this case. Among them 
6 right chirality Neel states and 6 ferromagnetic states. As will be seen below, the 
large number of degenerate classical ground states for J > finds its analogue in 
a large density of low- lying excitations of the quantum version of (Eg. ((75)) ). 



5.4.3 Quantum mechanical results. The insight into quantum mechanical 
properties of the system can be obtained through exact diagonalization of the 
Hamiltonian (Eg. ((75)) ). 

Numerical findings are presented in Fig. BTT b) and Tables [1] and [21 For J > the 
ground state exhibits the 120° Neel order with right chirality. This is illustrated in 
Fig. l21f b). where the planar spin-spin correlations are presented. Direct comparison 
with the correlations of the classical state, Fig. l^Tl a). shows that the exact quantum 
correlations, although smaller, have the same order of magnitude and sign as the 
classical ones. Especially, the relative values of correlations compare nicely to the 
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V7 
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120° 


-0.125 


0.25 


-0.125 


-0.125 


0.25 


N=24 


-0.096 


0.162 


-0.083 


-0.080 


0.156 


N=21 


-0.085 


0.135 


-0.071 


-0.067 





Tabic 1. Spin-spin correlations, (cr^^'cry^ + cr^*^ cr^-^^ ) , for J > as a function of distance (expressed in 
lattice units) between sites i and j. The first row presents classical predictions for the 120° Neel state 
(compare to Fig. I2l( l . 
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V3 
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V7 


120° 


-0.125 


0.25 


-0.125 


-0.125 


N=21 


-0.134 


0.237 


-0.117 


-0.116 


N=12 


-0.137 


0.251 


-0.125 





Table 2. The same as in Table [Ij but for J < 0. 



classical result. Interestingly, the 120° Neel order survives at finite temperatures, 
as is indicated by the results obtained for kT = 10~^J/2 (Fig. I^TTb)). At such 
temperatures about 800 low energy eigenstates contribute to the correlations. To 
quantify how exceptionally dense the excitation spectrum is, we note that for = 
21 there are about 2000 (800) excited states with energies less than 0.09J (0.05J) 
above the ground state. Most of them support the spin order of the ground state so 
that antiferromagnetic order persists at finite temperatures. Moreover, numerical 
simulations show that the spectrum becomes more dense as A'^ increases. 

The exact diagonalizations do not give a definite answer of whether the gap vanishes 
in the limit of infinite lattice. What can be found out from exact diagonalizations 
is that the gap, if any, should be smaller than about 10~^J/2. The appearance 
of this very small energy scale is surprising. The smallness of the gap and the 
large density of low-energy states resemble very much the behaviour of a quantum 
spin liquid of type II. The spin liquids of type II, however, possess extremely short 
range correlations, which is in striking opposition to the behaviour of the considered 
quantum magnet. For these reasons it was proposed to name this system a quantum 
spin-liquid crystal. 

The above results for J > are in a strong contrast to those for J < 0. In the 
latter case the ground state is the standard quantum antiferromagnet with 120° 
Neel order and left chirality. The spectrum is gapped, and the classical spin-spin 
correlations approximate well the quantum ones (Table [2]). In fact, the semiclassical 
theory works remarkably well even for system sizes as small as A^ = 12. The gap is 
of the order of | J|/2 in this case, so that there are at most few states with energies 
substantially below \ J\/2 for J < 0, as opposed to the huge number for J > 0. 

As it is evidenced from the above discussion, the theoretical studies of the quan- 
tum spin liquid crystals are limited to relatively small systems. We regard future 
experiments as the best possible verification of these predictions. 
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5.5 Realization of frustrated models in cold atom/ion systems 

Here we focus on control of interactions and recent proposals based on ultracold 
bosonic or fermionic atoms , cold gases of polar molecules jll5| , and trapped 
ions [4951149611497] . Realizations of frustrated models in ultracold atomic systems 
requires either the creation of a lattice with appropriate geometry, or the control of 
the effective atomic interactions. Geometry of the lattice, as we already mentioned, 
can be engineered using superlatticc techniques (see Fig. [51] and Rcf. [lllj for the 
details in the case of kagome lattices). 




Figure 22. Scheme of the proposed experimental set-up. Each arrow depicts a wave vector of a 
standing wave laser. The three vertical planes intersect at an angle of 120° . Dark (dark blue in the 
online version) spots in the right kagome figure indicate the potential lattice minima (from 



5.5.1 Simulators of spin systems with topological order. We discuss here 
the proposal of Duan, Dcmler and Lukin, for realization of quantum magnets in 
a system of cold atoms placed in an optical lattice |64] . For simplicity, we assume 
that the lattice is filled with bosons. It is made of standing laser beams producing 
the periodic potential sin^(A:^ -r) in the direction (fc is the beam wave vector). 
The index a = {Tj i} accounts for the fact that bosons arc trapped in two different 
internal states by independent lattice potentials. 

The system is governed by the following Hubbard Hamiltonian 
Hatoms = - ^ {tp^abl^bj^ + h.c.) + ^ ^ Uan^c{ni„ - 1) + U^l X! C^S) 

{ij)<y ia i 

where {ij) denotes nearest neighbors. hi„ {b\^) is the annihilation (creation) opera- 
tor of bosons in the i-ih lattice site. The Hamiltonian (Eq . ([76]) ). describes tunneling 
of bosons in different internal states between neighbouring lattice sites, and their 
on-site interactions. 

The Hamiltonian (Eq. ([76]) ) can be cast into the quantum magnet form if the follow- 
ing conditions are satisfied (i) tf^^ <^ Ua, U^i, which typically requires large enough 
V^a] (ii) + — 1- Then, in the lowest order nontrivial approximation in 
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tf.a/Uii, one gets 

where erf = n^i - erf = 6^^ + h\^bii^, and erf = -i{h\^hi ^ ^li^n)' ^^i^ 
X _ ^mT + ^Mi ^mT X _ 

The Hamihonian fEa. ([77)) ) corresponds to anisotropic Heisenberg spin-1/2 model, 
whose different variations can be studied due to the possibiUty of adjustment of i^o-, 
Ua and C/|| couphngs. This can be achieved, e.g., via manipulation of the lattice 
potentials imposed on atoms: V^a- 

To illustrate the possibilities offered by cold atomic systems, let us consider the 
limit of t^i/t^i 0, i.e., V^-\ <^ V"^;. Then, obviously the Hatoms Hamiltonian 
approximates the nearest-neighbor Ising model (Eg. ((66)) ). Depending on the 
couplings A^z can be either isotropic or anisotropic. Moreover, they can be all either 
positive or negative depending on the U^i/U^ ratio. Therefore, as discussed in Sec. 
15.11 both ferromagnetic and antifcrromagnctic Ising model can be simulated with 
cold bosons in a lattice. Additionally, as discussed in [51], the successful experi- 
mental realization of the Ising model in a system of cold atoms can be very useful 
in the experimental implementation of the one-way quantum computer, introduced 
by Raussendorf and Briegel [498,,4gS] . 

Perhaps the most unusual quantum magnet, that can be realized using the ideas 
presented above, corresponds to the anisotropic 2D spin-1/2 model on a hexagonal 
lattice. A hexagonal lattice can be created in a carefully designed setup of laser 
beams [64] . Assuming the same as above for derivation of Eq. (|77)) , and additionally 
J7|| « J7| w C/| ~ t/ the Hamiltonian (Eq. ([75)) ) can be transformed into 

Hatoms -X. ^'^i+^I^ E ^^''^'+^- E ^'^l' (78) 

X— links y — links z — links 

where the summation goes over nearest neighbours i and j placed on the x,y,z- 
links - see Fig. [T21 The couplings read = —t^^/{2U), where is the rate 
of tunneling of the atom being in the eigenstate of the Pauli operator cr^ with 
eigenvalue +1. 

The Hamiltonian (Eg. ljTS)) ) describes the Kitaev model |500j with topological or- 
der. Excitations in this model are both Abelian and non-Abelian anyons |5011I502) . 
having exotic fractional statistics. It is exactly solvable, and possesses other ex- 
citing features |500j . Its realization should provide an exceptional possibility for 
experimental observation of Abelian and non-Abelian anyons, and could serve for 
applications in quantum information as protected qubits. 



5.5.2 Frustrated models with polar molecules. Micheli et al. [115] have 
proposed a scheme to realize spin Hamiltonians with polar molecules. In fact, re- 
cently various schemes have been proposed for trapping different states of cold 
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polar molecules. New developments in this area can be found in Ref. [146j (for 
experiments, see |503j ). 

In Ref. |115| . the authors start with two polar molecules, trapped in an optical 
lattice. The outermost shell of an electron of a heteronuclear molecule represents 
the spin-1/2 system. The total Hamiltonian of a pair of heteronuclear molecules 
trapped in an optical lattice is given by 



where Hi„t and iJe^t represent respectively the internal and external dynamics 
of molecules. H^d is the Hamiltonian for dipole-dipole interactions between two 
molecules. The rotational excitation of each molecule is described by the Hamilto- 
nian 

Hm = BN^ + jN-S, (80) 

where A'^ is the dimensionless orbital angular momentum of the nuclei and S is the 
electronic spin. Here B is the rotational constant and 7 is the spin-rotation coupling 
constant. The rotational motion of molecules is coupled by the dipole-dipole inter- 
actions that are already present in the heteronuclear molecule. To obtain strong 
dipole-dipole interactions, a microwave field is introduced. If the field polarization 
is set to the z-axis of the two molecules, and the frequency is fixed as near resonant 
with the excited state potential, it leads to some spin pattern of ground states. By 
changing the frequency and field polarization, it is possible to obtain various spin 
models, e.g., Ising, Heisenberg, and Kitaev models, that have been discussed ear- 
lier (see Table I of |115j ). The advantage of this model is the strong dipole-dipole 
interactions which is due to the inherent properties of these molecules, and is also 
due to the introduction of microwave fields. The couplings are strong relative to 
decoherence rates. 



5.5.3 Ion-based quantum simulators of spin systems. It was recently 
proposed by Porras and Cirac that different spin systems can also be simulated with 
cold ions |495[|496j . This protocol, if implemented experimentally, should allow for 
studies of different quantum magnets in ion traps. As a result, the ions could then 
be used for investigation of condensed matter problems. In fact, it seems that ions 
are perfectly suitable for such experimental studies since (i) they can be trapped 
and cooled very efficiently; (ii) the position and internal state of every single ion 
can be measured and manipulated almost at will; (iii) the external parameters of 
that system can be well controlled and changed in real time. 

Below we will briefiy present the Porras and Cirac proposal. The spin 1/2 states are 
encoded in two internal hyperfine states of each ion. The local dynamics of these 
states is governed by the Hamiltonian 

Hm ~ ^ B°'cr'^, 

3 = 1. ..N 



76 



where N is the number of ions, and are "analogs" of the magnetic field 
induced by lasers resonant with the internal transition, whereas B^ is the energy 
gap between the two internal states. 

The ions are affected by both an external harmonic trapping potential with fre- 
quencies (ri2,,f2j,,rJz)and Coulomb interactions, which is described by the potential: 

where m is the mass of the ion and e is the electron charge. It is further assumed 
that the trap has a cigar-shaped geometry with axial z direction and radial x and 
y ones: Q,z ^ Q,x,^y- The competition between harmonic squeezing and Coulomb 
repulsion results in oscillations of ions around equilibrium positions. Dynamics of 

these oscillations is governed by the vibrational Hamiltonian: iJ„ ~ 'Y^'^Zi'^n '2m ~'~ 
V , where are momentum operators of the j-th ion. After expansion of V up to 
a quadratic order in ions displacement from equilibrium positions, the vibrational 
Hamiltonian can be diagonalizcd in a standard way in terms of collective modes 
(phonons): 

OL—x^y^z 



^ a 



where „ and aj^ „ are the phonon annihilation and creation operators and oJa^n 
are the collective mode frequencies. 

Finally, one needs a coupling between Hj^ and Hy, i.e., between "effective" spins 
and phonons. This is achieved by placing the ions in off-resonant standing wave 
beams. For instance, imposing standing light waves on ions one can realize the 
following Hamiltonian |496| : 



= -J^xE^^I tXt U.j-FyY.y^\ TXT \y,j, 

3 3 

where tr"! t)^ = | t)cf This Hamiltonian assumes that the laser beams push ions in 
the upper state only, which can be achieved by a proper adjustment of the relative 
phases of the beams. 

The total Hamiltonian of the system becomes then H = H„i + Hy + Hf. The 
final step to get the clear picture of the phonon-mcdiatcd interactions between 
"effective" spins relies on the unitary transformation H — > UHU^ with U specified 
in [5041l495ll496j . The results of this transformation will be illustrated on particular 
examples. 

Assuming that ^ and Fy = 0, one gets to the lowest order in rja = 

^ ^ ^ E ^fi'^'^l + E + + (81) 



where 



-^-i ^ \{z,) ^ {z,)r ^^^^ 
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with {zi) being the equilibrium position of the i-th ion. The first two terms in 
Eq. (|8f p act in the "spin" space and exhibit the quantum Ising model type of spin- 
spin and spin- field interactions. The fact that the interactions are long-range makes 
quantum properties even richer then those of a standard Ising model with only near- 
est neighbor terms. Therefore, one gets the possibility of studying the fascinating 
Ising-type model [496| in an ion chain if the corrections coming from the perturba- 
tion He |495il496j are negligible. To quantify influence of the perturbation He on 
spin dynamics, we note that the deviations of expectation values of M-spin(site) 
observables, induced by skipping the He term, scale as Mry^. Since the observ- 
ables of interest during analysis of quantum phase transitions correspond usually 
to M = 1 or 2 classes, the corrections coming from He should not cause major 
problems. 

Another possibility shows up when there are two forces acting in radial directions 
X and y [F^, Fy 7^ 0). Then one gets that the U transformed Hamiltonian is 

^ « ^ Y.^J^^r'^i + -^S^^J) + E + + (83) 

where are given by Eq. ((82| . This time, the first two terms correspond to the 
XY model with long-range interactions. The perturbation affecting spin dynamics. 
He, is now a little different than the one above, but as long as uj^ 7^ the errors 
induced by leaving the He out are of the same order. Hence, we conclude that ion 
chains can be used for simulation of yet another important spin model. 
Though the above discussion was focused on the one-dimensional geometry of the 
ion trap, a similar approach can be used for studies of ions in two-dimensional 
geometries, where they form Coulomb crystals |505[I506] . This time, two dimen- 
sional "effective" spin models should show up. Their experimental realization would 
greatly facilitate studies of frustrated quantum magnets. 

Finally, we remark that the scheme discussed in this section offers an access to 
measurements and manipulations of arbitrarily selected ion(s): an exciting oppor- 
tunity unavailable in traditional condensed matter systems. Interestingly, since the 
"spin" configuration is encoded in ions internal states, the fluorescence signal of 
the sample can contain enough information for detection of different phases. 



6 Ultracold spinor atomic gases 
6.1 Introduction 

Interactions in bosonic systems with spin degrees of freedom host a wide variety 
of exotic phases at zero temperature and a dynamics clearly differentiated from 
the one displayed by scalar gases. Spin effects are enhanced in the limit of small 
occupation number and strong interactions. This fact, makes the study of ultra- 
cold spinors in optical lattices of primordial importance and a very valuable tool 
to deepen our understanding of magnetic ordering and condensed matter related 
issues. 

On the other hand, there exists also an increasing interest in the study of spinor 
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fermions in optical lattices since these systems could serve as quantum simula- 
tors of fermionic Hubbard model and shed some light on the problem of high 
Tc-superconductivity. we shall address this point at the end of the chapter. 
The spinor degree of freedom on alkaline gases corresponds to the manifold of 
degenerate (in the absence of an external magnetic field) Zeeman hyperfine energy 
levels. The energy levels are described by the total angular momentum F = I + J, 
where I refers to the nuclear spin and J = L + S describes the total electronic 
angular momentum. In the ground state, alkaline atoms have a single electron in 
an ns orbital, their electronic angular momentum J = S = \/2 and, therefore, they 
have two possible hyperfine of Zeeman states. The nuclear spin I depends on each 
atomic species. ^'^^Cs atoms have a nuclear spin / = 7/2, while ^^Rb atoms have 
/ = 5/2, and *^Rb and ^'^Na have / — 3/2. In the case of Cesium, the hyperfine 
manifolds correspond to = 4,3; for '^^Rb they correspond to = 3,2 and *^Rb 
and ^•^Na have manifolds F = 2, 1. Each ground state manifold consists of all 
Zeeman states associated to a given total angular momentum, i.e. mp)} where 
mp = —F,..,F. In this context, we identify the total angular momentum of the 
atom F as the atomic spin. 

If atoms are magnetically trapped, the degeneracy is broken and the atoms min- 
imize their energy in the so-called "weak field" seeking states, characterised by 
a fixed hyperfine level \F,mF)- The simultaneous confining of different "low field 
seeking" states is usually unstable against spin-spin collisions. This is why magnetic 
trapping results in a frozen atomic spin "F" and spin projection "mi?" [508] . In 
that case, the bosonic quantum field operator ^' describing the creation of a boson 
in position r is a scalar in the spinor space with no dependence whatsoever on the 
hyperfine magnetic level. In turn, the mean field description of the ultracold bosonic 
gas is done via a scalar order parameter which has not explicit dependence on the 
ground state manifold \F,mp) in which the atoms were trapped. In some cases, 
despite the fact that the levels arc not anymore degenerate, it is possible to mag- 
netically trap simultaneously more than one hyperfine component. This is the case 
of ^'^Rb, where usually the singlet and the triplet scattering lengths are practically 
equal and thus spin exchange collisions are highly suppressed. It is then possible 
to magnetically trap simultaneously the \F = 2,mp — 2) and \F = l,mp = — 1) 
states. In such cases, one speaks of "multi-component" ultracold gases, but still the 
spin of the atoms remains "frozen" , although with two different possibilities |507| . 
On the contrary, if atoms are optically trapped with a far-off resonance laser, 
then all atoms, regardless of their hyperfine level are simultaneously trapped. The 
bosonic quantum field operator in such case is no longer an scalar and has to incor- 
porate this new degree of freedom "^mp ■ In turn, the corresponding order parameter 
describing mean field approach, becomes a vector whose components correspond 
to the different accessible hyperfine levels |508[I509| and which transforms under 
rotations in the spin space as a vector preserving the symmetries present in the 
corresponding spin space. 



6.2 Spinor interactions 

Like in the scalar case, ultracold atomic spinor interactions can be parametrised by 
two-body short range (s-wave) collisions. In the most general scenario, symmetry 
arguments impose that the collisions between two identical bosons in a hyperfine 
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spin level |F, ttif) are restricted to total spin S = 2F,2F — 2,...,0. The contact 
potential can be written as: 

V= 9sPs, (84) 

S=0,2,...,2F 

where Ps {S = 0, 2, 2F) denotes the projector onto the subspace with total spin 
S. The interaction strength gs characterising the contact potential interaction of 
the different S channels are given by gs = ^^^"^ ' where as is the corresponding 
scattering length and m the atomic mass. The different values of the various as 
will lead to distinct magnetic ordering. 

To understand the ground state properties as well as the dynamics, it is convenient 
to express the interaction potential V in terms of spin operators by using different 
operator identities, e.g. / = X]5=o 2 2F and Fi • F2 = ^ ^sPs where As = 
(1/2)[S'(5' + 1) - F{F + 1)]. For spin F = 1, the total spin is S* = 2, and using 
the above identities the contact potential can be written as: 

V=i) =co + C2Fi -Fs. (85) 

The terms with coefficients co and C2 describe spin-independent and spin-dependent 
binary elastic collisions respectively in the combined symmetric channels of total 
spin and 2, and are expressed in terms of the s-wavc scattering lengths ao and 02 
as: Co = 47rft^(ao + 2a2)/3m and C2 = 47r/i^(a2 — ao)/3m ^509, 510j. 
It is straightforward to generalize some spin identities to larger spins [509j by 
noting that {Fi ■ i^2)" = J2 ^^Ps- Then the interaction potential can be rewritten 

as V = X)n=o'^"(-^i ■ -^2)" where the c„ are linear combinations of the different 
scattering lengths as- 

For F = 2, there are 3 possible channels with total spin S = 0,2,4. The contact 
potential can be expressed as 

^(F=2) =co + ciPo + C2Fi •F2, (86) 

where the coefficients are given by cg — 47r/i^(3a4 + 402) /7to and Ci = 47rft^(3a2 — 
10a2 -I- 7ao)/7m and C2 = 47r?i^(a4 — a2)/lm. The projector Pq can be further 
expressed in terms of "singlet" pair operator as we shall see later. 



6.3 F — 1 and F — 2 spinor gases: Mean field regime 

The experimental achievement of Stenger et al. |511| in trapping ^'^Na by optical 
means in 1998, triggered the study of spinor ultracold gases. A mean field approach 
to describe a F = 1 condensate was developed by Ho [509| . and independently by 
Ohmi and Mashida |512| in the same year. Koashi and Ueda |510| . and Ciobanu, 
Yip, and Ho |514| have calculated the mean field phase diagram of F = 2 spinor 
condensates. More recently, Ueda and Koashi |515| have studied mean field theory 
for F = 2 atoms in presence also of a magnetic external field. The recent success 
on condensing Chromium atoms [145] have also initiated the study of mean field 
phases for spin F = 3 |516ll517j . Recent experiments concern studies of dynamics 
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of spinor BECs in traps [52ni[52Tl[5T9ll518ll528l[522] (for the theory, see [288]), and 
in lattices |524j . 

In scalar bosonic weakly interacting systems, where a mean field approach can 
be used, the ground state is found by approximating the bosonic operator of the 
corresponding Hamiltonian by its mean, and minimizing the energy functional un- 
der the constrain that the number of particles is fixed (grand canonical ensemble) 
d{H — ^N) /dj.! ~ 0. For spinor gases, one follows the same approach |509| and min- 
imization of the energy leads to different ground states depending on the values of 
the spin-spin coupling. 

The Hamiltonian of a trapped cloud of ultracold atoms with spin _F = 1 in second 
quantization reads 

+ + y*L*]F™fe • F,i^>iM>k} , (87) 

where 'l'm(r) (^'m) is the field operator that annihilates (creates) an atom in the 
m-th hyperfine state \F = 1, mp) at point r. The external trapping potential, Vext, 
is normally assumed to be spin- independent. 

There are two distinct ground state mean field phases for spin F ^ 1: 

• Ferromagnetic phase for F ^ 1. 

The system presents ferromagnetic order if C2 < (i.e. > ^2)- This configura- 
tion minimizes the energy by imposing that (F)'^ ~ 1. For a spin pointing along 
n ~ (cos a sin /?, sin a sin (3 , cos (3) the condensate order parameter is given by: 

/1\ / e-*"cos2(/3/2) \ 

e = e^^U{a, /3) U e^W ^(2) cos(/3/2) sin(/3/2) . (88) 
\0/ \ e*"sin2(/3/2) / 

This is, for example, the ground state of ^^Rb in the F = 1 manifold. 

• Polar phase for F = 1. 

The system presents polar (sometimes called antiferromagnetic) ordering if C2 > 
0. This is the case for ^^Na or '^^Rb. Minimisation is achieved by demanding that 
the expectation values of the spin component are zero along any direction, i.e. 
(F) = 0. However, variances are not equal to zero, indicating that the system 
does not posses rotational invariance. The general expression for the spinor in 
this case reads 

/0\ /-e— sin(/?)/v/(2)\ 

e = e"^ 1 = e"^ cos(/3) . (89) 

\0/ \ e"sin(/3)y(2) / 

Notice that this "antiferromagnetic" ordering do not refer to orient the spins an- 
tiparallely as it happens classically. 

For F = 2 (for experiments see [519] ) collisions can occur in one more channel 
corresponding to total spin S' = 4. As a result there is one more possible magnetic 
ground state, the so-called cyclic ground state. Now, to characterise the different 
ground states one should consider not only magnetization, i.e. (F), but also the 
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"spin singlet pair creation" expectation value |515|,l514j . The term proportional to 
Pq in the contact potential (see Eg. ([55)1 ) can be expressed by means of a "spin 
singlet pair operator" as Pq = (2/5)5' -i-S"- , where S+ = aJaQ/2 — a\a^_^ + a\a'_2, 

and aj^ {ua) creates (annihilates) a particle with spin projection a {S- ~ s\_). The 
operator S- applied on the vacuum creates, except for normalization, two bosons in 
a spin singlet state. Ground states can be classified according to their expectation 
values of the "magnetisation" and the "singlet pair creation" : 



m— — 2,n— — 2 



(90) 



Ferromagnetic phase for F = 2. _ 
The ferromagnetic phase is achieved for (F) ^ and {S+) = 0. The system 
prefers ferromagnetic order if ci < and ci — (C2/2O) < 0. The expectation 
value of the magnetisation can either be (F) = 2 or {F) = 1 but for both cases 
(5+) = 0. A representative of maximal spin projection, ground state is given by: 



' 




Vo/ 



(91) 



This corresponds, for instance, to the ground state F ~2 oi ^^Rb. 
• Polar condensate for F = 2. 

Polar (antiferromagnetic) ordering is described by a non zero singlet amplitude 
(5+) 7^ and zero magnetisation {F) = 0. A representative of the corresponding 
spinor order parameter is given by: 




1 




(92) 



The phase space boundaries for polar phase are given by C2 < and ci — (C2/2O) > 
0. 

• Cyclic phase for F = 2. _ 

Finally, cyclic ordering appears for ci, C2 > there (F) = and (5*+) = 0. 



6.3.1 F=l gases in optical lattices. Studies of F = 1 systems have already 
been carried out by Demler's group |1181I119] . They have derived an approximate 
phase diagram in the case of antiferromagnetic interactions. As in the standard 
Bose-Hubbard model, an F = 1 spinor gas undergoes supcrfluid to Mott insulator 
transition as tunneling is decreased. In the antiferromagnetic case in 2D and 3D, 
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the SF phase is polar, and so are the Mott states with an odd number TV of atoms 
per site (those states are also termed nematic). In the case of even N, for small 
tunneling the Mott states are singlets, and for moderate tunneling there occurs a 
first order transition to the nematic state. 

In ID there is furthermore the possibility of a dimcrizcd valence bond solid state, 
as in the Majumdar- Ghosh model [SO] (discussed also in subsection I5.2.2p . This 
possibility was studied by Yip |120| , who derived an effective spin Hamiltonian for 
the MI state with = 1, and used the variational ansatz interpolating between 
dimer and nematic states to argue that in a wide range of parameters the spinor 
^■^Na lattice gas should have the dimer ground state in ID, 2D, and 3D. This 
is a very interesting result, since dimer states have not been so far observed in 
experiments. This result has been supported by rigorous studies in Ref. [121| for spin 
systems with an even number of spins described by the same effective Hamiltonian: 
it was shown that while the ground state of the system has total spin Stot = 0, 
the first excited state has Stot = 2. Yip's results were recently confirmed by Rizzi 
et al. [525j . who numerically studied the SF - MI transition in the F = \ Bose- 
Hubbard model in ID, and found that the system is always in the dimerized state 
in low tunneling regime of the first MI lobe, where the effective spin model of 
Ref. |1181I119| works. Similar results were obtained by Porras et al. |526j . Thus, 
nematic order seems to be strictly speaking absent in ID in the thermodynamic 
limit. However, susceptibility to nematic ordering grows close to the border of the 
ferromagnetic phases, indicating that it may persist in finite systems. A completely 
new insight can be gained by looking at entanglement transport properties of = 1 
chains [527] , which seem to confirm existence of nematic and trimer regions for finite 
systems. 

Another interesting aspect, namely the possibility of controlling the order of the SF 
- MI transition by using appropriately polarized (lin-0-lin) laser fields to form the 
optical lattice was investigated in Refs. J28j529j. Such a laser configuration couples 
the states with mp = ±1, so that the system becomes effectively two-component. 



6.3.2 Bose-Hubbard model for spin 1 particles. The derivation of the 
Bose-Hubbard Hamiltonian for ultracold spinor gases is performed in the same 
way as in the scalar case. One has to add to the scalar Bosc-Hubbard model the 
spin dependent part of the interaction. Following the identities given in the previous 
part of the section, the Bosc-Hubbard Hamiltonian for spin 1 particles is obtained 
straightforwardly |119| : 

<ij>,<j i 
i i 

where h^i annihilates a boson in a hyperfine state m^? = cr at site i, rii denotes the 
number of particles at site i and = X^o-o-' ^aa'ba'i is the spin operator at site i 
{T^cr' being the usual spin matrices for a spin-1 particle) and < ij > denotes pairs 
of nearest neighbours in the lattice. The first two terms in the Hamiltonian represent 
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tunneling between nearest-neighbor sites and Hubbard interactions between bosons 
on the same site, respectively, as in the standard Bose-Hubbard model. The third 
term represents the energy associated with spin configurations within lattice sites 
which penalizes non zero spin configurations in each individual lattice. This term 
will induce distinct Mott phases that differ from each other in their spin correlations 
[119j . The appearance of spin mediated tunneling transitions in the optical lattice 
depends clearly on the ratio between the different energy scales appearing on the 
Bose Hubbard Hamiltonian. We shall consider in what follows the polar superfluid 
phase characterised by C2 > 0, and assume that spin- independent interactions Cq 
are larger than the spin-dependent ones C2 |119j . The supcrfluid-Mott transition 
depends then only on the ratio t/cQ. However, once the Mott regime is achieved, if 
C2 > t one should expect different magnetic ordering in the insulator phase due to 
the spin interactions. On the contrary, if t ^ C2 tunneling will manifest equal for 
all spin components and the gas will behave as a strongly correlated scalar gas. For 
small but finite tunneling t/co <C 1 it is possible to perform perturbation theory. 
To derive an effective Hamiltonian to second order in t/cQ. We split, as usually, 

the full Hamiltonian as H = Hq + Ht where Hq = (co/2)n(n - 1) -I- {c2/2)F'^ - 
2n — i^m describes the one site unperturbed Hamiltonian and Ht describes tunneling 
between two adjacent sites, i.e we consider a two site problem. To derive the effective 
Hamiltonian one looks how the energy of the unperturbed ground states \g, S) is 
lowered due to tunneling: £5 = — \i'^\^t\g^)\ -^j^gj-g labels all the (virtual) 
intermediate states and Ei,, Eg denote the unperturbed energies of the two-site 
state \g,S) (which are non-degenerate). The dependence of the energy shifts 
on the total spin S introduces nearest-neighbor spin-spin interactions in the lattice. 
It is sufficient to evaluate these shifts for only one value of ms because tunneling 
cannot mix states with different ms and the overlaps |(i^|-fft|g, S")! are rotationally 
invariant. The effective Hamiltonian can be written then as: 

Fi+F2=S 

H^,= E 'sPs, (94) 


where, now, the sum extends to all Fi + F2 values. The effective Hamiltonian for 
S' = 1 in second order on t can be expressed as a generalized quadratic Heisenberg 
Hamiltonian [119] : 

H,, = - Jo - Ji ^ FiFj -J2J2 (Fi • Fj)^ (95) 

{i,3) (iJ) 

where the explicit expressions for Jq, Ji, and J2 can be found in Ref. [119] . 
This Hamiltonian differs from the familiar Heisenberg Hamiltonian for spin 1/2 
particles due to the presence of the quadratic term on the spin interaction. Higher 
order terms on perturbation theory can give rise to terms with higher powers in 
(FjFj) but with much smaller coefficients |120] . We will see, however, in the next 
subsection, that higher orders in the Heisenberg interactions appear in a perturba- 
tive treatment to second order if larger spins are used. We should ignore now higher 
order perturbative terms (i.e. oc ) since their contribution are highly suppressed. 
We summarise the possible quantum phases in such a case. As reported by Imam- 
bekov et al. [119] . for an optical lattice in 2D and 3D, the phase diagram can be 
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Figure 23. General phase diagram for F = I bosons in 2D and 3D optical lattices (from [119] ). 

summarized as follows. The term Ji favors ferromagnetic order, while J2 enhances 
the singlet spin configuration on each bond. To solve this competition, there is the 
possibility of having nematic states that mix states with total spin S = and 5 = 2 
(and mp = 0) at each bond, but are product states; 

\N) ^Yl\F, = l,m,^0). (96) 

i 

This is a nematic state, with zero expectation value of each spin but not rotationally 
invariant. For an odd filling factor, the Mott insulating phase is always nematic. 
For even filling factors, there is always a spin-singlet phase, in which pairs of atoms 
at a site form singlets. See Fig. [221 

Apart from nematic and singlet states, in ID, an additional solution is the dimerized 
VBS state 

\D) - n 1-^^ = 1' ^'+1 = 1. ^» + F^+i = 0), (97) 

i—2n 

in which the translational invariancc is broken. 



6.3.3 F—2 gases in optical lattices. The mean field states of spinor F ~ 2 
gases have been for the first time investigated in Refs. [5141l510[l515j . It is worth 
noticing that mean field states are also valid for insulating Mott states with one 
atom per lattice site, provided all atoms are described by the same single-particle 
wave function attached to a given site. Refs. |510ll515j go one step further, and apart 
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from the mean field theory also consider the extreme case of quenched (immobile) 
F = 2 bosons in an optical lattice. In other words, these references characterize 
possible on-site states for TV bosons with total spin S in the absence of tunneling. 
After submission of the first version of this paper, Barnett, Turner and Demler 
have presented a beautiful and complete classification of the mean field phases for 
arbitrary F, based on the 19th century F. Klein's method of solving quintic poly- 
nomials by analysis of rotations of regular icosahedra [168j . We discuss their results 
in more details in the following. Also, very recently the effective spin Hamiltonians 
(in the first MI lobe), and quantum insulating phases of F = 2 bosons have been 
studied by Zhou amd Semcnoff [169] , using variational principle applied to product 
(Guztwiller ansatz, cf. [29j). dimer and trimer states. 



6.3.4 Bose-Hubbard model for F—2 particles. The derivation of the Bose- 
Hubbard Hamiltonian for spin F = 2 can be conveniently expressed |515j as: 

<ij>,a i 

(98) 

i i 

where the : : denotes normal ordering of the operators. 

The ratios between the various interactions, ci/cq and C2/C0, are fixed by the scat- 
tering lengths whereas the ratio t/co between tunneling and Hubbard interactions 
can be tuned by changing the lattice parameters The possible ground state 
phases for F = 2 have been studied in |122j using second order perturbation theory 
with t/co <C 1. At zeroth order, the Hamiltonian is a sum of independent single-site 
Hamiltonians: 

Ho^'fnin^l)+'^:F-F:+^S+S^. (99) 

Exact eigenstates of this Hamiltonian have been obtained in Ref. [515] . Since S± 
commute with the total spin operator, the energy eigenstates can be labeled with 
four quantum numbers as \N, Ns, F)mp- They have a 2F + 1-fold degeneracy as- 
sociated with the quantum number mp and their energies are given by: 

E = ^N{N - 1) + ^{F{F + 1)- 6N) + ^Ns{2N - 2Ns + 3). (100) 
2 2 5 

Assuming one atom per site, the effective Hamiltonian to second order in the per- 
turbation parameter t/co has the following form: 



p(y) p(v) p(y) 
50 52 54 



(101) 



Notice that £5 = —At^/gs indicates that the control and engineering of the magnetic 
properties of the system could be easily achieved by means of e.g. optical Feshbach 
resonances. The Hamiltonian fEq. pOip ) can be easily generalized to the whole 
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lattice, H = Hoi + X]<y> H^j''^ . It can also be transformed into a polynomial 
of fourth order in the Heisenberg interaction Yi ■ F ^ : 



« - E i^.,. + E 55SL-^(F, ^ F,) + »,F. ^ F, 



|2 



51 ^ ' ■'^ 102 

Teo lOe^ £4 \ . 3 + ^.p + 10^4 , 4_ (^02) 
204 204 72/^ 1020 ^ ' ^ 

Very recently, effective spin Hamiltonians (in the first MI lobe) and quantum insu- 
lating phases of = 2 have been studied in [169] . 

To study the ground state phase diagram of a one-dimensional F ~ 2 spinor gas 
assuming that there is the freedom to modify the value of (the zero magnetic field) 
scattering lengths |122| . one can use Matrix Product States (MPS) methods. MPS 
represent the ground state of a translationally invariant short range Hamiltonian 
exactly |311j , or nearly exactly (see subsection I3.10.2p . A straightforward way to 
do that [530j is to add to the bond Hamiltonian (jlOip c (> 0) times the identity 

operator /^'■'^ = -Pg on the bond, so that the new Hamiltonian H' = ijj''^^ +cl 
becomes positive definite, 

Hf^ = J2^sP^s'\ (103) 

s=o 

i.e., all As are non-negative. In particular, Ai = A3 = c and Aq = c — At^ / Qq, 
A2 = c — 4t^/(?2 and A4 = c — 4,t^ /g^. Since by definition, Ai = A3 are the largest 
parameters, the ground state will obviously belong to the symmetric subspace. 
Before we proceed, it is worth to discuss the mean field phase diagram obtained 
under the assumption that the ground state is a product state, |^') = |e,e,...) 
(see [51411168] ). There are 3 possible mean field (i.e., product) ground states, with 
|e) given (up to 50(3) rotations) by: 

• Ferromagnetic state, |e) = (1,0,0,0,0); possesses only the U{1) symmetry of 
rotations around the z-axis. and has maximal projection of the spin onto z axis. 

• nematic state, has the ry-degeneracy, |e) = (sin(?7)/%/2, 0, cos(?7), 0, sin(r/)/-\/2). 
This state is a MI version of the polar state in BEC; it has mean value of 
all components of the spin equal zero, but non vanishing singlet projection 
(singlet|e, e) 7^ 0; 

• Tetrahedratic (cyclic) state, |e) = (l/VS, 0,0, ^2/3, 0); this is a MI version of 
the cyclic state. The state may be uni- or biaxial, depending on whether the 
nematic tensor does, or does not have a pair of degenerated eigenvalues; it has 
vanishing of both, of mean values of all of the spin components, and of the singlet 
projection. 

The mean field phase diagram corresponds to: 

• a ferromagnetic state for A4 = 0, A2,Ao > 0, and for Aq = 0, provided A2 > 
I7A4/IO; 

• a nematic state for Aq = 0, provided 3A4/IOA2 < I7A4/IO; 

• a cyclic state for Aq = 0, provided A2 < 3A4/IO, and for A2 = 0. 

In ID, since the Hamiltonian is a sum of nearest neighbor bond Hamiltonians, 
we have X^fc fe' (^fc' ^fe • ■ ■ \^i\^k> ,Ck' ■ ■ ■) esc (efc|efe/)^~^ and thus in the limit of an 
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infinite chain the ground states are equally well described by product states (that 
will typically break the rotational symmetry). This means in this case we expect 
mean field product states to provide a very good approximation of the ground 
states with translational symmetry. Below we present a schematic classification of 
the possible ground states for the specific values |122] . 

(A) For A4 = A2 = Aq = 0, all symmetric states are ground states, i.e., in particular 
all product states |e, e . . .) with arbitrary |e). 

(B) For A4 = A2 = 0, Aq > 0, the ground states |e, e . . .) remind the cyclic states 
states of Ref. |118) (i.e., they correspond to translationally, but nor rotationally 
invariant product states), which now mix S = 2 and S = A contributions on 
each bond, and they have to fulfill the condition (singletje, e) = 0. Denoting by 
|e) = (e2, ei, cq, e_i, 6-2), this implies Cq — 2eie_i + 2e2e_2 = 0. These states form 
a much greater class than the cyclic ones, since they may have non-vanishing (and 
even maximal) components of the spin. Interestingly, the transition between the 
cyclic phase for A2 = 0, and the ferromagnetic phase for A4 = 0, occurs via such 
states, i.e., at the transition point the degeneracy of the ground states manifold 
explodes. 

(C) For A4 ~ and A2, Aq > 0, the ground states are ferromagnetic 
states |2)n|2)„ • • • |2)„, corresponding to a maximal projection of the local spin 
onto a given direction n = (sin(6') cos((/)), sin(6') sin(0), 003(6*)). Such vectors for 
F = 2 may be parametrized (in the basis of Fn with descending mp) as 
(X (z^^, 2z~^, V6, 2z, z^) with z = Izje**^, \z\ £ (—00,00). It should be stressed 
that ferromagnetic states are exact ground state in the entire part of the phase 
diagram whenever A4 = 0. 

(D) For Ao = and A4, A2 > 0, the ground states apparently favor antiferromag- 
netic order. This, however, can be misleading if A4 <C A2. In that the mean 
field diagram suggests, the ferromagnetic order might prevail. We have applied in 
ID a more general variational approach, going beyond mean field. We have looked 
for ground states by applying the variational principle to mean field (product) states 
|e, e . . .), Neel-type states |e, /, e, / . . .), and valence bond solid states with singlet 
states for distinct pairs (dimers) of neighboring atoms and translational dimer sym- 
metry. For the mean field case as discussed earlier the energy is either minimized by 
the ferromagnetic state |e) = |2)n (for A2 > I7A4/IO), by a nematic state |e) = |0)n 
(for 3A4/IOA2 < I7A4/IO; in this case the state is a combination of total spin 0, 2 
and 4), or, for A2 < 3A4/10, by a cyclic state, |e) = (e2, ei, cq, e„i, e_2) with 
eo ~ ^/V2, 62 = — e_2 = 1/2, ei = e-i = 0. Imposing Neel order with (e|/) 7^ 1 
always results in a larger energy, as Ai^s > A2,4, and the overlap with the singlet 
can be maximized already by restricting to product states. On the other hand, for 
the dimer state the energy per bond is given by ^Tr(iJ/^l (g) 1). One can apply 
MPS codes [122j to search numerically for the exact ground states in ID using the 
method of |311j to confirm the existence and phase borders of the ferromagnetic 
ground state, nematic and dimer regions. This particular phase diagram in the 
A2 , A4 phase space is displayed in Fig. [24l 

{E) For A2 = and A4, Aq > 0, as in the (D) case, mean field cyclic states are 
favorable over Neel states. One can compare them variationally to the analogues of 
the dimer states in the present case, i.e., configurations which have a state with total 
spin S" = 2 on distinct bonds. We call these state para-dimers. Now, contrary to the 
dimerized states in (D), where the states on the bond are unique, here states with 
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Figure 24. (Color online.) Sketch of the phase diagram, obtained by applying the variational 
principle in the A2, A4 phase space (for Aq = 0) to mean field, Neel, and dimer states with one 
atom per site. The scale is set by letting Ai = A3 = 1. Neel-type states are never favorable over 
nematic states. The ferromagnetic region (gray) was obtained numerically by imaginary time 
evolution of MPS (and comparing the results from runs with D = I and D = 5) in a chain of 50 
sites with open boundary conditions 13001 . (Note that here D is the dimension of the ancillary 
system that is used in the MPS method.) Of course, on the line (A4 = 0, A2) ferromagnetic states 
give always ground states. Dashed lines indicate the regions where the type of mean field state with 
lowest energy changes qualitatively. The red (dashed-dotted) line gives the values of (A2, A4) which 

can be obtained by changing the spin-independent scattering length co = {3g4 + 4g2)/7 through 
optical Feshbach resonances. [17511177] . Black and white circles indicate a change of co of 10% and 

100%, respectively. 



different ms=2 can form superpositions. It is easy to see, fiowever, tliat the best 
superposition should have maximal possible entanglement, in order to minimize 
interaction on the bonds that are not covered by the para-dimers |122] . 



6.3.5 Spinor Fermi gases in optical lattices. We move now to fermionic 
spinor gases. For F = 1/2, the quantum simulator of the fermionic Hubbard model 
has been proposed in [90] and recently studied with 3-component fermions in |531] . 

Liu et al. [H] proposed to use fermions with high F to realize spin-dependent Hub- 
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bard models, in which hopping parameters are spin-dependent. Such models lead 
to exotic kinds of superfluidity, such as to a phase in which SF and normal com- 
ponent coexist at zero temperature. W. Hofstetter and collaborators have written 
a series of papers, reviewed in [532] , on fermionic atoms with SU {N) symmetry in 
optical lattices. Such systems also have exotic superfluid and flavor-ordered ground 
states, and exhibit a very rich behavior in the presence of disorder. It is, of course, 
inevitable to ask which atoms can be used to realize high F fermonic spinor gases 
in optical lattices. The most commonly used alkali ^Li has hyperfine manifolds with 
F = 1/2 and F = 3/2. The latter is, obviously, a subject of two body losses, but as 
in the case of the F = 2 manifold of Rubidium, one can expect reasonably long life 
time in the lattice (especially in MI states with N = 1). Another commonly used 
fcrmion is a heavy alkali '^''K, which has manifolds F — 7/2 and F ~ 9/2. These 
fermions are particularly useful for spin-dependent Hubbard models [41| . There are 
several atoms whose lowest hyperfine manifold is = 3/2, i.e., in those ground 
states two body losses can be avoided: ^Be, ^^^Cs, or ^^^Ba, but so far only the 
bosonic Cesium BEC has been achieved |533] . On the other hand, recently a BEC 
of ^^^Yb atoms [534] . as well as degenerate gas of ^^^Yb fermions with F — 5/2, 
has been realized. Finally, fermionic Chromium has 4 hyperfine manifolds with 
F = 9/2, 7/2, 5/2,3/2, in the ascending order of energies, and after achieving BEC 
of the bosonic Chromium |145j . the prospect for achieving ultracold degenerate 
spinor Fermi gases are very good. 

Recently, there has been a lot of progress in understanding the special properties 
of F = 3/2 and F = 5/2 lattice Fermi gases [HllSipiMS]- C. Wu et al. [557] 
realized that the spin-3/2 fcrmion models with contact interactions have a generic 
SO{5) symmetry without any fine-tuning of parameters, and employed this fact to 
form a quantum Monte Carlo algorithm free of the sign-problem and to study novel 
competing orders in 1-dimensional optical traps and lattices |538] . In particular, 
the quartctting phase, a four-fermion counterpart of Cooper pairing, exists in a 
large portion of the phase diagram. The bosonisation approach was applied to ID 
systems with F = 3/2, 5/2, ... by Lecheminant et al. |207| . and exact Bethe ansatz 
in [208] . An overview of hidden symmetries and competing orders in spin 3/2 gases 
can be found in the excellent paper of C. Wu [539j . 



7 Ultracold atomic gases in "artificial" magnetic fields 
7.1 Introduction - Rapidly rotating ultracold gases 

Quantum systems in magnetic fields exhibit particularly interesting behavior, with 
fractional quantum Hall effect (FQHE) being a paradigm |540j . It is very well 
known that rapidly rotating trapped ultracold gases provide a possibility to realize 
"artificial" (standard, or better to say, Abelian) magnetic fields. 

In order to achieve such magnetic fields, one typically considers a quasi single 
component ( "spinless" ) 2D gas of N atoms in the XY plane rotating around the 
Z axis with frequency f2, and confined in a harmonic trap of frequency U!±. If 
the rotation is moderate f2 < in macroscopic atomic clouds, the Abrikosov 
vortex lattice is formed [541ll542j . As fl approaches uj±, the vortex lattice melts, 
and the system evolves through a sequence of states that have been identified in 
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literature [543|,ll25lll27| as highly correlated quantum liquids. 
The various regimes of rapidly rotating gases can be described in the terminology 
of fractional quantum Hall effect (FQHE) theory [540j . The crucial role is played 
by the direct analog of the Landau level filling factor in the FQHE, which can be 
related to the number of vortices Ny hy ly ~ N/Ny as defined in the BEC mean 
field description, which is valid for large systems and moderate rotations. 
In the first literature |125lll27j , the authors considered only the lowest Landau level 
(LLL) for strong enough rotation. Later, attention was paid to edge excitations and 
topological order [544| . Recently, correlated liquids at f = k/2 for k = 1,2,3, .. . 
for < — 6 — 10 have been discussed [T26.,545..5461l547j . These states resemble, 
to a great extent, the states from the Rezayi-Read (RR) hierarchy [548]: fc = 1 is 
the Laughlin state, fc = 2 is the Moore- Read paired state |5491l544j . etc. Rezayi et 
al. [550] have recently shown that the presence of a small amount of dipole-dipole 
interactions unambiguously makes the RR state with fc = 3 the ground state. This 
state is particularly interesting, since its excitations are both fractional, and non- 
Abelian. The validity of the LLL approximation for rotating gases is also discussed 
in |551| . Recently, a composite fermion (boson + one flux quantum) theory of 
rapidly rotating bosons has been formulated |552| . 

Most of the literature on ultracold rapidly rotating gases aims at considering rela- 
tively large systems and even the thermodynamic limit. Thus, in numerical simula- 
tions, cither periodic (torus) or spherical boundary conditions are used. Stability of 
the cloud in an harmonic trap requires il < oj±, since otherwise centrifugal forces 
drag the atoms away from the trap. Observation of the Laughlin states requires, 
on the other hand, to remain in the LLL. This in turn can be assured only if 
J7 — uj± = 0{1/N), i.e., it requires a very precise control of the delicate balance 
between D, and uj_i_. Unfortunately, in the case of contact (short range Van der 
Waals) interactions pseudo-hole excitations of the Laughlin state have vanishing 
interaction energy, similarly as the Laughlin state itself; they can differ only by 
angular momentum contribution (x O — — 0{1/N), i.e. they vanish at large 
N. Despite the progress in experimental studies of vortex lattices |5531l554j . and 
first steps toward reaching the LLL physics |555ll556ll557j . experiments have not 
yet reached this regime. 

The problems related to the short range nature of the Van der Waals forces can 
be overcome by using optical traps more stiff than harmonic, which is technically 
difficult, but in principle possible. Yet another promising idea is to use dipolar 
gases, i.e. gases that interact via magnetic or electric dipole moments (for a review 
see |144) ). Rotating dipolar bosonic gases are expected to exhibit exotic behaviour 
in the weakly interacting regime [149| , whereas fermionic dipolar gases have a finite 
gap for the = 1/3 Laughlin state [150] . The first observation of BEC of a dipolar 
gas of Chromium atoms with large magnetic dipole has been recently reported 
|179| . and several groups are trying to achieve an ultracold gas of heteronuclear 
molecules with large electric dipole moments j558j . The recent proposal of using 
polar molecules excited to the first rotational state is also very promising [115j . 
Completely different approach to create "artificial" magnetic fields employ the effect 
of electromagnetically induced transparency (EIT) and "slow light" ; it has been 
proposed by Juzeliunas and Ohberg [559II560II561II562] . This method is very flexible, 
and can be used to reach the integer quantum Hall effect [563] . and even regimes 
of FQHE in planar geometries [562] , or generate non- Abelian gauge fields |136| . 
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Nevertheless, atomic gases in optical lattices, or arrays of optical traps provide 
perhaps the best opportunity for FQHE. There are two ways of reaching the regime 
of Laughlin liquids. A very promising one, consists in using atoms in an array of 
rotating optical microtraps, either in an optical lattice [564| . or created by an array 
of rotating traps generated by an array of microlenses [565] , or an array of laser 
beams |566| . 

In such arrangements, it will be natural to study independent mcsoscopic or even 
microscopic systems of few atoms at each lattice site, i.e., avoid the problems of 
large N. The lattice will play, however, an important role. First of all, one can pre- 
pare equal number of atoms at each site using the Mott transition, and then turn 
the intersite hopping off. Second, the measurement and detection procedures could 
automatically be applied to many copies of the same meso-, or micro-system. Such 
experiments demand careful theoretical studies of few atom systems using possi- 
bly exact methods, such as exact diagonalisations of the Hamiltonian with open 
boundary conditions in the presence of the harmonic trap, or even a deformed trap. 
Such studies have been initiated recently. Popp et al. |128| analyzed in detail the 
possibilities of an adiabatic path to fractional quantum Hall states of few bosonic 
atoms. The papers of Barberan et al. |129| discuss the appearance of ordered struc- 
tures (vortices, vortex arrays, Wigner crystals) in systems of few bosonic atoms by 
looking at density and pair-correlation functions. In Fig. 1251 an example of such 
structures pair correlation function in the Laughlin states is shown. Similar studies 
(going even beyond LLL approximation) for bosonic and fcrmionic atoms have been 
reported in Ref. [567j . 

Another way to create highly-correlated liquids could be to mimic effects of mag- 
netic fields not by rotation, but by appropriately designed control of tunnelling in 
optical lattices. This is discussed in the next subsection. 



7.2 Lattice gases in "artificial" Abelian magnetic fields 

In order to create an effect of constant magnetic field in a 2D square lattice one 
needs to be able to control the phases of the hopping matrix elements. Single atom 
stationary Schrodinger equation, analyzed in the famous paper by Hofstadter [568j , 
reads: 

- ie"'^^= V(a;+a, y) - te'^^"V(a;-a, y) (104) 
-te^'T^^y'i{j{x, y+a) -~ ie'^'^"i/'(a;, y-a)^Ei\j{x, y), 

where ipi^^ y) is the wave function, a is the lattice constant, and A is the vector 
potential. The choice of A determines the magnetic field B, which in turn deter- 
mines the behavior of the system. For uniform field in the z direction, B = Be^, 
one may choose the potential A = (0, Bx, 0), and tunnelings in the y-direction 
acquire phases. This makes the problem effectively one-dimensional and Eq. p04p 
can be transformed into the Harper's equation [569j : 

g(m^X) — g(ra—\) -|- 2 cos(27rmQ! — v)g(m) ~ sg{m) (105) 

by using the ansatz ip{ma, no) = e"'"f/(m), where x = ma, y = na, and s = —E/t. 
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Figure 25. 3D (left) and eontour (right) plot of the pair correlation function p{r, ro) for the 
Laughlin states of A'' =6,7 and 8 atoms of mass m in a lattice site; ro = \/7V in harmonic oscillator 

units and A = y'H/2majX (from 11291 ). 

The eigenvalue problem for rational values of the magnetic flux a = per 
elementary plaquette becomes periodic. This results in a band spectrum whose 
bands form the famous Hofstadter butterfly (see Fig. Note that the regime of 
this spectrum requires finite values of a, i.e., magnetic fields B ^ which in 

the continuum limit a ^ become ultra-intense. 

Jaksch and ZoUer |130j , were the first who have recently proposed methods to realize 
such "artificial" magnetic field effects in lattice gases. Their method employs atoms 
with two internal states trapped in two lattices, laser assisted tunneling, lattice 
tilting (achieved by acceleration, or application of external static electric fields), 
and other experimentally accessible techniques. We discuss the generalization of 
this proposal to the case of non-Abelian fields in the next subsection. 
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Figure 26. The famous Hofstadtcr butterfly: energy bands plotted for rational values of the 

magnetic flux a. 



The physical features in "artificial" magnetic fields arc, as expected extremely rich. 
For instance, one may load the lattice with a cold BEC, or add disorder, and observe 
how the single atom properties of the spectrum changes in weakly interacting or 
weakly disordered systems. The modifications of the butterfly due to interactions 
and disorder, or both respectively, can be measured, as proposed in |130| . 

Several other groups made alternative proposals for "artificial electromagnetism 
for neutral atoms". E.J. Mueller jl31| generalized the method of |130j to atoms 
with 3 internal states in three lattices. The advantage of this approach is that it 
does not use lattice tilting, but on expense of more complicated laser configuration. 
Since the fields are "artificial" they do not even need to fulfill Maxwell equations, 
and can thus lead to physical realisations of otherwise unphysical effects. Mueller 
proposes for instance time dependent realisations of hopping matrix element on a 
ID ring that leads to "Escher staircase" : a single particle in such ring undergoes 
acceleration limited only by the Umklapp process: when the de Broglie wavelength 
of the particle becomes equal to the lattice constant, the matter wave is Bragg 
reflected off the lattice and reverses its direction. 

Sorensen et al. [132] have proposed yet another method that employs also time de- 
pendent hopping matrix elements along with a large oscillating quadrupolar poten- 
tial. These authors performed exact diagonalisations of the resulting Bose-Hubbard 
model in the limit of hard core bosons for N up to 5. In the limit of small a and 
low densities, the continuum limit may be applied and the model reduces then to 
a system of bosons in a magnetic field with contact interactions. For filling factor 
1/2 (total angular momentum L = 27V), the Laughlin wave function is then the 
exact ground state [448] . The numerical calculations show that the exact ground 
state remains very similar to the Laughlin state for a < 0.3. Preparation of such 
Laughlin state in the lattice may be achieved by creating first a Mott state with 
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one atom per site in a superlattice with the larger period. In this way low density 
is achieved, while the atoms are quenched and not affected by the "magnetic" field. 
By turning off the superlattice, the Mott state melts into a Laughlin liquid. 
Very recently, Palmer and Jaksch studied high field FQHE in optical lattices. They 
considered the value of a ~ ac = l/n, where n are small integers, and derived the 
corresponding effective Hamiltonian in the continuum limit, which for etc = 1/2 
reduces to a model similar to a bilayer FQHE system [570] . The corresponding 
ground state in homogeneous systems (in the absence of trapping potential) is, the, 
so called, 221 state, constructed as the Laughlin state for v = \/2 for particles 
from the same layer, and for v = \ for particles from different layers. Denoting the 
coordinates by Zi, Sj, the wave function for N = 2M atoms is 

. . . , ZM, Si, ... , sm) oc ]^(zt - Zj)^ ~ n^-^' ~ ^-J')- 

The angular momentum of this state is L = 2M{M — 1) + Af^, and filling factor 
ly = 2/3. Interestingly, the Hall current in such states in the presence of the trap 
might exhibit unexpected sign changes. 

Bukov and Demler [571j have recently studied "vortices" in a dice lattice, induced 
by "magnetic field" and have shown fascinating possibilities of creation of vortex- 
Peierls state (a bosonic analog of valence bond solids). Rotating optical lattices 
has been considered by several authors: Polini et al. |572j proposed a realization 
of a fully frustrated XY model with cold atoms in optical lattices, while Bhat et 
0.1- 1573] studied ground states of various symmetries for rotating 4x4 lattices. The 
somewhat related problem of vortex configuration in systems of rotating ultracold 
atoms in optical lattices has been studied by Wu et al. |574j , who finds that near the 
superfluid-Mott insulator transition, the vortex core has a tendency to approach 
the Mott insulator state. 



7.3 Lattice gases in "artificial" non-Abelian magnetic fields 

In a recent Letter |135| , it has been shown that by using atoms with more internal 
states, the application of state dependent, laser assisted tunneling (cf. [H]), and 
coherent transfer between internal states, one can generalize results of Jaksch and 
Zoller and create "artificial external magnetic fields" corresponding to non-Abelian 
U{n), SU{n), or even GL{n) gauge fields. In this case, the tunnehng amphtudes 
are replaced by (unitary) matrices whose product around a plaquette is non-trivial 
and its mean trace (Wilson loop) is not equal to n |134[|575j . 

We consider an atomic gas in a 3D optical lattice and assume that tunneling is 
completely suppressed in the z-direction, so that, effectively, we deal with an array 
of 2D lattice gases and we are able to restrict ourselves to one copy. The atoms 
occupy two internal hyperfine states \g)^ |e), and the optical potential traps them 
in the state \g) and |e) in every second column, i.e., for the y coordinate equal to 
. . . , n — 1, n + 1, . . . (. . . , n, n -I- 2, . . .). The resulting 2D-lattice has thus the spacing 
A/2 (A/4) in the x- (y-) direction. The tunneling rates in the x direction are due to 
kinetic energy; they are spatially homogeneous and assumed to be equal for both 
hyperfine states. The lattice is tilted in the y-direction, which introduces an energy 
shift A between neighboring columns. Tilting can be achieved by accelerating the 
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lattice, or by placing it in a static electric field. By doing this, standard tunneling 
rates due to kinetic energy are suppressed in the y direction. Instead, tunneling 
is laser assisted, and driven by two pairs of lasers resonant for Raman transitions 
between \g) and |e), i.e., n n± 1. This can be achieved because the offset energy 
for both transitions is different and equals ±A. Detunings of the lasers are chosen 
in such a way that the effect of tilting is cancelled in the rotating frame of reference. 
The lasers generate running waves in the ice-direction, so that the corresponding 
tunnehng rates acquire local phases exp{±iqx). 

In order to realize "artificial" non-Abelian fields in a similar scheme, one may use 
atoms with degenerate Zeeman sublevels in the hypcrfinc ground state manifolds, 
\gi), and \ei) with i = l,...,n, whose degeneracy is lifted in external magnetic 
fields. These states may be thought of as "colors" of the gauge fields. Promising 
fermionic candidates with these properties are heavy Alkali atoms, for instance, "^^K 
atoms in states = 9/2, mp = 9/2,7/2, . . ., and F = 7/2, m^^ = -7/2,-5/2, . . .; in 
particular, they allow for realizing "spin" dependent lattice potentials and hopping 

m- 

Having identified the "colors" , one modifies the scheme of Ref. |130| : laser assisted 
tunneling rates along the y-axis should depend on the internal state, although not 
necessarily in the sense of Ref. (41]. For a given link \gi) to \ei), tunneling should 
be described by a non-trivial unitary matrix Uy (x) being a member of the "color" 
group {U{n), SU{n), GL{n) etc.). For unitary groups, the tunneling matrix Uy{x) 
can be represented as exp(ia^j^(a;)). Here, a is real, and Ayix) is a Hermitian 
matrix from the gauge algebra, e.g., U{n) or SU{n). Since transitions from \gi) 
to \ei) correspond to different frequencies for each i, they are driven by different 
running wave lasers, and may attain different phase factors exp(±i(7ix). 
In order to create gauge potentials that cannot simply be reduced to two indepen- 
dent Abelian components, tunneling in the x-direction should be described by a 
tunneling matrix Ux, which fulfills t \Ux, Uy{x)\ ^ 0, so that a genuine non-Abelian 
character of the fields is assured. We stress that all elements of this scheme, as 
shown in Fig. 1271 are experimentally accessible. 

For the specific gauge fields considered in Ref. [135] one obtains the following U (2) 
generalization of the Harper's equation: 

axgim+l) — <7xg{m—l) + 2cos(27rmQ; — i')g(m) — sg(rn), (107) 

where g{m) is the two component wave function obtained by using the ansatz 
'ip[m.a, no) = e"^"g(77i), with x = ma, y = na, and e ~ —E/t; Ux is the Pauli 
matrix, and the "magnetic flux" matrix a — diag[Q!i, a2]. 

Given each ai = Pi/qi rational, the problem is Q-periodic (where Q equals the 
smallest common multiple of qi and (72)- The spectrum shows a band structure, 
and is bounded by two hyperplanes (Fig. I28p . It exhibits a very complex forma- 
tion of holes of finite measure and various sizes, which we name the Hofstadter 
"moth". Although a rigorous proof can not be provided, the "moth" reminds a 
fractal structure. Obviously, this fractal structure will be very sensitive to any sort 
of perturbation (finite size of the system, external trapping potential etc.) on very 
small scales. But, since the holes are true 3D objects with finite volume, the spec- 
trum will be more robust on a larger scale to perturbations than in case of the 
Hofstadter "butterfly". Very recently, the same gauge fields were used to inves- 
tigate metal-insulator transition for cold atoms |576j and integer quantum Hall 
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Figure 27. (Color online, see 11351 ). Optical lattice setup for U(2) gauge fields: Red and blue open 
semi-circles (closed semi-circles) denote atoms in states and |g2)i respectively (|ei) and leg)). 

Left) Hopping in the x-direction is laser assisted and allows for unitary exchange of colors; it is 
described by the same unitary hopping matrix Ux for both \gi) and \ei) states. Hopping along the 

j/-direction is also laser assisted and attains "spin dependent" phase factors. Right) Trapping 
potential in j/-direction. Adjacent sites are set off by an energy A due to the lattice acceleration, or 
a static inhomogeneous electric field. The lasers fiii are resonant for transitions \gii) <-> |e2i), while 
Q.2i are resonant for transitions between \eii) ■!-> \g2i) due to the offset of the lattice sites. Because 
of the spatial dependence of f2i 2 (running waves in ±x direction) the atoms hopping around the 
plaquette get the unitary transformation U = Uy{m)UxUy(m + l)f/J, where 
Uy(m) = exp(27rim diag[oi, ag]), as indicated in the left figure (from |135j ). 



effect [577) . 

The method presented in Ref. [135j can be easily generahzed to external gauge 
potentials of the form 

Ai{x,y) ai + bi{x/a) + Ci{y/a), 

with i = x,y, and ai,bi, Ci being essentially arbitrary n x n matrices. Furthermore, 
local disorder may be introduced in a controlled way that allows for small fluctu- 
ations of the matrices 0^,6^,0^. In particular, disorder can be made annealed, i.e. 
changing on a time scale comparable with the relevant time scales of the system, 
and, thus, mimics thermal fluctuations. It can be of significant amplitude, provided 
it does not drive the assisting lasers system out of resonance. Also, more compli- 
cated spatial dependences, e.g., piecewise linearity, of A are feasible by using static 
electric fields, laser induced potentials, etc. Additional lasers may introduce local, 
and in general time dependent unitaries. Such transformations would generate ar- 
bitrary local temporal components of the gauge potential, ^0(2^1 u)- Although for 
Yang-Mills fields in (2-|-l)D, this component may be gauged out adapting the Weyl, 
or strict temporal gauge |578| , the corresponding gauge transformations may intro- 
duce more complex spatial and temporal forms of the remaining two components 
of A 

In the limit of weak fields, the continuum limit may be used, and the single particle 
Hamiltonian reduces to that of a particle in the gauge field corresponding to the 
above potential, i.e.: 

The non-Abelian gauge fields and Hofstadter "moth" are interesting not only for 
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Figure 28. (Color online). The Hofstadter "moth" spectrum. Forbidden eigenenergies e are plotted 
versus ai = Pi/qi, S [0, = 0.5] (i = 1, 2), where qi < 41 and ai a2 (from 11351 ). 



fundamental reasons: they offer also a fascinating possibility of observing the non- 
Abelian Aharonov-Bohm effect and realizing non-Abelian atom interferometry and 
non-Abelian nonlinear atom optics |135ll579j . Phase shifts in non-Abelian interfer- 
ometers will correspond to matrices, and will depend not only on the perturbations 
of the trajectories of the interfering particles, but also on the specific locations of 
these perturbations. 

Another fascinating possibility concerns the possibility of realizing novel FQHE 
states, that would be associated with fractional non-Abelian excitations. For in- 
stance, moving of a pseudo-hole around origin should induce a non-Abelian Berry 
"phase" described by a non-trivial matrix [580] . Particularly interesting in this con- 
text would be generalizations of methods used in [133j . The studies of non-Abelian 
FQHE in such fields has just began. The resulting model |174j is related to, yet very 
different from, the bilayer of FQHE with tunneling between the layers |442[|570j . 

We expect that pseudo-hole in such model will transform according to the non- 
Abelian representations of the permutation group, or more precisely Artin's braid 
group, which in 2D considers particle exchange along the topologically distinct 
paths that avoid other particles (for a pedagogical review see |581) . for recent 
discussion of wave functions of non-Abelions see 15821 ). There has been a long 
lasting quest for non-Abelian anyons in the condensed matter literature, but no 
clear experimental observation so far. The most prominent candidates are electronic 
(fermionic) = 5/2 FQHE state. This state has been observed in experiments, and 
Moore and Read |549j (see also |583j '). and independently Greiter, Wen and Wilczek 
|584| proposed to explain it in terms of the "Pfafhan state" (see also Das Sarma 
and Pinczuk [585j ). Recently, however, Toke and Jain |586j proposed an alternative 
"composite fermions" model of the i/ = 5/2 state, which does not relate to non- 
Abelian statistics in any obvious manner. For bosons, a promising candidate is the 
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= 3/2 state from the, so eallcd, Read-Rezayi sequence of incompressible correlated 
liquids. This state seem to be a true ground state for the rapidly rotating gas of 
bosons interacting via contact (Van der Waals) forces with a moderate amount 
of dipolar interactions |550| . Such situation may be achieved, for instance, with 
Bose condensed Chromium, as in experiments of T. Pfau's group |145| . We hope 
that non-Abclian FQHE, due its profound and direct non-commutative character, 
will provide further, experimentally feasible examples of non-Abelian anyons. An 
alternative way to realize Pfaffian-like states and non-Abelions in ID ultracold 
gases has been very recently proposed in Ref . |587| . 



7.4 Ultracold gases and lattice gauge theories 

Another fascinating question concerns the possibility of using ultracold atoms for 
simulations of lattice gauge theories in (2-|-l)D. The main difference is that in lat- 
tice gauge theories (LGT) gauge fields are dynamical variables, whereas they are 
obviously not in the scheme of Ref. |135| . Moreover, the scheme is realized in real 
rather than in imaginary time. Nevertheless, the big advantage of the proposal is 
that given a gauge field configuration, the dynamics of matter fields in real time 
are given for free. By generating various configurations of gauge fields, we may try 
to "mimic" the Monte Carlo sampling of LGT in the limit in which gauge fields 
affect the matter fields, but not vice versa. Averaging over both, annealed disorder 
and quantum fiuctuations should approximate the statistical average in LGT. Such 
approach requires that generated configurations represent the characteristic or sta- 
tistically relevant ones of corresponding LGT phases. For instance, configurations 
in the confinement sector should exhibit an area law fulfilled by Wilson loops, ap- 
propriate distributions of centre vortices, Abelian magnetic monopoles, instantons, 
merons, calorons, etc. (for a recent reviews, see |588.,589j : see also [590j ). Although 
the gauge fields accessible in the scheme of Ref. |135j are limited, at least some of 
them share characteristics with LGT phases. 

The situation is much better when we turn to Abelian lattice gauge theories, since 
they reduce in some limits to models with ring exchange interactions, that involve 
product of operators over an elementary plaquette of the underlying lattice. The 
first proposal for such interactions has been formulated in Ref. |137j . The authors 
derived a spin model for a multicomponcnt Bose, or Fermi gas in the Mott limit 
with one particle per site. Due to the fact that in a triangular lattice tunneling 
from one site to another may occur directly, or via the remaining third site, they 
obtain 3-spin interactions in the third order of t/U expansion. This proposal has 
the disadvantage that it involves very small energy scales and long time scales, 
since in the Mott limit (t/U)^ < 1. 

Biichler et al. [33] consider a standard single component Bose-Hubbard model in 
a square (cubic) lattice coupled to lattice of diatomic "molecules" , trapped in the 
centre of each plaquette (see Fig. [29)) The coupling to the molecular state (of d- 
symmetry) takes the form: 

Hm = v ^ "^n"in + 5 ^ TO|!-](6i&3 - &2&4) + /i.e. . (109) 

plaquettes □ plaquettes □ 



Two atoms may perform a Raman transition to a molecular state with the coupling 
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Figure 29. The scheme of realizing ring exchange interactions (redrawn from 1431 ) 

g] V denotes detuning of this transition. Perturbative elimination of the molecules 
leads to the effective ring exchange Hamiltonian for bosons: 

Hre =K ^ {b\b2blbi + h%h\hi ~ ni7i2 - 713^4) . (HO) 

plaquettes □ 

More precisely, the model can be realised using atoms with two internal states: 
one state trapped in the square (cubic) lattice and described by the simple Bose- 
Hubbard model, and a second one trapped in the centres of the plaquettes in a site 
potential that is not symmetric, but rather has a point symmetry of the lattice. 
"Molecules" are build from two atoms in that second internal state, and have the d~ 
symmetry. The model described by the Bose-Hubbard Hamiltonian (Eq.([T4])) plus 
Hm is a promising candidate for a deconfined quantum critical point in 2D [591j . 
This is because it undergoes (most probably) a quantum phase transition between 
the bosonic superfluid phase (which occurs when 1/ 3> g, and t ^ jv^ and which 
breaks the ?7(l)-symmetry), and a molecular "density wave" (stripes) phase (which 
occurs in the opposite limit, when < and \v\ 3> g^ and breaks the translational 
symmetry). Both in 2D and 3D this model describes a C/(l) gauge theory, and is 
likely to exhibit a variety of quantum phases, such as C/(l) deconfined insulator, 
f/(l) deconfined phase, etc. [TOgilM^ . 

Just before this review was submitted, Tewari et al. [138] proposed to use a dipolar 
Bose gas to simulate the compact J7(l) lattice gauge theory. Such model is well 
described by an extended Bose-Hubbard Hamiltonian, as discussed in [147j . Tewari 
et al. propose to use 2D kagome, or 3D pyrochlore lattices, assuring isotropy of 
dipolar interactions between nearest neighbors on the lattice. In the limit in which 
the on-site interactions U are comparable to the n.n. interactions V ^ and both are 
strong, the physics is dominated by the configurations where the fiuctuations of the 
number of atoms in all elementary plaquettes is zero. In the lowest relevant order 
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of perturbation theory {{t/U)^) one obtains a hamiltonian with "ring exchange" 
term on a dual hexagonal lattice. This model reduces directly to the U{1) lattice 
gauge theory, and allows to realize various fractionalized topological phases. More- 
over, the authors propose methods of detecting signatures of the emergent U{1) 
Coulomb phase, terming it as "emergence of artificial light in an optical lattice". 
It is finally worth stressing that by combining, vortex lines with fermions, one can 
realize perhaps the first connections of ultracold atoms with string theory [593j . 



8 Quantum information with ultracold gases 
8.1 Introduction 

There has recently been a real explosion of interest in the studies of the interface 
between quantum information (QI) and many body systems, both in condensed 
matter, as well as in the physics of ultracold atomic gases. First, in Ref. [31] a 
proposal of using ultracold atomic gases in optical lattices for QI was formulated. 
This has stimulated interest in distributed quantum information processing. Second, 
in Refs. |594 , 595, 596J. the first connections between entanglement and quantum 
phase transitions (QPT) were discussed. This discussion has opened fundamental 
questions of deeper understanding of QPT's as well as practical questions of em- 
ploying QI ideas in simulation codes for many body quantum systems (see subsec. 

ran. 

In this section we will (partially) review these two major themes of such studies. 
One of them uses realizable many body systems as viable systems for quantum 
computation and other quantum information tasks. In the other theme, the aim is to 
use quantum information concepts such as entanglement, to understand condensed 
matter and many body phenomena, like for instance quantum phase transitions. 
This section has a slightly different character than the rest of the paper, in that 
it does not make distinction between condensed matter and atomic systems, but 
rather considers general properties of quantum states. It is also organized differ- 
ently. In order to make the review self-contained, we begin with subsec. 18.21 by 
presenting necessary basic notions and tools of quantum information theory in 
general, and entanglement theory in particular. In the next subsection (subsec. 
18. 3p . we consider quantum phase transitions from the quantum information per- 
spective. In subsec. [8.41 we discuss possibilities of performing quantum computation 
using atoms in optical traps. The last subsection 18.51 is devoted to the concept of 
measurement-based quantum computation, both in ordered and disordered lattices. 



8.2 Entanglement: A formal definition and some preliminaries 

The concept of entanglement can be traced back to manuscripts of E. Schrodinger 
around 1932, when it was realised that the quantum mechanical formalism allows 
the existence of pure states of two systems, for which the information about the 
whole system is not the sum of the information about the separate systems. With 
the advent of quantum information, a more operational definition of entanglement 
was required. Such a definition was formulated by Werner in 1989 |597j . Suppose 
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that there are two systems belonging to two observers called Alice (A) and Bob 
(B). They are, in principle, at distant locations, which, operationally means that 
they can only act on their local parts. Under such conditions, the most general 
state that they are able to prepare is of the form 

QA ® Qb- 

However, they may additionally be allowed to communicate classical information 
between them. In that case, the most general state that they are able to prepare is 
of the form 

Y.P^Sa®Qb, (111) 

i 

where the {pi\ are a set of probabilities. The set of operations that Alice and Bob 
are allowed to perform is then called local operations and classical communication 
(LOCC). The states that are of the form in Eq. are called separable states, 

and they (and only they) can be prepared under LOCC only. However, there are 
states that cannot be written in the form in Eq. (jllip . and consequently it is not 
possible to prepare them by LOCC, an example being the singlet state. These states 
are called entangled states. 

A pure state separable if and only if its local subsystems are also pure, 

that is when it can be written as ® |'!/'2)b- The case of mixed states is, 

however, nontrivial. In fact, so far an operational necessary and sufficient criterion 
for detecting entanglement in mixed states has not been found (see |598] for a 
tutorial). Using a semidefinite programming, it has been shown that separability 
can be tested in a finite number of steps, although there is no limit on the number of 
steps needed |599j . See also Ref. [600] . However, there exist many efficient criteria, 
that are either necessary or sufficient. One of them is described below. 



8.2.1 The partial transposition criterion for detecting entanglement. 
The partial transposition criterion is a very useful criterion for detecting entangle- 
ment. We will require the notion of partial transposition of a density matrix qab- If 
partial transposition of qab is taken over the subsystem A, denoted by Qab, then 
its elements are defined as 

{PAB)if^J'^ {PAB)]ii.iuj (112) 

where the subscripts are for Alice's subsystem and the subscripts /i, are for 
Bob's subsystem. Note that although the partial transposition depends upon the 
choice of the basis in which qab is written, the eigenvalues do not depend on 
the basis. We say that a state has positive partial transposition (PPT) whenever 
> 0, i.e. all the eigenvalues of g^^ are nonnegative. Otherwise the state has 
negative partial transposition (NPT). Peres [601] noticed that for separable states, 
the eigenvalues of the partially transposed density matrix are always positive. That 
is, all separable states are PPT, and there exist entangled states which are NPT. 
Failing to be PPT is thus a signature of entanglement. An NPT state is often said 
to violate the PPT criterion. One of the most important results of the theory of 
entanglement has been proven by the Horodecki family, who showed that in 
and (8) C'^, a state is separable if and only if it is PPT [602| . This is no longer 
true in higher dimensions. 
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8.2.2 Entanglement measures. The partial transposition criterion provides 
a tool to check whether the given state is entangled or not. Now, we will discuss 
ways to find out how much entanglement a given state has, once we know that it is 
entangled. This quantification is necessary, at least partly because entanglement is 
viewed as a resource in quantum information theory. There are several complemen- 
tary ways to quantify entanglement (see [60311604060506060607060806090610061111612] 
and references therein). We will present here three possible ways to do so. 

Entanglement of formation. Consider a bipartite pure state \tpAB) shared be- 
tween Alice and Bob. It was shown by Bennett et al. [613j . that given nE[ipAB) 
copies of the singlet state shared between Alice and Bob, they can by LOCC, 
transform them into n copies of the state \i}}ab)^ if n is large, where: 

^(V'ab) = SiQA) = S{qb), (113) 

with QA and qb being the local density matrices of \iIjab)- Here S{p) is the von 
Neumann entropy of p, given by: 

S{p) = -trplog2p. 

If the state |^a_b) describes a C*^ ® C^^ system, then E ranges from to logjd, 
where d = minjiii, ^2}- Note that E is vanishing for separable (pure) states. A state 
is called maximally entangled in C'''^ (8) C'^^, if its E is \0g2d. The singlet state is a 
maximally entangled state in (S) and its entanglement is unity. This amount 
of entanglement of a singlet is also called an ebit. A maximally entangled state in 
C'^'^ (g) C^^ has \0g2d ebits. An example of a maximally entangled state in C'*^ (E) C*^ 
is: 

1 

where \ai)j^ (|/3i)_b) are orthogonal states in Alice's (Bob's) subsystem. With the 
above terminology, we can therefore say that the state \iPab) has E{iPab) ebits. 
Since E{'iPab) is the number of singlets required to prepare a copy of the state 
\iPab), it is called the "entanglement of formation" of IV'As)- We are therefore 
using the amount of entanglement of the singlet state as our unit in quantifying 
entanglement. 

For (bipartite) pure states, the entanglement of formation is essentially the only 
(asymptotic) measure of entanglement, as the singlet state can be locally obtained 
from a pure state at the same asymptotic rate |613j . 

After obtaining the definition of entanglement of formation for pure states, let us 
now define this measure for mixed states j603j . Any mixed state can be expressed 
as a mixture (convex combination) of pure states. That is, a mixed state g can 
always be expressed as: 

i 

where {pi} are a set of probabilities. The states are not necessarily orthogonal. 
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It is then tempting to define the entanglement of formation of a mixed state 

QAB = X^P* Wab) {^ab \ > 

i 

as the average 

i 

However, a mixed state can be decomposed into pure states in an infinite number 
of ways. The entanglement of formation Ep of qab is defined as the minimum over 
all such averages. Precisely, |603j : 

Ep{qab) min^p,£;(V'AB); 

i 

where the minimum is over all decompositions of qab into pure states. 
The procedure for calculating entanglement of formation of any state in ® 
was given in Refs. [61411615] . For states in higher dimensions, the entanglement of 
formation (or even more difficult asymptotic entanglement cost for many copies) 
has been calculated for only a few rare instances using certain symmetries (see for 
instance Ref. [STCllCTTllCTBlingilS^ 1 . 

Concurrence. As stated before, the entanglement of formation of an arbitrary 
(possibly, mixed) state in €? ® has been calculated. It is given by: 

where 

C = C(p) = max{0, Ai — A2 — A3 — A4} , 

and H{-) is the binary entropy function, defined for < x < 1 as H{x) = —x log2 x— 
(1 — x) log2(l — x). The Ai are the eigenvalues (with Ai being the greatest) of the 

Hermitian matrix {^q'^qq'^^ . where q = Oy ® ayQ*ay ® dy, with the complex 
conjugation over q being taken in the eigenbasis. 

In ® C^, the quantity C(gAB), called the concurrence, is defined for arbitrary 
states, and moreover, the entanglement of formation and concurrence are mono- 
tonically nondecreasing functions of each other, with both quantities ranging from 
to 1. Therefore, in (X) C^, we will also use the concurrence as a measure of 
entanglement. 

Logarithmic negativity. The partial transposition criterion can be used to define 
a useful measure of entanglement, that is easily computable. To define the measure, 
let us first introduce a quantity called negativity. The negativity N^qab) of a 
bipartite state qab is defined as the absolute value of the sum of the negative 
eigenvalues of The logarithmic negativity (LN) is defined as [621] 



EnIqab) = \og2{2N{gAB) + 1). 
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Note that for two qubit states, g/'g has at most one negative eigenvalue |622| . 
Moreover, it foUows from our discussion in 18.2.11 that for states in ® and in 
(8) C'^, a positive LN imphes that the state is entangled, while En = implies 
that the state is separable. 

The logarithmic negativity has an operational interpretation in terms of another 
measure related to the entanglement of formation [623| . and moreover, is an entan- 
glement monotone under deterministic LOCC [624j . However, it is not convex. 



8.3 Entanglement and phase transitions 

Over the past few years, the question whether quantum phase transitions (QPT) 
can be understood from the perspective of quantum information has been posed, 
and at least partially positively answered. 

QPTs are nonanalyticities in the ground state energy of physical systems. A typical 
situation is the following (see e.g. [154] ). Consider a physical system that is defined 
on some lattice, and that is described by a Hamiltonian of the form 

H = Ho + XHi, 

where [Hq, Hi] ^ 0. The parameter A is to be visualized as an external parameter, 
which can be changed by the experimenter. In the limit of an infinite lattice, if as the 
system crosses a certain value of A = Ac, the ground state has a qualitative change, 
and there is an associated nonanalyticity in the ground state energy at A = Ac, 
then the system is said to have a QPT at that point. Since the transition is in the 
ground state, so that the system is at zero temperature, the transition is driven 
by quantum fluctuations (and not by thermal fluctuations like in classical phase 
transitions), and this is the reason why it is called a "quantum" phase transition. 

The study of entanglement in strongly correlated systems was initiated in e.g. '626| 
[5?7l.628. 630, 635, 638,.639.,640..642. 641 , 636, 637, 634, 625, 629, 631 , 633, 632, 644, 643] 
(and references therein). It was realized by Osborne and Nielsen |594[|595j and 
Osterloh et al. |596| that the distinctly quantum phenomenon of QPTs can be re- 
lated to the quantum phenomenon of entanglement. The scaling of entanglement as 
an infinite chain of spin-half particles, described by the XY model, goes through a 
QPT was considered. Ideally, one would like to consider the scaling of the multipar- 
ticle entanglement of the whole infinite spin system in the ground state. However, 
the subject of multiparticle entanglement is not fully developed at present. As a 
way out, the authors in Rcfs. [594,595. 1596] considered the nearest neighbor entan- 
glement of the ground state of the spin chain. This required the consideration of 
entanglement of (mixed) states in C^, which is relatively well-developed. We 
consider this topic in subsection 18.3. II 

As we just commented, for higher number of parties, the theory of entanglement 
is not so well developed, even for pure states, at least not quantitatively. However, 
the entanglement of pure states in arbitrary bipartite dimensions is quite well- 
understood, as we have discussed in 18.2.21 Consequently, and since the ground 
state is a pure state, one may consider the scaling of entanglement of the ground 
state for a bipartite split consisting on a certain number of sites n and the rest of 
the lattice, as n grows. The entanglement of one site to the rest of the chain was 
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considered in e.g. Ref. [638|,l639ll595j . and the scaling in this scenario was initiated 
by Vidal et al. [ISB]. We discuss this topic in lSX^ 

Yet another type of scaUng was considered by Verstraete et al. |302| by using the 
concept of "locahzable entanglement" , which will be defined and then its scaling 
discussed in 18.3.31 

We omit several topics in the following, such as the topic of scaling of certain mul- 
tipartite entanglement measures (see e.g. |645[I646[I647[I648[I649| . and references 
therein) . There is also an interesting string of research that deals with criticalities 
in the Hubbard model (see e.g. [650ll6511l646ll652] . and references therein). Dorner 
et al. [65j consider a string of neutral atoms in a one dimensional beam splitter con- 
figuration, where the longitudinal motion is controlled by a moving optical lattice 
potential; they show that it is possible to create maximally entangled states in this 
setup, by crossing a QPT. Lambert et al. [653 considers an infinite collection of 
two-level atoms interacting via the Dickc Hamiltonian with a single bosonic mode, 
and describes the behaviour of entanglement when the system is near its QPT. 

Most of the papers that deal with the scaling of entanglement near criticality deal 
with asymptotic measures of entanglement, or those related to such measures. The 
entanglement of formation (and concurrence), for example, is defined in terms of the 
the entropy of the reduced density matrix of a pure state, which has meaning only 
in the asymptotic regime. Exceptions include the work by Eisert and Cramer [654| . 
which considers single copy entanglement of the ground state of critical quantum 
spin chains. In the last part of this subsection, a possible way to use (bipartite) 
entanglement to observe a type of phase transitions (called dynamical phase tran- 
sitions) in the evolved state of a strongly correlated system is discussed [655] . 



8.3.1 Scaling of entanglement in the reduced density matrix. The first 
system for which the scaling of entanglement near a QPT was considered is a system 
of spin-half particles on a chain, described by the XY Hamiltonian: 

i i 

with J > 0, /i > 0, and < 7 < 1. This system is known to undergo a QPT as the 
parameter X = passes over A = Ac = 1 (see e.g. j6561l658[[657ti659n660ii661ll6621 

HM]). 

In Refs. |5941l595ll596| . it was proposed that this QPT may be analysed by looking 
at the behavior of the entanglement in the two site states (the state obtained by 
tracing over all except two sites of the chain) of the ground state of the system, 
as it goes through the QPT (see Fig. [5D|) . In particular, in the case of the Ising 
Hamiltonian (7 = 1) on an infinite chain, the A-dcrivative of the nearest neighbor 
concurrence diverges as [596]: 

dxCi = -^ln|A- Acl, 



(up to an additive constant) as A approaches Ac (See Fig. [50)1 . The subindex 1 of 
Ci indicates the fact that the nearest neighbor entanglement is calculated. 
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X 

Figure 30. (Color online). Scaling of nearest neighbour entanglement of the ground state of the 
transverse Ising Hamiltonian on a chain (from [663]). The inset at top right shows the nearest 
neighbour entanglement, while the main figure shows its derivative, plotted against the parameter 
A. In the main figure (top right inset), say at A = 1 (say at A = 1.25), the curves from the top are 

respectively for 7, 9, 11, 41 sites (the quality of the figure improves in the online version). 
Bottom left inset: The position Amin of the minimum of d\Ci in the main figure changes with the 
total number of sites. The current inset shows the plot of the logarithm of A„jin — 1 against the 
logarithm of the total number of sites (denoted in this figure as lnf2). 

The Refs. [594115951 1596] stimulated further studies of two-site entanglement to 
characterise QPTs. Bosc and Cliattopadhyay |640] have considered it in the case of 
some frustrated spin models in one dimension, ladder, and two dimensions. Two- 
site entanglement in the XXZ chain was studied by Gu et al. |664] . The Lipkin- 
Meshkov-Glick model has been considered in Refs. [665[l666ll667] . Stauber and 
Guinea [668] have looked at the entanglement of spins at the boundary of a spin 
chain governed by the transverse Ising model, as well as that of two spins coupled to 
a critical reservoir. The anisotropic ferromagnetic Heisenberg chain in the presence 
of domain walls was considered by Alcarez et al. [669| . 

Yang [670] obtained the following apparently counterintuitive result. An exactly 
solvable three-spin interaction Hamiltonian was considered, and it was shown that 
there can appear a discontinuity of the first derivative of the two-site concurrence, 
which has no connection to any QPT (see also Refs. [665[I671| ). It was further 
shown that this discontinuity originates from finding the maximum in the defini- 
tion of concurrence. Subsequently, Wu et al. |672j . found necessary and sufficient 
conditions under which discontinuities in the two-site entanglement (as quantified 
by concurrence, or negativity) can be an indicator of a QPT in physical systems 
governed by Hamiltonians containing only two-body interactions. 



8.3.2 Entanglement entropy: Scaling of spin block entanglement. A 

different approach than the two-site entanglement was considered by Vidal et al. 
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[298j . The idea is to consider a block of spins and consider its entanglement with 
the rest of the system, for the ground state, as the block size increases. Since the 
ground state is a pure state, a good measure of entanglement of a block of spins with 
the rest of the system, is the entropy of the block of spins. This has been termed 
the entanglement entropy (EE) of the model. The models that were considered in 
Ref. |298j (see also |673| ) are the one dimensional XY and XXZ models, where the 
latter is described by the following Hamiltonian: 



It was shown that away from the criticalities, EE saturates to a constant as a func- 
tion of the block size. However, near the criticalities, EE diverges logarithmically 
with block size. Similar behaviour was seen to be true for the AfHeck-Kennedy- 
Lieb-Tasaki valence bond solid models by Fan et al. [674] . 

Techniques from conformal field theory were used by Korepin [650] for evaluating 
the EE in the limit of large block size for infinite spin chains. Its et al. j675| have 
evaluated this limiting EE for an infinite chain governed by the XY model. 

EE was also considered by Wellard and Oriis |676j for the case of QPTs in antifer- 
romagnetic (Ising) planar cubic lattices in the presence of homogeneous transverse 
and longitudinal magnetic fields, and by Latorre et al. Ref. [677| for the Lipkin- 
Meshkov-Glick model. 

The relation between entanglement entropy and area for the cubic harmonic lattice 



has been discussed in Refs. |6781I679| (for fermionic systems, see |680ll681j ) (cf. 
[3T41I682] ). 



In Ref. |683] , an unbounded increase of entanglement entropy was obtained for a 
lattice of spins that interact via some long-range Ising type interaction. 



8.3.3 Localizable entanglement and its scaling. For a given multiparty 
state, £'i,2,...,jv, of N parties, the localizable entanglement (LE) [302] is the max- 
imum average entanglement that can be made to be shared between two pre- 
determined parties (say, 1 and 2), by measurements at the rest of the parties. 
We still have to state which entanglement measure we consider in finding the av- 
erage entanglement. Similar considerations were used in the context of violation of 
Bell inequalities ^EM in Refs. [SHSISHiinST] . In Ref. [3Dg, the LE of the ground 
state of certain spin-half systems were considered, where concurrence was used as 
the measure of entanglement in the definition of average entanglement. 

Using the concept of LE, a related "entanglement length" , ^e, was defined in [301j . 
which gives the typical length scale (of a spin-half chain, e.g.) at which it is possible 
to create maximally entangled states by performing measurements at the rest of 
the parties. Let us make the definition more specific in the case of a spin chain. Let 
Ei^i^n be the LE that can be made to share between the ith and the {i + n)th spins 
of the chain. If Ei^i^n decays exponentially with n, then the entanglement length 
of the chain is defined by: 



Hxxz 
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If the decay of Ei^i^n is non-exponential, then the entanglement length is infinite. 
A similar definition was given by Aharonov [688] . considering a transition from 
quantum to classical physics. The definition of entanglement length parallels that 
of correlation length, for which if the correlation {OiOi+n) between two obscrvables 
Oi and Oi+n decays exponentially, then the correlation length ^ is defined by: 



A diverging correlation length implies a diverging entanglement length. However, 
the converse is not true: for example, the gapped spin system governed by the 
modified Affleck-Kennedy-Lieb-Tasaki model, which has a finite correlation length, 
was proven to have a diverging entanglement length by Verstraete et al. |301j . 



8.3.4 Critical behaviour in the evolved state. Phase transitions are usually 
viewed for ground states, at least when considered from the perspective of quantum 
information theory. In contrast, Ref. |655| considers a criticality of nearest-neighbor 
entanglement in the evolution of an infinite spin chain described by the asymmetric 
XY model in a time-dependent transverse field. 

The initial state of the evolution is taken to be the canonical equilibrium state at 
zero temperature (which, in this case, is the ground state). The transverse field is 
then suddenly turn off at time equal to zero. The system is thus initially disturbed, 
and its properties are then studied at later times. The nearest-neighbor entangle- 
ment (quantified by the logarithmic negativity) in the evolved state at a fixed time 
shows a criticality (which has been called a dynamical phase transition (DPT)) 
with respect to the transverse external field. The region of the initial transverse 
field for which the entanglement is nonvanishing (vanishing), at a fixed time, is 
referred to as the "entangled phase" ("separable phase") (see Fig. [?T|) . 

Interestingly, the nature of the DPT depends on whether we are near, or far from 
the time of initial disturbance. Moreover, for values of the initial transverse field 
near the criticalities, as well as in the separable phase, and for short times, the 
nearest-neighbor entanglement shows nonmonotonicity with respect to tempera- 
ture. Accordingly, the criticalities are referred to as "critical regions" , signifying 
that "critical" effects persist for a small region around the critical value of the 
transverse field. 

The study of nonmonotonicity of entanglement with respect to temperature is 
interesting, as preservation of entanglement in a hostile environment is one of the 
main challenges in quantum computation and quantum information in general. A 
common belief is that temperature is a form of noise, i.e. it destroys the subtle 
quantum correlations. However, this seems not to be universally true. 

Nonmonotonicity of entanglement with respect to temperature was also found by 
Arnesen et al. |632| . in the canonical equilibrium state of the one-dimensional 
isotropic antiferromagnetic Heisenberg model. Similar behavior was obtained 
by Scheel et al. j689| . in the Jaynes-Cummings model. The idea that dissipa- 
tion can assist generation of entanglement has also been put forward in Refs. 



It should be mentioned here that the asymptotics of temperature correlations in 




[MilMIlIMaiMaEMlIMlEM] . 
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Figure 31. Dynamical phase transitions: a "river" of separable states between the two regions of 
control parameter space where nearest neighbour entaglement is non-zero in an asymmetric XY 
model in a time-dependent transverse field, (from 16551 ). Entanglement (logarithmic negativity) is 
plotted against the initial transverse field (a) and (real) time (t). A similar behaviour is absent in 

magnetization 16551 . 

the isotropic (i.e 7 = 0) XY model was considered by Its et al. [697j (see also 
[6981l6991l7001l701| ) . Also, previous studies of the quantum dynamics of spin models 
after a rapid change of the field include Refs. [6g7l[6Ml[66m[66Tl[66l [3l [70l[70^ . 
while effects of a sudden switching of the interaction in arrays of oscillators were 
studied in Ref. [TDi] , 



8.4 Quantum computing with lattice gases 

The computers that have become a part of our daily life are very efficient and useful 
machines. However, there arc certain problems that may not be efficiently solved 
with such computers. For example, it is believed that the problem of finding prime 
factors of a given integer is exponentially inefficient on a classical computer, since 
the number of steps required to find the factors scales exponentially with the size 
of the input integer. 

Classical computers follow the principles of classical physics, but it is possible to 
consider computers that are ruled by the laws of quantum mechanics. There are 
at least two important motivations for considering computing at a quantum level: 
(1) The classical laws of physics governing the functioning of classical computers 
are an approximation that suits our purpose for the current size of microchips. If 
size is going to keep decreasing, so that computational speed increases, quantum 
effects will unavoidably emerge; (2) In 1994, Shor discovered a quantum algorithm 
that can efficiently factorize integers into their prime factors |23| . 
There are several proposals for implementing quantum computing [705j . A very 
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Figure 32. Schematic representation of the gate operation. Solid and dashed Unes represent 
potentials felt by the atoms in the different internal states, respectively. By shifting the potentials, 
the pair of atoms in the middle experiences a coUisional shift (from 1311 ) 



promising one was put forward by Cirac and ZoUcr in 1995 that uses cold trapped 
ions interacting with laser beams [53] . For the effects of environmental decoherence 
on such an ion trap quantum computer and possible ways to get over them, see 
e.g. [706l[707l[708ll709[l710[|7TT] . and references therein. Following the Cirac-ZoUer 
proposal, a fundamental quantum logic gate, the CNOT gate, was experimentally 
demonstrated by Monroe et al. [712] , where the two qubit state was encoded in the 
two internal states and two external states of a single trapped ion. For recent devel- 
opments in ion trap computers, see |7131I7141|7151I7161I717| . In particular, Ref. [714j 
has reahzed the CNOT gate by two ''°Ca+ ions, held in a linear Paul trap. Jaksch 
et al. [31j proposed to use atoms in an optical lattice for quantum computation. 
They use cold controlled collisions between two atoms in moving trap potentials 
for implementing the following two-qubit gate: 

|a)i|a>2 -> |a)i|a)2, 

|a)i|fo)2^e-^^-|a>i|^)2, 
|6)i|a)2^ |&)i|a)2, 

\b)l\b)2^\b)i\b)2. (114) 

The two qubits are encoded in the internal states \a) and |6) of the two atoms. 
The two atoms are labelled as 1 and 2 and each atom is trapped in the ground 
state of a potential well. The phase (j^ab will be a fixed number depending on the 
implementation. Initially, the potential wells are at a distance that is sufficiently 
large distance so that there is no interaction between the atoms. They are then 
moved along certain trajectories toward each other, so that the atoms interact for 
a certain time, and then finally the wells are restored to their initial positions. The 
potential wells are moved adiabatically, so that the atoms remain in the respective 
ground states of the wells. The implementation requires an interaction that induces 
internal state dependent potentials, so that the atoms in different internal states 
feel different potentials (see Fig. [5^ . 
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The effective Hamiltonian governing the system is given by: 

i 

ij 

Qi (a|) and bi are the annihilation (creation) operators of the internal levels 
I a) and \b) for an atom in the ground state of the potential well at position i. By 
assuming that (i) there is no sloshing motion of the atoms in the wells when the 
latter are moved, (ii) the potentials are moved so that the atoms remain in the 
ground state (adiabaticity condition), and (iii) the rms velocity of the atoms in the 
vibrational ground state is sufficiently small so that zero energy s-wave scattering 
approximation is valid, then it is possible to show that the internal states of the 
atoms transform in the way described in Eq. (|114p . 

There are several other papers that deal with quantum computation with atoms 
in optical lattices. In particular, Brennen et al. ^718j have used dipole-dipolc inter- 
actions to implement the CNOT gate. For further work, see e.g. 1730,741,32,732', 

\m\mmmmmmi\m ■ 



8.5 Generation of entanglement: The one-way quantum computer 

There has been several studies for generating entanglement in different physical 
systems. We want to discuss generation of entanglement from the perspective of 
performing quantum computation. We shall discuss entanglement generation in 
ordered and disordered lattices separately. 



8.5.1 The one-way quantum computer. The one-way quantum computer 
(also referred to as "measurement-based quantum computation" ) was proposed by 
Raussendorf and Briegel |743j . who showed that arbitrary quantum computation 
can be simulated via single-particle measurements on a specially prepared quantum 
state termed as "cluster state". Cluster states can be created efficiently in any 
ultracold atomic systems with Ising-type interactions, and allow the possibility of 
preparing the system in a specific initial product state. A proposal to implement 
cluster states uses four atomic levels trapped in appropriately tuned optical lattices 

[7I3j (cf. UMEM)- 

To define cluster states |747] , consider an arbitrary graph made up of N "vertices" 
(which will be identified with the qubits making up the quantum computer), and a 
certain number of "edges" (which will be identified with Ising interactions between 
the qubits) connecting the vertices (see Fig. [33]). The terminology used here is 
from graph theory (see e.g. [748] ) . and thereby cluster states are also called "graph 
states" . Suppose that the qubits are spin-half particles, and that each of them is 
prepared in the state \+) = (|0) + |1))/a/2, where |0) and |1) are respectively the 
spin up and down states in the z direction. We then apply an Ising interaction 
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Figure 33. An example of a graph. In this case, it has five vertices and four edges. The vertices 
represent qubits, and the edges represent Ising interactions between them. 



governed by the Hamiltonian 

between the qubits, where the summation runs over aU pairs of sites that are 
connected by an edge in the graph. After a certain finite time (= fj), the state 

|0<,)=J]C^,^.|+)«^ (115) 

is created. Here CSij is the controUcd-PHASE gate, defined as 

|0)|0)-|0)|0), 
|0)|1)-|0)|1), 
|1>|0)-|1>|0>, 

ll>ll>-*ll>ll>, 

acting on the spins i and j. The product on the right hand side of Eq. pisp . 
again runs over all pairs (i, j) of sites that arc connected by an edge in the graph. 
The state \4>c) is called the cluster state corresponding to the given graph. These 
states have several interesting properties. Perhaps the most interesting property is 
that for sufficiently large graphs, the corresponding cluster states can be used as 
a substrate for quantum computational tasks. More specifically, any quantum gate 
can be implemented with this resource, simply by performing measurements on 
the different spins of a suitably chosen graph, where it is assumed that the choice 
of a measurement basis on a particular spin may depend on outcomes of previous 
measurements on other spins. As an example, the Hadamard gate (H), defined as 

|0)^(|0) + |1))/V2, 

|1) - (|0) - |l))/^/2, 

can be implemented in the following way. Suppose that at a certain point of com- 
putation, we obtain the state \'ipin)i which we want to transform into Hj^/^in). 
This can be realized by using a graph of five vertices in a chain as shown in Fig. 
[34l j743p749[750j . The input state is initially at the extreme left. The rest of the spins 
are initially in the state |+), C2I+) = |+). Ising interaction with nearest neighbor 
interactions (as depicted in the figure by the edges) for the requisite amount of 
time (equivalent to the controUcd-PHASE gates, as described in the text) results 
in a (entangled) state of the five spins. Subsequently, measurements are performed 
on all spins except that on the extreme right, in the eigenbasis of the operator 
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indicated over the corresponding vertices. This resuhs in the state Hj?/';,!) at the 
extreme right spin. A five-spin graph and corresponding measurements implement- 



(7,, 



+ 



+ 



+ 



Figure 34. A graph to implement the Hadamard gate. 

ing an arbitrary singie-qubit gate, and a 15-spin graph implementing the CNOT 
gate is given in the same papers. 

There is a large body of work in this direction that deal with several important 
issues ( see e.g. |751l752ll753l754y755l756y757l758l75Hy76()l7H4ll76TT762y765l7f)6y7()8l 
[767l[7551[7M[770l[77Tl[77l[77ll774ll^ ) . In particular, Diir et al. [Hp deals with 
the entanglement properties of the time-evolved states in a model with long-range 
Ising-type interactions, in the thermodynamic limit. 



8.5.2 Disordered lattice. Interestingly, it is also possible to perform quantum 
computing in systems with quenched disorder. This may sound contradictory, but 
we want to employ here the possibility of creating controlled disorder in atomic 
gases in optical lattices. Here, by "controlled" , we mean that a particular realization 
of disorder remains fixed for times much longer than over which we monitor the 
evolution of the system (the disorder is "quenched"). This allows, as discussed in 
Section m to study Anderson and Bose glasses in a Bose gas [B^, or spin glasses 
with short range interactions in Fermi-Bose, or Bose-Bose mixtures [51]. Using 
linear chains of trapped ions |495| , or dipolar atomic gases }144ll775j , it is possible 
to realize complex spin systems with long range interactions that may serve as a 
model for classical and quantum neural networks [776j (cf. [777] ). 
Disordered systems offer at least two possible advantages for quantum computing. 
First, since they have typically a large number of different metastable (free) energy 
minima, such states might be used to store information distributed over the whole 
system, similarly to classical neural network (NN) models [778] . The information is 
thus naturally stored in a redundant way, like in error correcting schemes |779II780| . 
Second, in disordered systems with long range interactions, the stored information 
is robust: metastable states have quite large basins of attraction in the thermody- 
namical sense. 

The aim of using complex disordered systems is therefore to be able to perform 
distributed quantum computing in a way that is resistant to noise. In Refs. |781[ 
1782] . the following questions in the nan- distributed case were addressed: Q (i) Can 
one generate entanglement in such systems that would survive quenched averaging 
over long times? (ii) Can one realize quantum gates with reasonable fidelity? Both 



^ Since the submission of an earlier version of this review, the distributed case has also been considered 
[776ir783] . 
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questions were answered affirmatively for botli short and long range disordered 
systems. Here we present as an example, the second question for the Hadamard 
gate [78T] . 

We assume that the computation is performed in a spin lattice created with cold 
ions in an optical lattice, of which particles 1 and 2 are part of. At a certain time, 
particle 1 is in an arbitrary state ji/'iTi), while 2 is in an eigenstate of cr^ with 
eigenvalue +1 e.g. |+). Then we let 1 and 2 evolve according to the Edwards- 
Anderson Hamiltonian 

He- A ~ —J(Tla2 

for a suitable duration of time. Here J is a random quenched coupling, that in the 
case of the Edwards- Anderson model, can be regarded a Gaussian random variable. 
After a suitable duration of evolution, a measurement is performed on particle 1 
(in a suitable basis). For the case when J is a Gaussian variable with mean= 5 
and variance=l, the state in particle 2 attains the Hadamard rotated state H \'ipin), 
with quenched averaged fidelity greater than 0.85. One can increase such fidelity by 
increasing the number of spins, and employing measurements assisted by outcomes 
of previously performed measurements. Note, that the fidelity of the Hadamard 
rotated state using the classical information obtained only from the measurement 
of particle 1, is only 2/3 [7M] . 



9 Summary 

Summarizing, we have tried to survey the state-of-the-art on the active research 
field of ultracold atomic gases in optical lattices. On one hand, we have aimed at 
giving a general outline of new topics and open problems that can be addressed with 
these systems. Our list, however, is by no means exhaustive. On the other hand, 
we have briefiy reviewed some of the best established methods and techniques to 
study strongly correlated systems in the framework of ultracold gases. In spite of 
the fact that we have mostly focused on the connections of ultracold gases with 
condensed matter topics, we envision and hope that, in the next years, the frontiers 
of this field will greatly expand. 

Finally, we hope we have at least partially convey our enthusiasm for ultracold 
gases physics. A concise version of our summary is contained in the motto, which 
has in fact been used by one of us (M.L.) at the end of his conference presentations 
and seminars. This motto should express a pure joy of curiosity and discovery (even 
if it was rc-discovery) that we and many others have experienced while working on 
ultracold atomic gases over the last years. This is the most important message: 
Ultracold atom physics is FA-A-ANTASTIC!!! 
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Appendix A: Effective Hamiltonian to second order 

The Bose Hubbard Hamiltonian or Fermi-Bose Hamiltonian can be split into 
H = Hq + Hint, where Hq denotes the on-site Hamiltonian obtained in the limit of 
zero tunneling, and the hopping term is denoted by Hint- The Hamiltonian Hq can 
be easily diagonalised and has well defined eigenstates, that are grouped in blocks 
(or manifolds). Each manifold is well separated from each other. Typically differ- 
ent manifolds are separated by terms of the order of U, being U the interaction 
between two bosons. We denote by Pa, the projector on each block space, where a 
is the block index, and the z-th state in any block is denoted by \a,i). Note that 
PaHoPfs = holds for a ^ f3. The Hint part of the Hamiltonian will introduce cou- 
plings between the bare block a and /3, i.e., PaHintPp can be different from zero for 
a ^ p. Following [230] , one can construct an effective Hamiltonian, _ffe/ / , from H 
such that it describes the slow, low-energy perturbation-induced tunneling strictly 
within each manifold of the unperturbed block states and has the same eigenval- 
ues as H . Tunneling processes between different manifolds are thus neglected. We 
demand that the effective hamiltonian i?e// fulfills that: 

(i) is by construction hcrmitian, with the same eigenvalues and the same degen- 
eracies as i/. To achieve that, one defines T := e*'^, with S hermitian, S = S\ 
and chosen such that: 



(iii) As the first two conditions still allow for an infinite number of unitary transfor- 
mations (all UT are still possible, U being any unitary transformation acting 
only within the manifolds), the following additional condition is imposed: 



H,ff = THT^ . 
(ii) Heff does not couple states from different manifolds: 



(Al) 



PaHeffP0 = O, a ^(3 



(A2) 



PaSPa — for any a. 



(A3) 
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Expanding the first condition using the Baker-HausdorfF formula, one obtains: 
H,ff^H+[iS,H] + ^[tS, [iS,H]] 

+ ^[iS,[tS,[iS,H]]] + ... . (A4) 

Making a power-series ansatz for in the perturbative parameter t S, 

S^tSi+ t^S2 + t^S3 + ..., (A5) 
and employing H = Hq + tHint one obtains from (|A4p to second order 



-^e// — Hq + t{[iSi,Ho] + Hint) 



iff 



( [iS2, Ho] + [tSi,HM] + \ [iSi, [iSi,Ho]] ] . (A6) 

This is a power series for H^ff, with its moments denoted by H^^j, H^j^.. to first 
order, i.e. -ffe// = Hq + tH^ij and S = tSi. Using the second and third conditions, 
as well as PaHoPp=0 and the expression for -ff^// (|A6p . one finds: 

{a,i\iSi\p,j) {Epj - E^i) + {a,i\H,nt\P,j) = (A7) 

^{a,l\^S^\P,]) = \ e^.-e,, . (A8) 

a — p 

Thus, the effective Hamiltonian within the a -manifold, depends only on Hint and 
not on 5*1, i.e., {a,i\Hlj j\a, j) = {a,i\Hint\a, j) . A general result for any n is that 
{a,i\H^j j\a, j) is independent of 5*,!. Based on the third condition, and on the 
observation that Sn enters the expression for only in the commutator with 

Hq, which is diagonal in the manifold index. 

Thus, when continuing to second order, the term [iiS'2,i/o] hi the expression for 
H^fj can be dropped. Of the two remaining terms defining H^^^ in (|A6|) . the 
second one can be simplified by observing, that according to (jA7p the operator 
[iSi,Ho\ is purely non-diagonal in the manifold index, with values opposite to 
those of the non-diagonal part of Hint- Thus, \ [iSi, [iS'i,i?o]] ~ —\ ^^i>i't] • 
Now inserting the identity between the operators in the still untreated second term 
in H^ff, [iSi,Hint] ■, one sees that due to 5*1 being non-diagonal in a, again only 

the non-diagonal part of Hint can contribute: [iSi,H,nt] = [iSi,H^^^]. Therefore, 
one has: 

Hi J f = [zS^, H^ni] + l[tS^, [is^Hol] = ^ , ifr„1] . ( A9) 
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Collecting all terms relevant for (a, i\Hef / |a, j) to second order in t, and introducing 
the notation Q^i := J2k,-f=^a g^fJ^'^'j , one finds: 

{a,i\Heff\a,j) = E^iSij + t{a,i\Hint\a,j) 

-— {{a,i\H,nt [Qai + Qaj] H,nt\a, j)) , (AlO) 

where the identity operator has been inserted in the final expression for H^j.j.- in 
formula (jA9[) . and then evaluated using formula (jASp . which naturally leads one to 
define the operator Qai as above. Note that this construction can be generalized 
to arbitrary orders in t in a straightforward manner, as detailed in [230j . 



Appendix B: Size of the occupation-reduced Hilbert space 

We derive here Eq. ((26|) . i.e., an expression for the size of the Hilbert space in the 
system subjected to the restriction that at most K bosons can populate a single lat- 
tice site. We assume that there are N bosons placed in M lattice sites, therefore K < 
N and KM > N. By definition H5(7V,M)|^ = Ef,=o ' ' ' Ef„=o '5"i+-+nM,Ar. 
Using the following Kronecker delta representation, Si^m = ^ d(^e*'^'~™''^, one 
gets 

where the right hand side is obtained after integration on a unit circle z — cxp^iip) 
in the complex plane. Then, using the Leibnitz formula one obtains 

which can be further reduced to ([26]). 
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